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Liquid crystal elastomers realize a fascinating new form of soft matter that is a composite of 
a conventional crosslinked polymer gel (rubber) and a liquid crystal. These solid liquid crystal 
amalgams, quite similarly to their (conventional, fluid) liquid crystal counterparts, can sponta- 
neously partially break translational and/or orientational symmetries, accompanied by novel soft 
Goldstone modes. As a consequence, these materials can exhibit unconventional elasticity charac- 
terized by symmetry-enforced vanishing of some elastic moduli. Thus, a proper description of such 
solids requires an essential modification of the classical elasticity theory. In this work, we develop 
a rotationally invariant, nonlinear theory of elasticity for the nematic phase of ideal liquid crystal 
elastomers. We show that it is characterized by soft modes, corresponding to a combination of 
long wavelength shear deformations of the solid network and rotations of the nematic director field. 
We study thermal fluctuations of these soft modes in the presence of network heterogeneities and 
show that they lead to a large variety of anomalous elastic properties, such as singular length-scale 
dependent shear elastic moduli, a divergent elastic constant for splay distortion of the nematic di- 
rector, long-scale incompressibility, universal Poisson ratios and a nonlinear stress-strain relation for 
arbitrary small strains. These long-scale elastic properties are universal, controlled by a nontrivial 
zero-temperature fixed point and constitute a qualitative breakdown of the classical elasticity the- 
ory in nematic elastomers. Thus, nematic elastomers realize a stable "critical phase", characterized 
by universal power-law correlations, akin to a critical point of a continuous phase transition, but 
extending over an entire phase. 



I. INTRODUCTION 
A. Classical and generalized elasticity 

The classical theory of elasticity [J, |3| is a well- 
established discipUne developed in the mid-nineteen cen- 
tury by Hooke, Cauchy, Poisson, Green and many others. 
It models conventional solid materials as continuous me- 
dia with an energetically preferred, stress-free, reference 
configuration. In this theory, the elastic energy of a solid 
is a functional of the deformation tensor, defined rela- 
tive to the aforementioned stress-free reference state. As 
shown by Cauchy, the conservation of angular momen- 
tum requires the stress tensor in conventional solids to be 
symmetric, encoding the underlying physics that torque 
cannot be transmitted through the media without the 
transmission of force. 

The development of solid-state physics provided a mi- 
croscopic understanding of inter-atomic cohesive forces 
inside crystalline solids, and allowed a first-principle 
derivation of their elasticity, as well as the computation 
of the associated elastic constants. 

In another major development, complementing this mi- 
croscopic description, it was appreciated that the rigidity 
of crystalline solids, which distinguishes them from liq- 
uids, arises as a consequence of spontaneous breaking of 
translational symmetry. Transcending microscopies, the 
functional form of the associated elastic energy is com- 
pletely determined by the symmetry, i.e., the space sym- 
metry group, of the resulting ordered crystalline state. 



From this perspective, the spontaneous nature of crys- 
talline ordering requires the existence of the low-energy, 
long-wavelength excitations, i.e., phonons, that are the 
Goldstone modes Q associated with the spontaneous 
breaking of translational symmetry Q . 



This intimate relation between symmetry breaking, 
generalized rigidity and elasticity is one of the most 
important cornerstones of modern condensed matter 
physics 0]. It extends to a wide class of system with 
translational, orientational, and other more exotic or- 
ders, such as magnets, superfluids, and liquid crystals 
. Each system is characterized by its own set of Gold- 
stone modes and generalized rigidity, corresponding to 
the symmetry spontaneously broken by the particular or- 
dered state. For example, a nematic liquid crystal breaks 
the rotational, but not translational, symmetry, and, as 
a result, is able to transmit a torque, but not a force. It 
thereby admits an asymmetric stress tensor, i.e., contains 
both a symmetric and an antisymmetric part, in strong 
contrast to an ordinary isotropic liquid. A smectic liquid 
crystal, on the other hand, breaks translational symme- 
try in one direction, along which it is able to transmit 
force. The term (generalized) elasticity, therefore gains 
much wider meaning in the context of modern condensed 
matter physics, applying to all thermodynamic phases 
that exhibit a spontaneous long-range order. 
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B. Classical rubber elasticity 

Despite the tremendous success of this symmetry 
breaking description there are classes of sohds whose elas- 
ticity do not fit into this paradigm. One example is a 
structural glass, which, despite of not breaking any spa- 
tial symmetries, exhibits a macroscopic isotropic solid- 
like elasticity on experimental time scales [6\. 

Elastomers (rubber) constitute another class of ma- 
terials that break no spatial symmetries but neverthe- 
less, macroscopically are true three-dimensional solids, 
in the sense that they do not flow and have a finite shear 
modulus. Microscopically, elastomers are polymer net- 
works formed by crosslinking polymer melts. Strain de- 
formations decrease the entropy of constituent fluctuat- 
ing polymer chains, and therefore cost elastic free energy. 
The shear modulus of elastomers is proportional to tem- 
perature 0] , and is typically 4 — 5 orders of magnitude 
smaller than that of crystalline solids. Due to the flexibil- 
ity of polymer chains, elastomers can support extremely 
large (typically several hundred percent) elastic shear de- 
formations. 

The elasticity of isotropic rubber is usually understood 
within the framework of classical theory of rubber elas- 
ticity [tI, developed long time ago by Kuhn, Wall, Flory, 
Treloar and others. This theory treats, in an obviously 
erroneous way, the crosslinks of polymer networks as sta- 
tionary objects that nevertheless deform affinely with the 
network when it is subjected to an external strain defor- 
mation. The total entropy of the network is then simply 
the sum of those of individual chains. This theory ex- 
plains the stress-strain relation of rubber for small strain 
reasonably well but fails at larger deformations. A num- 
ber of extensions of the classical rubber theory 0] have 
been developed subsequently, which attempt to include 
non-Gaussian chain statistics, irreversibility, short-scale 
nematic and crystalline ordering, and in particular poly- 
mer entanglements [1, [l^l- The success of these ef- 
forts however are rather difficult to assess, largely due 
to the fact that they do not make prediction other than 
stress-strain relations. Thus, despite their technological 
importance and abundant existence, our understanding 
of elastomers remains rather incomplete, and in part is 
due to the lack of proper description of their random 
structures 

A recent exciting development is the observation and 
the subsequent resolution [12] of internal inconsistency 
of the classical rubber theory. Namely, that the entropy 
of fluctuating crosslinks, ignored by the classical theory, 
is comparable to that of the polymer chains. Taking 
into account phonon fluctuations on scales longer than 
the network mesh size, while incorporating rubber's near 
incompressibility, one finds quantitative agreement with 
the stress-strain relation of rubber in the large deforma- 
tion regime jl2| . without appealing to any other more 
complicated and poorly characterized mechanism. 



C. Nematic liquid-crystalline elastomers 

Another fascinating, relatively new class of elastic ma- 
terials that falls outside the scope of classical crystal elas- 
ticity is that of fiquid crystaUine elastomers [sLfia. [T^. [l5| . 
These are crosslinked networks of liquid crystalline poly- 




ABC 

FIG. 1: An illustration of (A) main-chain liquid crystalline 
polymer and (B), (C) side-chain liquid crystalline polymers. 

mers that incorporate mesogenic groups, i.e., chemical 
groups with a sufficient anisotropy to exhibit liquid crys- 
talline orders. These mesogens may either be parts of the 
polymer backbone, separated by flexible chemical spacer 
groups, as in the case of main-chain liquid crystalline 
polymers, or be attached to the main polymer backbone, 
as in the case of side-chain liquid crystalline polymers, 
illustrated in Fig. [TJ Quite similarly to their liquid crys- 
tal counterparts, liquid crystalline elastomers may also 
spontaneously partially break translational and/or rota- 
tional symmetry, when the anisotropic interaction be- 
tween mesogenic groups is sufficiently strong. The re- 
sulting phases have long-range liquid crystalline order in 
the background of a random, statistically isotropic poly- 
mer network. The coupling between network elasticity 
and liquid crystalline order leads to a rich phenomenol- 
ogy @i [l3| and an exciting possibility of controlling one 
using the other. 

Nematic elastomers [isl . are nematically-ordered 
rubbery materials, prepared by crosslinking nematic 
polymer melts or solutions under appropriate conditions. 
In an idealized approach, one would crosslink a homo- 
geneous polymer melt in the isotropic phase, and subse- 
quently lower the temperature so that the nematic order 
appears spontaneously. We shall refer to such a fictitious 
nematic elastomer as ideal. In reality, however, network 
heterogeneity interferes with the nematic ordering and 
a system prepared this way only exhibits a finite-range 
nematic order at low temperature. The precise mecha- 
nism leading to suppression of long-range nematic order 
is, however, not yet completely understood [l6|. To ob- 
tain a single-domain nematic elastomer, the crosslinking 
is instead done slightly in the nematic phase or under a 
weak strain 17]. Hence a weak nematic order is perma- 
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nently imprinted into the system, leading to a small free 
energy cost for the Goldstone modes of the ideal systems, 
a property referred to as semi-softness [31 ■ In this work, 
however, we shall only consider ideal homogeneous nc- 
matic elastomers crosslinked under isotropic conditions. 
We shall consider the effects of network heterogeneity to- 
ward the end of this paper. 




Isotropic Nematic 

FIG. 2: A cartoon of a liquid crystalline polymer coil under- 
going the isotropic-nematic transition, driven by the nematic 
ordering of the mesogenic groups (not shown in the figure) . In 
the isotropic phase the mesogens are randomly oriented in the 
space, while in the nematic phase they are align on average. 
In an elastomer, this is accompanied by a uniaxial distortion 
of polymer coils that in turn leads to a spontaneous uniaxial 
deformation of the underlying polymer network, i.e., the solid 
itself. 

The interplay between liquid crystalline order and net- 
work elasticity becomes most interesting when the liquid 
crystalline order appears spontaneously, and the Gold- 
stone theorem dictates that there must exist correspond- 
ing soft modes whose elastic energy vanishes in the long 
wavelength limit. This fascinating possibility has been 
explored in most detail in nematic elastomers. Due to 
the coupling between the mesogenic groups and the poly- 
mer backbones, the isotropic-nematic (I-N) transition in 
a liquid crystalline polymer melt is accompanied by a 
spontaneous uniaxial deformation of all polymer coils, il- 
lustrated in Fig. II Cl in a crosslinked network of polymers, 
such a shape change of all polymer coils implies the same 
shape change of the whole network, i.e. a spontaneous 
uniaxial stretch of the elastomer [13, [3 • Since the ro- 
tational symmetry is broken spontaneously, there exists 
a continuous manifold of degenerate ground states inside 
the nematic phase, corresponding to distinct choices of 
the uniaxial axis in the statistically isotropic solid. It is 
important to note that these degenerate states are not 
simply related to each other by global rotations of the 
sample. Only a special combination of shear deformation 
of the solid and rotation of nematic director can take the 
system from one such ground state to another. These 
soft modes, first predicted by Golubovic and Lubensky 
[20| based purely on symmetry consideration, are respon- 
sible for most of the unusual properties of nematic elas- 
tomers. As noted by these authors, the effective elastic 
free energy of the nematic elastomer resembles that of a 



uniaxial solid, but with a strict vanishing of one of the 
five elastic constants, enforced by the spontaneous break- 
ing of rotational symmetry and encoding the presence of 
the corresponding Goldstone mode. This then leads to a 
striking elastic response, where under a shear transverse 
to the nematic axis, the stress-strain relation of an ideal 
nematic elastomer exhibits an extended plateau with a 
vanishing stress [21, 22. l23j. 

Complementary to this pioneering work, a very suc- 
cessful molecular level elasticity theory of ideal nematic 
elastomers, referred to as the neo-classical theory, was 
developed by Warner and Terentjev [l^, Essentially 
a spontaneously anisotropic generalization of the classi- 
cal theory of rubber elasticity , it successfully captures 
the soft modes of ideal nematic elastomers. Character- 
ized by only two parameters, the shear modulus /i and 
the nemato-elastic coupling constant, which determines 
the spontaneous distortion ratio, the neo-classical theory 
provides a minimal model remarkably useful for practical 
calculations [l3[. It is important to note that this elastic 
energy cannot be expressed in terms of any symmetric 
strain tensor, linear or nonlinear. This reflects the facts 
that nematic elastomers are not ordinary solids as un- 
derstood by Cauchy and his contemporaries, and that a 
correct elastic description of these materials necessarily 
leads to physics beyond the framework of classical elas- 
ticity theory. 

It is worthy to note that theoretical attempts to go 
beyond the classical elasticity theory date back to e arly 
nineteenth century, by Voigt and Cosserat brothers [24| . 
These authors considered fictitious elastic media, i.e., 
the so-called Cosserat medium, that consists of three- 
dimensional rigid constituents with both translational 
and rotational degree of freedom, quite analogous to the 
modern realization in nematic elastomers. It was under- 
stood quite early that these media can transmit torque 
independent of force, and therefore admit asymmetric 
stress tensors. Receiving little attention for a long time, 
research in Cosserat elasticity was revived about half a 
century later in the engineering community under the 
new name of micropolar media. However, almost all of 
theses studies focus on isotropic micropolar elastic me- 
dia, lacking the strain-rotation coupling. As discussed 
above, this coupling is the most interesting ingredient of 
the elasticity of micro-polar media and underlies all of 
their unusual properties. 



D. Thermal fluctuations, heterogeneity and 
anomalous elasticity: critical phases 

In spite of its remarkable success [l^, the principle 
shortcomings of the neo-classical theory are evident. As 
already discussed above, being built on the classical the- 
ory of rubber elasticity it misses a substantial entropic 
contribution associated with fluctuations of the network 
crosslinks [l^]. Neither does it address network hetero- 
geneity, i.e., elastomers quenched random structure. As 
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have been recently demonstrated in a series of publica- 
tions [2^, [2^, [13, [2^ , thermal fluctuations and network 
random heterogeneity have even more drastic effects on 
the macroscopic elasticity of nematic elastomers, as com- 
pared with the conventional isotropic rubber. This is 
mainly due to the aforementioned existence of soft Gold- 
stone modes, together with the resul ting strict vanish- 
ing of one of the shear moduli [l^, Hif. As we will 
show in detail below, the presence of local fluctuations 
of these soft modes, relevant elastic nonlinearities lead 
to a qualitative modification of the macroscopic elasticity 
of nematic elastomers. As a result, at long scales these 
materials exhibit a variety of anomalous elastic proper- 
ties. Namely, all shear moduli, as well as the splay con- 
stant for nematic distortion, become singular functions of 
the probing length-scale and macroscopic strain deforma- 
tion, and either diverge or vanish in the thermodynamic 
limit. Related to these, the anomalous solid exhibits a 
long-scale incompressibility, universal Poisson ratios, as 
well as nonlinear stress-strain relation for arbitrary small 
strains. These properties are universal, controlled by a 
non-Gaussian, finite-temperature (for idealized homoge- 
neous elastomers) fixed point, as illustrated in the verti- 
cal temperature axis in Fig. [3l and constitute a qualita- 
tive breakdown of classical elasticity theory. Similar phe- 
nomena have been previously shown to appear in other 
"soft" solids, such as smectic and columnar liq uid crys- 
tals [13, [U, [13] , polymerized membranes [33l [11] and 
putative spontaneous vortex lattices in ferromagnetic su- 
perconductors 35, 36]. 

These effects are similar but qualitatively stronger in 
(more realistic) elastomers with random heterogeneity, 
controlled by a zero-temperature, finite-disorder fixed 
point, as illustrated in the horizontal axis of Fig. [3l 
It characterizes the "rough" ground state of a hetero- 
geneous elastomer that exhibits large sample-to-sample 
fluctuations relative to an idealized homogeneous elas- 
tomer MM. 



The physical mechanism of anomalous elasticity in ne- 
matic elastomers can be visualized by considering the 
low-energy elastic excitations, illustrated in Fig. 01 A 
rotation of the nematic director and a simultaneous uni- 
form shear deformation, shown in Fig. |4(3, is a soft mode 
and costs zero elastic free energy (37,] . It is important to 
note the total effect of this deformation is not a simulta- 
neous global rotation of the solid and the nematic direc- 
tor. Therefore as shown in Fig.|4p, long wavelength fluc- 
tuations of these soft deformations are only penalized by 
the Frank free energy for distortion of nematic director. 
As we will show, at finite temperature [25l. [27l. [2^ . and /or 
in the presence of elastomer's random heterogeneity[26j, 
these fluctuations diverge with the system size, and hence 
necessitate a nonperturbative treatment of relevant elas- 
tic nonlinearities. All modes of the elastic deforma- 
tion that are coupled to the soft modes are strongly 
affected by these fluctuations, leading to the aforemen- 
tioned anomalous elasticity. By contrast, because the 
volume remains fixed in soft deformations, the bulk mode 




Fixed Point 



Disorder 



FIG. 3: A schematic renormalization group flows for an 
anomalous solid, such as a nematic elastomer. The elasticity 
of an ideal homogeneous system is controlled by a finite tem- 
perature non-Gaussian fixed point on the vertical temperature 
axis. This fixed point is unstable with respect to infinitesimal 
heterogeneity. The elasticity of a heterogeneous system is con- 
trolled by a zero-temperature, finite disorder non-Gaussian 
fixed point, that endows the nematic gel with its anomalous 
universal elastic properties. 



is not coupled to the soft mode. Consequently the bulk 
modulus of nematic elastomers is only weakly affected by 
fluctuations. 



A 





FIG. 4: Local fiuctuations of soft deformations is the mecha- 
nism of anomalous elasticity in nematic elastomers. (A): the 
reference state of an ideal nematic elastomer. The ellipse de- 
notes the uniaxial nematic order. The light horizontal lines 
are visualization of local translation parallel to the director, 
i.e., the u^ phonon field. (B): a soft volume- preserving defor- 
mation, which costs zero elastic energy in an ideal elastomer. 
It is a simple shear deformation combined with a uniform ro- 
tation of nematic order. (C): long wavelength fiuctuations of 
such soft deformations are the dominant low-energy fluctu- 
ations. Their energy cost is the Frank free energy for splay 
distortion of nematic director field. 



It is worth stressing that the renormalization group 
behavior characteristic of nematic elastomers (and other 
soft solids), where, as displayed in Fig. [31 the low- 
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temperature ordered phase is controlled by a non- 
Gaussian fixed point, is quite exotic from a field-theoretic 
point of view. A more standard scenario is an ordered 
phase controlled by a Gaussian (noninteracting) fixed 
point, with fluctuations and nonlinearities only qualita- 
tively important when the system is fine-tuned to a crit- 
ical point, controlling a continuous transition out of the 
ordered state. In contrast, for nematic elastomers and 
other soft solids [13, M, Hi [H, H [H, HH , fluctuations 
and nonlinearities are important throughout the ordered 
phase, i.e., no fine-tuning is necessary. Consequently, 
the whole ordered phase displays the anomalous elastic- 
ity with nontrivial and universal power-law correlations, 
similar to the behavior near a critical point of more con- 
ventional systems. It is thus appropriate to refer to a 
nematic elastomer and other similar soft solids, like the 
smectic and columnar liquid crystals and flat phase of 
tensionless polymerized membranes, as a "critical phase" 




E. Outline 

In this paper we give a detailed and pedagogical report 
of our recent works on the nonlinear elasticity and fluctu- 
ations in homogeneous as well as heterogeneous nematic 
elastomers. Partial results of this study were published 
earlier in brief communications [25, 26]. The remainder 
of this paper is organized as follows. In Sec. |TT1 we in- 
troduce basic concepts for the development of elasticity 
theory for nematic elastomers and discuss the underly- 
ing symmetries as well as the corresponding Goldstone 
soft modes. Based on these we construct a general elas- 
tic free energy for ideal nematic elastomers, in terms of a 
nonlinear strain tensor which completely encodes the soft 
modes. We also discuss the relation between this elastic 
free energy and the neo-classical theory. 

In Sec. IIII[ we study the elasticity of nematic elas- 
tomers in the context of micropolar elasticity theory. We 
demonstrate that the conservation of angular momentum 
ensures that such medium is characterized by a generi- 
cally asymmetric stress tensor, which encodes angular 
momentum transfer between the translational and the 
internal nematic degree of freedom. By enforcing me- 
chanical equilibrium, we derive the stress tensor and the 
couple-stress tensor of a nematic elastomer as appropri- 
ate derivatives of the elastic free energy. We also explic- 
itly calculate the Cauchy stress tensor and the couple- 
stress tensor for two elastic energy models. 

In Sec. lIVi we formulate the macroscopic elasticity the- 
ory of nematic elastomers at a finite temperature using 
the constant-strain and the constant-stress thermody- 
namic ensembles. These two ensembles are related by a 
Legendre transformation, and contains the same macro- 
scopic physics. We then proceed to construct an effective 
strain-only description of a nematic elastomer by expand- 
ing the elastic free energy in terms of the nonlinear, fully 
rotationally invariant Lagrangian strain tensor of the ne- 



matic state. We derive five Ward identities imposed by 
the underlying rotational symmetry. These identities re- 
late the coefficients of all anharmonic terms to that of 
harmonic ones, as well as require a strict vanishing of 
one of the five elastic harmonic elastic moduli that char- 
acterize a uniaxial solid. 

In Sec. |Vl we use this effective theory to analyze the 
effects of thermal fluctuations on the long-scale elastic 
properties of ideal, homogeneous nematic elastomers. We 
first study the fluctuations at the harmonic level and 
demonstrate that on sufficiently long length-scales, elas- 
tic nonlinearities always become qualitatively important 
for three and lower dimensional elastomers. We then em- 
ploy the machinery of renormalization group transforma- 
tion (RG) to systematically study the long-scale effects 
of thermal fluctuations in the nonlinear theory. Incorpo- 
rating the most relevant nonlinearities, we compute the 
resulting long-scale, universal anomalous elasticity of ho- 
mogeneous nematic elastomers. 

In Sections I VII and IVlIl we study the effects of network 
heterogeneity in nematic elastomers. We first identify the 
most relevant type of quenched disorder and show that 
at a harmonic level it leads to distortions relative to the 
ideal nematic reference state, that diverge for dimensions 
five and below. To characterize long-scale effects of such 
distortions we treat the nonlinear elasticity using renor- 
malization group technique and show that the long-scale 
elasticity of a heterogeneous nematic elastomer is in- 
deed controlled by an aforementioned zero-temperature, 
finite-disorder fixed point. We calculate the renormal- 
ized correlation functions and show that all elastic con- 
stants, except the bulk modulus, are singular functions 
of probing length-scales. Consequently, the stress-strain 
relation is strictly nonlinear for arbitrary small deforma- 
tion. We further show that the spatial correlation of 
the non-affine deformation field acquires anomalous scal- 
ing at long length-scales. Finally we conclude this paper 
with a discussion of experiments and future research di- 
rections. 



II. SYMMETRIES, SOFT MODES, AND 
NEMATO-ELASTIC THEORY 

A. Basic ingredients 

An ideal nematic elastomer is a homogeneous elastic 
medium with internal orientational degree of freedom 
(characterized by the nematic order parameter), which 
may spontaneously break the underlying, internal rota- 
tional symmetry. Similarly to a conventional (fluid) liq- 
uid crystal it thus exhibits a high-temperature isotropic 
phase and a low-temperature nematic phase, which are 
separated by a thermodynamic phase transition [39j . In- 
side the isotropic phase and for small distortions, its elas- 
ticity is described by the conventional elasticity theory of 
isotropic solids P, 0] , characterized by two Lame coeffi- 
cients. 
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A conformation of such d-dimensional ordinary elastic 
medium is described by a d-dimensional vector-valued 
function r{X), which specifies the position of each mass 
point labeled by X, that is also a d-dimensional vector. 
Through out the paper, we will always choose the label- 
ing X such that it coincides with the average position 
of the corresponding mass point in the high-temperature 
isotropic phase of the system. We shall refer to X as 
the isotropic referential, or Lagrangian, coordinate [40| . 
The coordinate space where the vector X "lives" shall be 
referred to as the isotropic reference space and the cor- 
responding state r = X a.s the isotropic reference state 
(IRS). The vector r{X), denoting the spatial position of 
a mass point X in an arbitrary deformed state will be 
referred to as the current or Eulerian coordinate and the 
space in which it "lives" as the embedding space. We shall 
always use symbol a, 6, etc, with a vector on the top for 
all vectors, and use boldface characters such as Q and A 
for matrices and rank-two tensors, use symbols a, 6, etc, 
for unit vectors with magnitude one. 

Using a fixed orthonormal coordinate system 
{ei, 62, . . . , ed = z}, we may decompose the vectors 



as X = '}2a=i^a.ea, and 



Z]j=i^Qea. The elastic 
deformation can be locally characterized by a c? x c? 
matrix function 



^^^a{X) 



drj 
dX'a' 



(2.1) 



which is usually referred to as the deformation gradient 
or as the Cauchy deformation tensor. Note that Kia is a 
"mixed tensor" with one index a in the reference space 
and the other index i in the embedding space. They 
transform differently under rotations in each space. 

Starting from the isotropic phase, the isotropic- 
nematic phase transition is characterized by a develop- 
ment of the nematic order, accompanied by a sponta- 
neous uniaxial strain deformation. The nematic order is 
characterized by a second rank symmetric traceless ten- 
sor field Q(J?), that naturally vanishes in the isotropic 
reference state (IRS). In a low-temperature (uniaxial) ne- 
matic reference state (NRS), it is given by 



Qo = 5 ( fiafiQ - I 



(2.2) 



where the unit vector no is the nematic director that 
gives the uniaxial nematic direction in the NRS and S 
is the magnitude of Q that characterizing the strength 
of nematic order; I is the identity matrix tensor. Fur- 
thermore, the nematic reference state is related to the 
isotropic reference state via a spontaneous uniaxial de- 
formation 



Afio = C± I + (ClI - C-i) n^n^, 



(2.3) 



where Q and C_l are ratios of a deformation in the direc- 
tions parallel and perpendicular to the nematic director 



no. Because typical nematic elastomers (and more gen- 
erally rubber) are characterized by a bulk modulus that 
is much larger than the shear modulus, the spontaneous 
deformation is nearly volume preserving, satisfying a con- 
straint det Ao ~ 1 , that imposes a relation 



CiiCi 



1. 



(2.4) 



by 



The corresponding conformation vector (X) , defined 



(2.5) 



specifying average positions of mass points in the NRS is 
given by 

i^{X) = A;i„ • X - (a I + (Cii - a)^io "o) • X. (2.6) 

To simplify the notations we shall always use x to denote 
the position r^{X) of mass point X in NRS, i.e., 



(2.7) 



and will refer to it as the nematic referential (Lagrangian) 
coordinates and to the space in which it lives as the ne- 
matic reference space. 

An arbitrary state of a nematic elastomer is charac- 
terized by a pair of fields (r(X),Q(X)), which specify 
the geometric deformation as well as the nematic order, 
respectively. Note that we can equivalently treat r and 
Q as functions of the nematic Lagrangian coordinate 
in which case we use the notation (r(a;),Q(x)) [4l|. In 
the nematic phase, it is often convenient to use the defor- 
mation gradient defined relative to the nematic reference 
state: 



A - — 

dxi ' 



(2.8) 



Using Eqs. (|2.5l2.6l2.7p . we readily find the relation be- 
tween A and A: 



A,:„ = 



dr 



dxj 

— Xij (Afjfj)ja, (2-9) 



dXa dXj dXa 



where Ang is the spontaneous deformation, Eq. 
The relations between the IRS and NRS is summarized 
schematically by the upper half of Fig. [5l 



Symmetries and nemato-elastic Goldstone 
modes 



Since the rotational symmetry is spontaneously bro- 
ken in the nematic phase, the Goldstone theorem dic- 
tates that in the nematic phase, the ground state must 
be infinitely degenerate and there must be a soft mode 
associated with this degeneracy. This soft mode was first 
predicted by Golubovic and Lubensky [20] on symmetry 
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IRS NRS 




FIG. 5; Schematic illustration of the soft modes. The 
isotropic reference state (upper left) is spontaneously de- 
formed into the nematic reference state (upper right). The 
ellipse in the center defines the director no. Alternatively, 
one may first rotate the isotropic reference state in upper left 
by a global rotation O into a distinct isotropic reference state 
in the lower left and then have it undergo the I-N transition 
and the corresponding deformation into a physically distinct 
nematic reference state (lower right). The two different ne- 
matic reference states are energetically equivalent, and are 
connected by a soft (Goldstone mode) deformation, as given 
by Eq. IgTSj. 



grounds. A thorough analysis was then given by Olm- 
sted [i^ and in Ref. [1^, and so we only give a brief 
discussion here. As shown in Fig. O starting from the 
IRS (upper left) , there is a continuous manifold of differ- 
ent nematic reference states available to the system. One 
representative is in the upper right, characterized by a di- 
rector fiQ and a spontaneous deformation A^^, Eq. Ij2.3p . 
Energetically equivalent nematic reference states can be 
obtained by first rotating the isotropic reference state by 
an arbitrary orthogonal matrix O (lower left), 

X' = 0-X, (2.10) 

and only then have it undergo an I-N transition to a 
nematic state characterized by a distinct director n (lower 
right), with a nematic order parameter 

Q = 5 (^hh-h^ . (2.11) 

This transition is accompanied by a spontaneous defor- 
mation 

A,, = (ai + (Cii-a)^"). (2T2) 

It is clear from Fig. [S]that starting from one NRS with a 
director fiQ (upper right), the other nematic ground state 
(lower right) with a director n is accessible via a rotation 
of a nematic director, combined with a deformation (right 
arrow in Fig. 

ho^h, AcAf^AftOAri (2.13) 



We note (as long as n is not related to Uq by a rota- 
tion O) this is not a simultaneous rigid body rotation of 
the deformed elastomer and the nematic director in the 
upper right. Instead, the above combination of a defor- 
mation and a director rotation takes the system from one 
ground state to an elastically distinct but energetically 
equivalent one, and therefore is the Goldstone soft mode 
associated with the spontaneously broken rotational sym- 
metry in a solid. Clearly then, in an ideal case the defor- 
mation in Eq. (|2.13p costs zero elastic free energy. This 
soft mode, fully characterized by an arbitrary rotation O 
and the nematic director n, is one of the most striking 
consequence of a spontaneous symmetry breaking, and 
provides the key to many unusual properties of nematic 
elastomers. The main purpose of the present paper is to 
study the local fluctuations of this soft mode and their 
effects on the long length-scale elastic properties of ne- 
matic elastomers. 

An arbitrary uniform deformation of the system away 
from the NRS is described by a triad {no,n,\), where 
no and n are the nematic director of the initial and fi- 
nal state respectively, while A, defined in Eq. (|2.8p . is 
the deformation gradient defined relative to the NRS. 
The elastic free energy density f{ho,n,X) with strains 
measured in the nematic referential coordinates x must 
reflect the symmetries discussed above. 

Namely, / must be invariant both under an arbitrary 
global rotation of the nematic reference state and under 
that of the current state, i.e., it satisfies 

/(no, n, A) = /(Go • no, O • n, OAOj), (2.14) 

where Oo and O are two arbitrary orthogonal matri- 
ces. More importantly, the elastic energy / has to van- 
ish for an arbitrary soft deformation Agm, as given by 
Eq. (i2T3)l . i.e., 

/(no,^,A,,OA7i)=0, (2.15) 

with a special case of the director in the final state the 
same as the initial state given by 

/(no,"o,A,aoOArJ)=0. (2.16) 

More generally, the free energy should be independent 
of whether deformation {ho^h,X) is preceded by a soft 
deformation (no, no, A^j^O Ar^"'^) . Thus the elastic free 
energy density of a ncmatic elastomer must satisfy 

/(no,n, A) = /(no,n,AAft„OArJ). (2.17) 

It clear that Eq. (|2.17p contains Eq. (|2.15p (and hence 
Eq. (j2.16p ) as a special case. 

Let us examine the symmetry properties of the neo- 
classical theory [l^ [13] for ideal, single domain ncmatic 
elastomers. In this theory, the elastic free energy density 
of a uniform volume preserving deformation A is given 
by 

/ = i/xTrlo A^^r^A, detA=l (2.18) 
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where the tensors Iq and 1, defined in Eqs. (|2.22p . are the 
so-called step length tensors and are functions of the ne- 
matic director fiQ in the initial and n in the final states, 
respectively. It is easy to show that indeed the neo- 
classical theory, Eq. (|2.18[) satisfies properties (|2. 14112. 17)) . 
However, it also exhibits another symmetry, namely 



/(no, n, A) = f{no, n, AnO A,^ ^ A). 



(2.19) 



Physically this corresponds to an invariance of the free 
energy under a soft deformation applied after an arbi- 
trary deformation. This property does not seem to follow 
from any underlying symmetry of a nematic elastomer 
and is therefore not expected to hold in real systems. 
Thus, this unphysical constraint and associated limita- 
tions of the neo-classical theory should be experimentally 
detectable. 



C. Invariant strains 

It is a generic feature of nonlinear elasticity theories, 
that there is no unique definition of a strain tensor. For 
ordinary elastic media, where there is no independent 
orientational degree of freedom, infinite number of strain 
tensors can be defined in terms of the deformation gra- 
dient A (H: 



- ((A^A)^ 

m \ 



(2.20a) 
(2.20b) 



where m is an arbitrary integer. A common feature 
of these strain tensors is that they vanish for an arbi- 
trary global rotation. Furthermore, u^'s (v^) are ten- 
sors (scalars) in the reference space, and scalars (tensors) 
in the embedding space. 

For isotropic solids, e.g., ordinary rubber or glass, ei- 
ther u's or v's provide an equally good description of 
elasticity, with elastic free energy expressible as a func- 
tion of any of the tensors above. For anisotropic solids, 
e.g., crystals, tensors u^, provide a natural elastic de- 
scription, while tensors are inadequate to capture the 
anisotropy of the reference state. 

In the case of nematic elastomers, a priori, none of 
these strain tensor provide a complete description since 
they do not contain any information about the nematic 
director (but see Sec. IIV C|) . However, as we now show 
they can be generalized to incorporate the nematic direc- 
tors of the nematic reference and deformed states. Let us 
look at the following generalized symmetric strain tensor 
in the embedding space: 



V 



1 



AloA^ 1 



(2.21) 



where 



lo - Aft„ AT = Ci I + (Cy - Ci) noho, (2.22a) 

T 



Aft A, 



(2.22b) 



are the so-called "step length tensors" in the initial and 
final state respectively. It is easy to see that V is a 
generalization of V2 [44| as defined in Eqs. (|2.20p : in the 
isotropic limit, i.e., 1 = Iq = I, 



AA^ 



I 



V2. 



(2.23) 



The generalized strain V has many nice symmetry 
properties. It is invariant with respect to arbitrary rota- 
tion in the reference space: 



V(Oo • no, n, AOj) = V(no, n. A). 



(2.24) 



It transforms as a second rank tensor with respect to 
rotations in the embedding space: 



V(no, O • n, OA) = OV(no, h, X)0^ . 



(2.25) 



More importantly, it automatically vanishes for arbitrary 
soft deformations Eq. (|2.13p . i.e., satisfies Eq. (|2.15p and 
Eq. (|2.17p . much like the free energy density does: 



V(no,n,AfiOA7i) = 0; 



(2.26a) 



V(no,n,A) = V(no,n,AAft„OAr^;). (2.26b) 

Because of these symmetry properties, we shall call V 
an invariant strain. A scalar function of the tensor V 
automatically satisfies both Eq. p.l4p and Eq. p.lTp . 




Ground State Manifold 

FIG. 6: A schematic of the invariant strain, with the blue 
circle denoting the ground state manifold, no and n are the 
nematic director in the initial and final states respectively. 
Point A^o denotes the nematic reference state (initial state), 
TV the nematic ground state with director n, and A'^' the final 
strained state. The generalized strain tensor V measures the 
deformation of the system state relative to the whole ground 
state manifold, instead of any particular ground state, i.e., 
independent of the choice of A'o. This is the reason that it 
vanishes identically for an arbitrary soft deformation. 



To clarify the origin of the symmetry properties of V, 
we express it in terms of the deformation gradient A 
defined relative to the isotropic reference state. Using 
Eq. dlJl) and Eq. (|2.22ap in Eq. ((2:2T|l we find: 



2V = Ar 



AA^ - AftAi 



(2.27) 
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Remembering that A^i is the energetically preferred spon- 
taneous deformation measured relative to the isotropic 
reference state, if the nematic director of the current state 
is n. Therefore, the tensor inside the square bracket in 
Eq. (|2.27p describes the nonlinear strain V2, measured 
relative to the isotropic reference state, minus its opti- 
mal value for given nematic director n. It is same as the 
tensor 6v defined in reference 29], Sec. IV B (Theory 
with strain and director). In other word, V measures 
the tensorial "distance" between the deformed state and 
the ground state manifold, as illustrated schematically in 
Fig. [HI This is the fundamental reason why V vanishes 
for an arbitrary soft deformation. 

We note in passing that all other nonlinear strain ten- 
sors u„j and v„i, as defined in Eq. (j2.20|l . can be gener- 
alized to the case of nematic elastomers. Both properties 
Eq. (|2.14[) and Eq. (|2.17p are satisfied by the generaliza- 
tion of v,„ for arbitrary integer m. On the other hand, 
analogous generalization of does not satisfy the prop- 
erty Eq. (|2.17p and therefore does not provide an ade- 
quate description of the elasticity of nematic elastomers. 



D. Modes of deformation 

For convenience of later discussion, we shall let the di- 
mension d be an arbitrary integer larger than two. Let 
{n, 171,1, . . .} be an orthonormal basis of the d dimensional 
embedding space, where h is the nematic director in the 
deformed state ^45i] . This coordinate frame therefore ro- 
tates with the nematic director n, in strong contrast with 
the frame {ei, . . . ,ed — z} that we introduced earlier. It 
is convenient to introduce two complementary projection 
tensors in this space: 



n m T 



nn, 
I P. 



that naturally satisfy the following properties: 



P P 

p±p-L 

n n 

I 



P^ 

n 1 

P^Pn = 0, 

P„, + P;^. 



(2.28a) 
(2.28b) 



(2.29a) 
(2.29b) 
(2.29c) 
(2.29d) 



P„ projects out the one dimensional subspace along the 
nematic director n, while P^ projects out the d — 1 di- 
mensional subspace perpendicular to n. 

The strain tensor V can then be decomposed into dif- 
ferent parts using these projectors: 



P„VP„ 


— V V 


(2.30a) 


PnVP;^ 


= nVn±, 


(2.30b) 




= Vn± n. 


(2.30c) 


P^VP^ 

n n 




(2.30d) 



n 

m 

T 



Vnn 




Vn. 





FIG. 7: Decomposition of the invariant strain tensor V into 
three parts: Vnn, V„i, and Vx. 



Using this basis, as illustrated in Fig. [7l we decom- 
pose the strain V into an isotropic part and some other 
traceless parts according to 

1 



v--(Tvv)i-f y„„j 



where J is a d x d traceless matrix 

1 

- 





Vn±\ 






\Vn± 


j~ 







(2.31) 



(2.32) 



Ij^ a (d — 1) X (d — 1) identity matrix, and a scalar Vnn 
and a (d — 1) x (d— 1) traceless matrix Vj^ are given by 



-(TrV) 



id -I) 



(2.33a) 
(TrV_L)I_L. (2.33b) 



To illustrate the geometric meaning of the different 
components of the tensor V in Eq. I2.31i let us consider 
a special case of small, spatially dependent deformation 
from the nematic reference state. To lowest order ap- 
proximation, we need not distinguish the Lagrangian and 
Eulcrian coordinates. We define a phonon field m, a lin- 
earized symmetric strain tensor e, and a linearized anti- 
symmetric strain a in a standard way: 



il{x) — f{x) — X, 
dr-i 

■^m = T\ — Sia + daUi 

dXa 

aab = ^ (daUb - dbUa) , 



(2.34a) 
(2.34b) 

(2.34c) 
(2.34d) 



where we have adopted the nematic Lagrangian coordi- 
nates and have used the shorthand da for d/dxa- We 
also choose fio = i-d = z, and consider small fluctuation 
of the nematic director: 



This decomposition is illustrated in Fig. [T] 



n — fiQ + 6h — z + 6h, \Sn\ <C 1. (2.35) 
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A) Bulk Mode B 
z 



B) Soft Mode lOzi 

y 




C) Longitudnial Shear)a,L D) Transverse Shear |a, 

FIG. 8: Schematics for a bulk, soft, longitudinal shear, and 
transverse shear modes, as well as associated rigidities. Differ- 
ent components of the strain tensors that correspond to these 
deformation modes are listed in Table |ll] At lowest order, the 
bulk, longitudinal shear and transverse shear modes are the 
same as ordinary uniaxial solids. The soft mode mixes simple 
shear with rotation of the nematic director and thereby qual- 
itatively distinguishes a nematic elastomer from an ordinary 
uniaxial solid. 



The full expressions for components of V in terms of 
these variables can be obtained by using Eqs. ()2.34|) and 
Eq. (|2.35p in Eq. (|2.2ip . Up to linear order in u and dh, 
this leads to: 



TrV w Tie = V 



(2.36) 



Consequently TrV describes the isotropic compres- 
sion/dilation of the solid, i.e., the hulk mode, illustrated 
in Fig. [8l\. We also find that 

= K„, -^TrV 



■Tie 



dzUy. - -(V • m), 
a 



(2.37) 



where we have used the fact that z k, h sX the lowest 
order. Therefore Vnn describes a dilation (compression) 
along the h axis combined with an isotropic compres- 
sion (dilation) transverse to this axis (perpendicular sub- 
space), with the total volume conserved, as illustrated in 
Fig. [8p. We shall call this mode the longitudinal shear 
mode. The strain component V_l is, up to linear order 
of the phonon field, given by 



k=l 



d-l 



(Vi • U^) 5t^{2.i%) 



where indices i j and k are restricted to the {d— 1) dimen- 
sional transverse subspace. Therefore Vj^ describes shear 
deformation in the d— 1 dimensional subspace perpendic- 
ular to n, as illustrated in Fig. [8jl). This mode breaks the 
uniaxial symmetry and we will refer to it as the transverse 
shear mode. Up to a linear order, the above three modes, 
i.e., the bulk, the longitudinal shear, and the transverse 
shear only depend on the symmetric linearized strain e 
and are the same for nematic elastomers as for ordinary 
uniaxial solids. 

However, the strain operator Vn± mixes the symmetric 
strain e with the antisymmetric strain a and the director 
fluctuation Sn. Upto linear order in u and Sn, it is given 
by 



(a„i - Sfii) 



(?.39) 



where r = J- is a dimensionless ratio characterizing the 
spontaneous anisotropy of the nematic phase. In the limit 
of a vanishing nematic order, C-L ^ C||: ^ ~* Ij ^nd Vni 
consequently reduces to the linearized symmetric strain, 
Sni- For a non- vanishing nematic order, however, the 
strain Vni involves a symmetric strain Sni, a rotation of 
the solid a„i , as well as a rotation of the nematic director 
Sn. Strain components Vni_ qualitatively distinguish a 
nematic elastomer from an ordinary uniaxial solid. It is 
clear from Eq. (|2.39p that Vn± describes a simple shear 
in a plane containing n and the axis, compensated 
by an infinitesimal rotation of the nematic director as 
well as a global rotation of the solid, and vanishes for 
an infinitesimal soft deformation, where the cancellation 
between these three contributions is complete. This is 
illustrated in Fig. [8(3. 



E. Nemato-elastic model 

The invariant strain tensor V measures deformations 
relative to the uniaxial ground state manifold of an ideal 
nematic elastomer, and therefore vanishes for arbitrary 
soft deformation. Thus, for a uniform deformation, the 
elastic energy of an ideal nematic elastomer can be ex- 
pressed as a scalar function of V. Indeed, in the neo- 
classical theory, the elastic free energy per unit volume 
[i^ . Eq. (P7T5)1 . can be written as 

/nco- ^MTrloA^riA^^iTrV+^/i, (2.40) 

where we have used the cyclic property of the trace opera- 
tion and det 1 = det Iq that follows from Eq. (|2.22p . Using 
Eq. (|2.2ip . the incompressibility constraint det A = 1 can 
also be written in terms of V as 



det(I + 2V) = 1. 



(2.41) 



Because of the highly nonlinear incompressibility con- 
straint Eq. (I2.4ip . the neo-classical theory Eq. (|2.40p is 
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not particularly convenient for our goal of studying long 
wavelength fluctuations of a nematic elastomer. Fur- 
thermore, as we have noted at the end of Sec lllB| the 
neo-classical theory exhibits extra unphysical symmetry, 
Eq. (|2.19[) , and is thus not the most general theory satis- 
fying physical elastomer symmetries Eqs. p.l4p . (j2.17|) . 
We therefore relax the incompressible constraint and ex- 
pand the elastic free energy in powers of the invariant 
strain V. This removes the aforementioned unphysical 
symmetry of the neo-classical theory. Linear terms can- 
not appear since the nematic reference state is a mini- 
mum of the elastic free energy. In accordance with the 
uniaxial symmetry in the nematic phase, five different 
quadratic terms can appear. Ignoring all higher order 
terms, the elastic energy density can be expressed as [13] 

fe = iB.yl + A,iK„(TrV^) + iA(TrV^)2 

+ MTrVi+Ai„±K?^ (2.42) 

=== i B (Tr V)2 + C (Tr V) -f i VL 

+ /iTrVi+Ai„±(14±)', (2.43) 

where we have used two different decompositions of V, 
related to each other by Eqs. (|2.33p . Note that we have 
used the coordinate system {n, m,l, . . .} as illustrated in 
Fig. [7l The total elastic free energy is then given by /e 
integrated over the nematic referential volume element 

Fe^ J fed''x = J J-^Jed\ (2.44) 

where 

[)t ■ 

J — det A — det — — ^ 

is the Jacobian factor relating the Eulerian volume ele- 
ment to the Lagrangian volume element. For most rub- 
bery materials the Jacobian factor is very close to unity 

(up to 10-5). 

The two sets of parameter in Eq. ^(I^ and Eq. ^(I^ 
are related via: 

B, = B + 2(l-i)C+(l-^)VL, (2.45a) 
a a 

Xni_ = S + (l-^)C-^(d-l)Mi, (2.45b) 
A ^ (2.45c) 

Inverting the relations (|2.45p . we find 

B = ^ {B, + {d-l){2C+{d-l)X + 2fi)()%A6a.) 
d^ 

ML = -2C + A+-^, (2.46b) 
C = i (B;, + (d - 2)C - (d - 1)A - 2/i) . (2.46c) 



Following our discussions of different deformation 
modes in the preceding subsection, we shall call the elas- 
tic constants B, fij^ and fi in Eq. ([2:431) the bulk modulus, 
the longitudinal shear modulus, and the transverse shear 
modulus, respectively. Gencrically, rubber is character- 
ized by a bulk modulus B that is four to five orders of 
magnitude larger than all other elastic moduli. Further- 
more, as we shall see in Sec. |V] and Sec. IVIl ideal ne- 
matic elastomers are actually strictly incompressible due 
to strong fluctuations that singularly renormalize elastic 
moduli, driving all (/xl, and C) but the bulk modulus 
B to zero at long length-scales. This essential differ- 
ence between the bulk modulus B and other moduli (/iL, 
fi, C) indicates that Eq. ([23T|) and Eq. f^A3^ provide 
a more natural parameterization to the invariant strain 
V and the elastic energy, respectively. By contrast, if 
one insists on using the parameterization in Eq. (|2.42[) 
with independent parameters B^, Xz±, A, /i, one would 
find, from Eqs. f^A5\i that aU renormalized elastic con- 
stants except /X flow to the same limiting value, i.e., to 
the renormalized non-universal bulk modulus B. The 
essential physics (i.e., strict incomprcssibility, universal 
Poisson ratios, etc.) is thus missed in this latter pa- 
rameterization, hidden in the differences between these 
renormalized constants and B, which go to zero in the 
long wavelength limit. 

We note that by virtue of the invariant strain V, the 
elastic energy (I2.43|) vanishes identically for an arbitrary 
soft deformation, i.e., it satisfies Eqs. (|2.15l - (|2.17p . On 
the other hand, since the elastic energy is expanded up 
to the second order of the strain V, it is valid for small 
value of V, i.e., the strained state should be not far away 
from the ground state manifold. Because of the uniax- 
ial nematic order, the elastic energy ()2.43p of a nematic 
elastomer bears a resemblance to a conventional uniax- 
ial solid, that to a quadratic order in the nonlinear La- 
grange strain tensor u (using the decomposition intro- 
duced above) is given by [1] 

fe = ^Ci + C4 ii„„Tr u_L + i C2 (Tr m_l)^ 
+ CsTiul + C^iu^nf. (2.47) 

Although both elastic energies, Eq. (|2.42p and Eq. (|2.47p 
are characterized by five independent moduli, there is 
an essential qualitative difference between a conventional 
uniaxial solid and an ideal nematic elastomer: while the 
Lagrange strain u is a tensor in the reference space and 
only depends on the geometric deformation r{x), the in- 
variant strain V is a tensor in the embedding space and 
depends both on the geometric deformation and on the 
nematic director rotation, and more importantly, van- 
ishes for arbitrary soft deformations. 

For non- uniform deformations, the elastic energy 
Eq. (|2.42p and Eq. (I2.43P must be augmented by the free 
energy cost for distortion of the nematic order, that in 
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3d is given by the Frank free energy 



fd 



fd (TX 



i i^i (V • + 1 (n • V X nf 
1 



(2.48) 



where Ki, K2 and are the splay, twist, and bend- 
ing constants, respectively. Again, we have defined the 
free energy density fd as the Frank free energy per unit 
volume measured in the nematic referential coordinate 
system. Here it is important to note that the gradient 
operator V is defined in the Eulerian coordinates r, i.e., 
Vi = d/dri, because the nematic director interaction is 
dominated by the liquid fraction [4§| . 

For small deformation, we can take n « z -\- 6n and 
make linear approximation of the invariant strain V in 
Eq. (|2.43p . in terms of linearized strains e, a and the 
director fluctuation 5h. The elastic energy Eq. (|2.43p 
then reduces to 



/e - iB(Tr£)2 + C(Tr£)e, 



1 ~2 

2 A^L 



i-L Tr e]_ + /i„_L 



where 



Tr£ 



-aa ^zz 1 / ,^ZZ) 
1=1 



r - 1 



d-1 



(2.49) 

2 



(a^i - Shi) 



Szz - ^Tr£, 
a 



1 



- 1 



(2.50) 

(2.51) 
(2.52) 

(2.53) 



The last term in Eq. (|2.49p is the most interesting one. 
It explicitly demonstrates our earlier general symmetry 
assertion that a simple shear e„i can be completely com- 
pensated by a combination of the elastomer and the ne- 
matic director rotations, a„i and Sh respectively, such 
that the net elastic energy cost is zero. An infinite set of 
such soft elastic distortions is specified by a constraint 



1 



(a„i - 6n,) ^ 0. 



(2.54) 



This complete compensation is a consequence of the 
spontaneous symmetry broken nematic elastomer state 
and is the lowest order manifestation of the soft defor- 
mation Eq. (1133]). 

If one is primarily interested in the elastic degrees of 
freedom, then director fluctuations Sri can be integrated 
out of the Eq. (|2.49p , thereby obtaining an effective elas- 
tic description purely in terms of the strain e. It is easy 





a(l+exx) 

FIG. 9; A simple experiment on a uniaxial nematic elastomer. 



to show that at the quadratic order this leads to a rota- 
tion tensor a and Sn that automatically adjust to e to 
satisfy the soft mode constraint Eq. (|2.54p . Thus it is 
quite clear that the resulting effective elastic free energy 
density has the form Eq. (|2.49p . but with the last fin± 
shear term identically vanishing, a result that was first 
derived (in a quite a different way) by Golubovic and 
Lubensky [2§| . We thus observe that the elasticity of 
uniaxial nematic elastomers resembles that of ordinary 
uniaxial solids, but with key distinction that the trans- 
verse shear elastic modulus C5 strictly vanishes, which is 
guaranteed by the underlying rotational symmetry spon- 
taneously broken in the nematic state. 

To illustrate the basic structure of the harmonic theory 
Eq. (|2.49p . let us consider a simple mechanical experi- 
ment illustrated in Fig. [HI With the initial nematic order 
(uniaxial axis) along z, we study the elastic response to 
a small strain deformation Exx — £• To avoid the elastic 
stripe instability |50|], we take e to be negative, as shown 
in Fig. [HI To find the equilibrium state we minimize the 
elastic free energy Eq. (|2.49p with respect to all strain 
components except Exx- In strong contrast to a stretch 
along ix axis, a compression along is stable against an 
inhomogeneous soft deformation, and therefore Sn and 
strain components Eni as well as a„i remain zero. We 
can therefore ignore the term in Eq. (|2.49p . 

For concreteness we focus on three dimensions, and 
simplify the analysis by taking the bulk modulus B to 
infinity, a reasonable approximation for elastomers that 
are essentially incompressible. Minimizing the elastic en- 



ergy Eq. arm 

we obtain: 



over Ey 



and Ezz for fixed e^ 



-2 Id 



(2.55a) 

(2.55b) 
(2.55c) 



which satisfy the traceless incompressibility condition 

3 



Tr£ = ^. 



0. 



(2.56) 



As expected, in this simplest (fiuctuation-free) har- 
monic treatment the compressional elastic response co- 
incides with that of an ordinary uniaxial solid. Because 
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in the incompressible limit the cross-coupling C drops 
out, the ratios between different strain components (i.e., 
Poisson ratios) only depend on a single ratio As 
we shall see later in this paper, long wavelength fluctua- 
tions drive this ratio to a universal constant, set by the 
ratio between the renormalized shear moduli /i^ and /i^ . 
The simple experiment illustrated in Fig. [9] allows a di- 
rect test of the universal Poisson ratios predicted by our 
theory. 

Although our focus here (and throughout the 
manuscript) is on the most experimentally relevant 
uniaxial nematic state, our phenomenological model, 
Eq. p.42p or Eq. p.43p . can be easily generahzed to 
a biaxial nematic elastomer. Indeed we can define the 
invariant strain tensor V by Eq. (j2.2ip for arbitrary m- 
axial nematic order simply by using corresponding m- 
axial step-length tensors Iq and 1. However, such non- 
uniaxial nematic elastomer generically has many more 
elastic constants in its free energy. For example, a three- 
dimensional biaxial nematic elastomer is characterized 
by biaxial step length tensors Iq and 1. All rotational 
symmetries are broken and the corresponding quadratic 
elastic free energy density has 9 independent parameters: 

f biaxial = C*! F,j„ -1- 6*2 + C3 

+ C4 Vnn Vmm + C5 Vnn Vu + Cg Vmm Vll 

+ CrV^„, + C^V^i+C^Vli, (2.57) 

where {n, fn, V\ are three mutually orthogonal principle 
axes of the biaxial nematic order in the deformed state. 
Similarly a two-dimensional ideal nematic elastomer is 
characterized by 4 elastic constants and the following 
elastic free energy density: 

fbiaxial — C*! ^fn + + ^3 VnnVmm + C4 Vmm 

(2.58) 

where {to, n} are the principle axes of the nematic order. 

We conclude this section by contrasting the differences 
between the neo-classical elastic model, Eq. (|2.18[) and 
the nemato-elastic model, Eq. (|2.43p . While both the- 
ories can be formulated in terms of the invariant strain 
V and both capture the full nematic elastomer symme- 
try, the neo-classical theory has the advantage of simplic- 
ity (characterized by only two independent parameters, 
the shear modulus jji and the dimensionless ratio C||/C-i)j 
and can be easily used to describe large deformations. 
On the down side, as mentioned earlier, the neo-classical 
model is over-constrained, containing a larger symmetry, 
Eq. ((2?T9)) than that exhibited by a generic nematic elas- 
tomer. The neo-classical theory is also not well-suited for 
studying long wavelength fluctuations. 

On the other hand, although the Laudau theory 
Eq. (|2.43p does not have this redundant symmetry, char- 
acterized by five elastic moduli, it can only be used for 
small deformations near the ground state manifold. In 
the limit of infinite bulk modulus, the moduli B and 
C drop out and the Landau theory considerably simpli- 
fies to three independent parameters. Furthermore, it 



allows for a systematic study of long wavelength fluctu- 
ations and heterogeneities as well universal properties of 
nematic elastomers. 



III. COSSERAT, MICROPOLAR, AND 
ASYMMETRIC ELASTICITY 

A. Micropolar elastic media 

Nematic, as well as other liquid crystalline elastomers, 
are examples of a wider class of micropolar elastic media 
[51 )1 that are elastic solids with both translational and 
rotational degrees of freedom. A proper description of 
micropolar elasticity therefore falls outside the domain of 
the classical elasticity theory . Related to the concept 
of micropolar elastic media, there is another concept of 
micropolar fluids, i.e., fluids with broken rotational sym- 
metry. Examples of micropolar fluids include nematic, 
hexatic, and ferroelectric fluid phases, all characterized 
by an asymmetric stress tensors. In between micropolar 
solids and micropolar fluids, there are also systems with 
orientational order and partial translational order. For 
example, in a smectic liquid crystal, in addition to the 
mass-density wave position, whose fluctuations are speci- 
fied by a phonon field u{x), the orientation of the nematic 
director field n{x) must be specified to fully character- 
ized the state. In the limit of infinitely strong anchoring 
(which can be shown to be always the case in the long 
wavelength limit), the director h{x) becomes locked to 
the layer normal, thereby allowing an effective descrip- 
tion in terms of purely translational elastic degrees of 
freedom. Although a three-dimensional material, a smec- 
tic liquid crystal is not a three-dimensional solid, since it 
has rigidity only in one direction alone the layer normal. 

Historically [24,], the elasticity of micropolar elastic 
solids was first studied by W. Voigt and the Cosserat 
brothers back in early 1900's, and is often referred to 
as Cosserat media. In the standard micropolar elastic- 
ity theory, a deformation of the solid is described by the 
displacement field x{X) as well as an independent ro- 
tation vector </)(X), with latter specifying the orienta- 
tion of the "molecules" (which are usually modeled as 
isotropic rigid bodies) relative to some reference config- 
uration. A distinguishing macroscopic property of such 
media is that they can transmit torque independent of 
force. Not being widely noticed for a long time, this sub- 
ject was revived about half a century later, most notably 
by Nowacki, Eringen and others [5ll |. in the mechanical 
engineering and applied mathematics communities, and 
was renamed as micropolar elastic media. One major 
motivation underlying these efforts seems to study the 
elasticity of solids with microstructures. The theories 
developed by these authors are quite formal and with 
presentation that parallels the classic work of elasticity 
by Love [l|. They focus on conservation laws, as well as 
general equations of motion, without paying much atten- 
tion to the microscopic physical mechanism that leads to 



14 



the unconventional elasticity. The success of these the- 
oretical efforts on micropolar elasticity is limited by the 
lacking of well controlled experiments on relevant ma- 
terials. We note that only quite recently did physicists 
start to successfully synthesize these elastic solids with 
internal degree of freedom. Furthermore, most of the re- 
search on the micropolar elasticity was concentrated on 
the isotropic state of the micropolar media, where the 
coupling between the relative rotation and shear defor- 
mation, e.g., the last term in Eq. (|2.49p or Eq. (|2.43p . 
does not appear. As we will show in subsequent sections, 
this nemato-elastic coupling is in fact the most interest- 
ing and most important feature of micropolar elasticity. 

For nematic elastomers, the internal orientational de- 
gree of freedom is characterized by a unit vector field n(x) 
with two degree of freedom, not a rotation vector 4>{x) 
with three degree of freedom. In this section, we shall ex- 
amine the elasticity of nematic elastomers from the view 
point of micropolar elasticity theory, with the conserva- 
tion of angular momentum as the central principle. Our 
discussion shall thereby clarify various aspects of qual- 
itative difference between ordinary solids and nematic 
elastomers. We will also derive equations that must be 
satisfied by the stress tensor and the coupled-stress ten- 
sor in elastic equilibrium, obtained from the variation of 
the elastic energy with respect to a virtual deformation 
and rotation. From these equations (working in 3d) we 
will find explicit expressions for the true stress tensor and 
the couple-stress tensor for a uniformly deformed nematic 
elastomer. 

Following the standard convention of the nonlinear 
elasticity literature, in this section we shall use a slightly 
different notation than the rest of the paper. We will use 
X to denote the referential (Lagrangian) coordinate and 
X the target space (Eulerian) coordinate [52|. A state 
of an ordinary elastic media is then given by a function 
x{X,t) where t is the time. Let dX, dA, dV be the line, 
surface, and volume elements of the Lagrangian coordi- 
nates. Then the corresponding quantities in the Eulerian 
coordinates dx, da, dv are given by standard relations 

dx = XdX, (3.1a) 
da = J\^^ ■ dA, (3.1b) 
dv = JdV, (3.1c) 

where J — \dx/dX\ is the corresponding Jacobian factor. 

A local physical quantity g, e.g., a certain extensive 
quantity per unit mass, can be treated as a function of 
either the Lagrangian or the Eulerian coordinates, i.e., 

g{x,t)=gix{X,t),t). (3.2) 

We use g to denote a material time derivative, i.e., the 
time derivative of g with the Lagrangian coordinate X 
fixed, and use dg/dt for the time derivative with the Eu- 
lerian coordinate x fixed. Clearly these two derivatives 



are related to each other in a standard way 

9 = ^^+v-Vg, (3.3) 

where 

v = j^x{X,t)^S{X,t); (3.4) 

is the velocity field. 

As the elastic medium deforms, the region Vl{t) that it 
occupies in the Eulerian coordinate space changes with 
time. The corresponding region fio in the Lagrangian 
coordinate space remains constant in time. Taking the 
mass density po j measured in the Lagrangian coordinates 
to be spatially uniform, we have: 

4 / gpodV^ [ gpodV, (3.5) 
at J no J no 

since the material time derivative commuting with inte- 
gral over the Lagrangian coordinates. The mass density, 
p, measured using the Eulerian coordinates is related to 
Po by 

pdv — pj dV — padV — >p=J^^po- (3-6) 

We can use this to express Eq. (|3.5p in terms of the Eu- 
lerian coordinates: 

/" gpdv^f gpdv^ [ (^ + v-Vg] pdv. 

"-t Jn{t) Jn{t) Jn{t) \ot J 

(3.7) 

Eqs. (13.71) and (|3.5p are referred to as the transport the- 
orem [53| in elasticity literature. 

B. Conservation of angular momentum in 
conventional elastic media 

Let X (Eulerian coordinate) be an arbitrary point in- 
side an elastic media and da ~ dda a surface element 
passing through x, where a is the unit normal vector. 
The stress vector field t{x, da) is defined as the elastic 
force that is exerted on one side of the surface element 
da from the other side. The famous Cauchy's theorem 
[l], [5^ then dictates that the stress vector t depends lin- 
early on the vector da and there exists a stress tensor 
field T(x), referred to as the Cauchy stress, or the true 
stress, given by 

t{x, da) — T(a;) • da, ti ~ Tijaj. (3-8) 

As we will see, for ordinary liquids or solids, conservation 
of angular momentum requires the Cauchy stress tensor 
to be symmetric, i.e., = Tji. We can also express the 
stress vector in terms of the Lagrangian surface element 
dA: 

t{x, da) = T(x) • da = S(X) • dA, (3.9) 
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through the so-called nominal stress tensor S [53j , related 
by Eq. ((XTb)) to T via: 

S = J TA-'^, T = S\^. (3.10) 

The nominal stress is generally not symmetric. 

For a deformed body, Newton's second law dictates 
that the rate of change of total linear momentum is given 
by the total force acting on the body: 

/ vpdv = / fpdv + [ T(f) • da, (3.11) 

™ Jn(t) Jnit) Jdn{t) 

which consists of the bulk and surface contributions, with 
/ the body force per unit of mass acting at a point x. 
Applying the transport theorem on the left hand side, 
as well as the Gauss theorem on the right hand side, we 
find the following equation of motion for an (ordinary or 
micropolar) elastic solid: 

/ pildv^ [ (pf{x)+\/-T)dv. (3.12) 
Jn{t) Jnit) ^ ' 

where V • T is a vector with components dTij/dxj. 
Since this must be true for an arbitrary volume element, 
Eq. (|3.12p must be satisfied by equality of integrands and 
thus leads to a local field equation of motion 

pi/ = p/(f)+V-T, (3.13) 

referred to as Cauchy's first law of motion. For an elastic 
body in equilibrium and under a vanishing body force, 
we recover the well known result that static stress field 
is divergenceless in the absence of external forces, i.e.. 

The angular momentum of an ordinary solid is due 
to the center of mass motion of its constituents (atoms, 
molecules, colloids) and is given by 

Le = pdvx X V, (3.14) 
Jn 

where volume v is not to be confused with velocity v. 
The time derivative of the angular momentum is given 
by the total torque acting on the elastic body: 

[ A - - f ^ - ^ 
— Lf, = — / pav X X V = / pdv XXV 

at dt Jfj 

= J dv X X {pf + V -T), (3.15) 

where we used the transport theorem, Eq. (j3.7[) and 
Cauchy's first law of motion, Eq. p.l3p . Defining a 
pseudo- vector Ta associated with the antisymmetric part 
of T by: 

{TA)i = iijkTjk, (3.16) 

it is easy to see that 



Using this identity together with Gauss's theorem in 
Eq. ((XT5)) wc find 

= I dv{j)xx J + Ta^ + I xxT-dd. (3.18) 

For a vanishing body force /, the bulk contribution to 
angular momentum must be locally conserved, with 
only an external surface stress contributing to a change 
in Le- From above expression this is clearly true if and 
only if Ta = 0. We therefore find Cauchy's second law of 
elasticity, namely, that in an ordinary solid the Cauchy 
stress tensor T must be symmetric. 

C. Conservation of angular momentum in nematic 
elastomers 

A nematic elastomer is an anisotropic micropolar elas- 
tic medium, and its complete description is character- 
ized by the mapping x{X) and the nematic director field 
n(X). The time dependent variation of the nematic di- 
rector gives an additional intrinsic contribution, Li to 
the total angular momentum, 

Ltot^Le + L,. (3.19) 

The extrinsic part Le, defined in Eq. (|3.14p . gives a con- 
tribution associated with the motion of the center of mass 
of the molecules. The interaction between the transla- 
tional and orientational degrees of freedom allows an ex- 
change between the intrinsic and extrinsic angular mo- 
menta, and thus, neither is separately conserved [s^ . 
The total angular momentum, Ltot is of course always 
conserved. Going back to Eq. p.lSp . we see that for ne- 
matic elastomers, the antisymmetric part of the Cauchy 
stress tensor, Ta is generally nonzero, describing the pro- 
duction of the extrinsic angular momentum due to the 
coupling between translational and orientational degrees 
of freedom. 

To study the intrinsic angular momentum, we need 
to introduce two physical quantities special to microp- 
olar media. Adopting the terminology from the elas- 
ticity literature [5l|, the body torque, G, is the exter- 
nal torque (e.g., due to an external magnetic field) per 
unit of volume acting on the nematic degrees of freedom. 
The couple-stress vector c{x, da) is defined as the torque 
transmitted through an infinitesimal surface element da, 
solely due to the distortion of nematic director across the 
surface. Clearly this torque acts on the orientational de- 
grees of freedom, i.e., the nematic director [55|. One can 
extend the Cauchy's theorem and show that the couple- 
stress vector also depends linearly on da and therefore 
can be written as the contraction between da and the 
couple-stress tensor C{x): 



[x X (V • T)], = di {e.jkXjTu) + {Ta\. (3.17) 



c{x,dd) ^ C{x) ■ da. (3.20) 
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The couple-stress tensor describes the torque transmitted 
through a unit surface area due to a distortion of the 
nematic order, instead of a geometric deformation [s^ . 

From a general principle, the rate of change of the total 
angular momentum is given by the total torque, which is 
given by combination of the body force and body torque 
acting on the bulk, and by the stress and couple-stress 
vectors acting on the boundary: 



—rJ^tot 
at 



pdv (x X f + g) + [ (xxT + C)- da 
^ ' Jon 

(3.21) 



Subtracting Eq. ((3?T8l) from Eq. ((3?2T|) . wc find the rate 
of change of the intrinsic angular momentum 



= j pdv {g ~fj^ + j C da. (3.22 



In the absence of body force and torque as well as the 
boundary stress and couple-stress, Eqs. (|3.18p and p.22|) 
reduce to 

dLe_ ^ = I pfiy ^ j pfiy BieijkTjk, (3.23) 



dt 



demonstrating that the antisymmetric part of the Cauchy 
stress tensor describes angular momentum transfer be- 
tween translational and rotational degrees of freedom of 
a micropolar elastic medium. 



D. Equilibrium conditions, stress and couple-stress 

The total free energy of an ideal nematic elastomer 
consists of the elastic energy as a function of A and h, 
as well as the Frank energy Eq. p.48p . External body 
forces, such as those due to e.g., a gravitational or mag- 
netic fields, can also be easily included. For example the 
free energy associated with the interaction a magnetic 
field h and the nematic director field h(x) is given by 



Fm = / fmJ-'dv = -:r7a / (^(f) • h)' J' ' dv (3 .24) 

that induces body torque. We note that the energy den- 
sity fm is defined with respect to unit volume in La- 
grangian coordinate and the integration is over the Eu- 
lerian coordinates x. 

To derive the equilibrium conditions, let us apply ex- 
ternal force ii .^j (x) daj and torque ii %j (X) daj on the 
surface element da at the boundary. We consider an in- 
finitesimal variation of the equilibrium state, i.e., trans- 
lation of the center-of-mass Eulerian coordinate as well 
as rotation of the nematic director field as follows. 



X 

n{x) 
5fi{x) 



x' = X + 5x, 
n {x ) = n{x) ■ 
5ljj{x) X h{x). 



(3.25a) 
5n{x), (3.25b) 
(3.25c) 



Operationally, we first rotate the nematic director by 
5lij{x) without deformation, and then translate each point 
X (together with the rotated director fi + 5n) to x' with- 
out rotation of the director field. The external traction 
(stI on the boundary introduces two following terms in 
the total free energy: 



boundary 



an 



(x) 5x^ + (f) buJiix)) d4^,.26) 



The force couples to the translational degree of freedom 
8x while the torque couples to the rotational degree of 
freedom bh. The total free energy of nematic elastomer 
under external boundary traction is then given by 

Ftot^ \ {fe + fd + frn - l^{x)n^) J^^dv + bFboundary, 

Jn 

(3.27) 

where for convenience, we included a Lagrange multiplier 
term to ensure that the constraint that n is a unit vector 
field is satisfied. 

In equilibrium, the first-order variation of the total free 
energy must vanish. Since the referential volume element 
J^^dv — dV in Eq. (j3.27p is invariant under the variation 
of the current coordinates x x' , the variation of the 
total elastic energy is given by 



6Fu 







{Sfe + Sfd + 5fm - v{x) flibhi) J ^dv 



in 
5F, 



boundary ■ 



(3.28) 



Let us look at the variation of the different part of 
the bulk free energy separately. The elastic free energy 
density fe[\,n\ is a function of the deformation gradient 
^ia — dxi/dXa and the nematic director n. (It is also a 
function of the initial director no, which does not change 
in the variation.) Under the variation Eqs. p.25p the 
deformation gradient changes as follows: 



dx', 



b^ia — dabXi — 



dxj dbxt 

dXa dXj 



+ ^o.bx^, (3.29) 
= Xjadjbxi, (3.30) 



where to avoid confusion, we use Roman indices i, j, fc, 
etc., to label Eulerian coordinates and use Greek sub- 
scripts a, (3, etc., to label Lagrangian coordinates. The 
first-order variation of the elastic free energy is therefore 
given by 



bFe = 



dfe 



bhi 



dfe 



J 



dh 

-idfe 

dfii 

1 a/, 



J 



bfi 



Xjadjbxi ] J ^dv 



dXj \ dXia 



dn 



dX,,, 



■ \j(ySxi da j 



dv 



(3.31) 
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The Frank free energy density /^[n, Vn] is a function 
of both the director h and its gradient Vn. Under the 
variation (|3.25p . the gradient Vn changes as fohows: 



dhi 



dx'- 



dx^ d 
dx'j dxk 



{■hi + Sfii) , 



(3.32) 



which, using Eqs. p.25p . to the lowest order gives 

djSfii — djSxkdkfii (3.33) 



on the surface element with normal Cj, it has to be bal- 
anced by the corresponding component of elastic force, 
i.e., the Cauchy stress tensor. We immediately see that 
the remaining terms in the surface integrand must be the 
Cauchy stress tensor Ty , giving us an important consti- 
tutive relation for the stress tensor field: 



ddj hi 



(3.38) 



The first-order variation of the Frank free energy is there- 
fore 



For the convenience of later discussion, we separate the 
stress tensor into two parts. 



djSf 



rri ry-re 



rpa 



(3.39) 



n dvj 



where 



'nj^ 



dhi 
dxi 



(3.34) 



(3.35) 



where due to the elastic deformation, and Tf^ due to 

■/ . . . •' 

a distortion of the nematic director field are given by 



7-1 9fe . 



'jk 



(3.40a) 
(3.40b) 



Finally, the magnetic free energy density fmin], as 
given by Eq. (j3.24p . is a function of the nematic director 
only. The variation of F„i is therefore given by 



SF„ 



■Shi J ^dv 



For a conventional solid, there is no Frank free energy 
fd, and therefore the stress tensor is solely given by 
in Eq. (j3.40ap . Furthermore, rotational symmetry in the 
embedding space implies that the elastic energy must be 
a function of the metric tensor defined as 



In dh, 

= —^aHi I {n • h)Shi J~^dv. 
Jn 



(3.36) 



9af3 



dx dx 



dXa dX 
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(3.41) 



Su bstitu ting Eq. ^M, Eq. ([331, Eq. (p^S]) and 
Eq. p.26p into Eq. (|3.28p . and performing integration 
by parts appropriately, we arrive at 



which then gives for the Cauchy stress tensor of an ordi- 
nary solid 



an 



daj 



1 




In 





dfe ^ ^ 
dXia 



dfd 
dfd 



Vik 



dVjk 



Vik 



£^ij (x) 

Sxi dv 



1 df 

Tij = Tij — 2 J — \irv\ 



Sxi 



dgap 



(3.42) 



da-i 



an 



d fd 

J~^— — Shi — ^ij{x)Suji{x) 



n 

dh 
dh 

0. 



— v{x)h 



drin 
dh^dfd 
dhi dhi 



Jdj I J 



-idfd 



which (as expected on general grounds discussed above) 
is obviously symmetric. 

Let us now calculate the stress tensor Tf^ for an ideal 
nematic elastomer. For the neo-classical elastic energy, 
Eq. (|2.40p . to impose the incompressibility constraint, we 
introduce a Lagrange multiplier p, physically correspond- 
ing to the pressure. We then find that the elastic stress 
tensor is given by 



Shi dv 



T" = n\io \^ -pi. 



(3.43) 



(3 



At mechanical equilibrium, each volume and surface in- 
tegral in the left hand side must vanish separately. 

Let us first look at the surface integral Eq. (|3.37p . 
which is linear in Sx. Since by definition ^ij is the z— th 
components of the external force per unit area acting 



37) which is clearly asymmetric, unless 1 commutes with 
A Iq a""" . The pressure term is always isotropic and sym- 
metric regardless of the deformation, as one would ex- 
pect. 

A different form of the stress tensor Tf. is obtained 



for our phcnomcnological elastic model, Eq. (|2.43p . A 
standard but somewhat tedious calculation gives 
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(2V + 1) 



{BTtY + C K«) I + (C Tr V + K„)(™ 

I 



(3.44) 



which is also generically asymmetric. 

Let us next turn to Tf^ in Eq. (|3.40b|) . derived from the 
Frank free energy. It is non- vanishing only for an inhomo- 
geneous configuration of the director field. Except for the 
isotropic pressure term, which can be associated instead 
with Tlj , Tfj is identical to the Ericksen stress tensor 
studied extensively in nematic liquid crystals literature. 
Using Frank free energy, Eq. (|2.48p inside Eq. (j3.40b[) we 
find 

Tfj = - Ki{V ■ n)dinj - K2{n ■ V x n)ejkinidink 



K3nj{dink){n ■ \7nk), 



(3.45) 



that is also generically asymmetric, unless all Frank elas- 
tic constants are the same. 

Let us now look at the bulk term (volume integral) 
in Eq. p.37p associated with the variation Sx . It is 
clear that the integrand describes the total force per unit 
volume acting on the element dv. Indeed the integrand 
is identical to —djTij /Sxi, with Tij the Cauchy stress 
tensor given by Eq. p.38p . Mechanical equilibrium then 
implies 



T = 0, 



(3.46) 



corresponding to a force balance equation of ordinary 
solids, in the case of a vanishing body force. 

Similarly, the integrand in the surface integral of 
Eq. (|3.37p involving orientational degrees of freedom 
must also vanish. Since the tensor is the i— th compo- 
nent of the external torque per unit area acting on the the 
surface element with normal e^, it must be balanced by 
the corresponding component of the couple-stress term, 
defined in Eq. (0201) • Using Eq. (|3.25cp to relate Sn to 
SlS, we see that the couple-stress tensor Cij must be given 
by 



Cij(x) = J EjM- Tlfc. 



(3.47) 



Substituting the Frank free energy Eq. (|2.48p into 
Eq. ((OT)) we find 

Cy = -Ki{V ■ n)eijknk- K2{n-V y. h){5ij - niUj) 
+ K^e^kinjnk{n-V)ni. (3.48) 

Similar to the stress tensor T^^ , the couple-stress tensor 
Cij is nonzero only if the nematic director field is non- 
uniform. 

Finally, we look at the integrand of the bulk term in 
Eq. (|3.37p linear in 5n. Defining a local molecular field 
acting on the nematic director by 



hi{x) 



dfe df, dU , jaJj-l|Z^i)(3.49) 



dfij dhj dhj 



bulk torque balance requires a vanishing of this inte- 
grand. This then gives a torque balance equation 



hi{x) + v{x)fii = 0, 



(3.50) 



which imposes an equilibrium orientational constraint 
that the local molecular field must be parallel to the ne- 
matic director. We do not analyze the molecular field h 
in detail in this work. 



IV. NONLINEAR ELASTIC MODEL OF 
HOMOGENEOUS ^^M.ATIC ELASTOMERS 

A. Stress versus strain ensembles 

The main objective of this paper is to study the ef- 
fects of long wavelength fluctuations on the macroscopic 
elasticity, e.g., the stress-strain relation, of an ideal ne- 
matic elastomer. To formulate the problem, we need to 
consider the relation between fluctuations and finite tem- 
perature elasticity in some detail. At finite temperature, 
mass points in an elastic media fluctuate around their 
equilibrium positions, regardless of whether or not the 
system is coupled to an external traction. On the other 
hand, an external traction causes a macroscopic deforma- 
tion, which changes equilibrium positions of mass points. 
Furthermore, the spectrum of the fluctuations may also 
be modified by a macroscopic deformation. 

The central quantity we want to calculate is the elas- 
tic free energy, as usual obtained as a sum over all elas- 
tic configurations with an appropriate Gibbs-Boltzmann 
weight, determined by the elastic energy discussed in the 
previous section. It then determines a macroscopic elas- 
ticity, renormalized by long wavelength fluctuations. The 
free energy can be computed either as a function of the 
imposed macroscopic strain tensor and the nematic di- 
rector, or as a function of the external traction. The 
choice of the macroscopic control variable (strain or trac- 
tion) amounts to a choice of a statistical ensembles within 
which to study thermal fluctuations. As usual, we expect 
that in the thermodynamic limit, macroscopic quantities 
are independent of the choice of the ensemble. 

Let us first look at the constant strain ensemble by 
fixing the macroscopic deformation gradient Aq and the 
nematic director fiQ. For simplicity, we shall only study 
the case where Aq and no are uniform. We shall switch 
back to the notations we introduced in Sec. |TT1 using x 
for the nematic referential coordinates and r{x) for the 
deformed (target) coordinates. To compute the thermo- 
dynamics (e.g., the partition function) we need to sum 
over all fluctuations of the position field r{x) and the 
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nematic director field n(x) around their equilibrium val- 
ues. Naturally, we parameterize these two fields in terms 
of the phonon displacement field, u{x) and the nematic 
director fluctuation, Sn(x), defined by 

■r{x) = 'ro{x) + u{x) ^ Xo ■ X + u{x), (4.1) 
n{x) = hQ + Sh{x). (4.2) 

Expressing the Landau elastic energy, Eq. (|2.43|) in terms 
of these fields, we obtain the elastic "Hamiltonian" func- 
tional H[Xo,ho,u,Sh]. Since u(x) and 5h(x) are local 
fluctuations, they have well-defined Fourier transforma- 
tions [sgj . The corresponding thermodynamic potential, 
given by 

F[Ao,no] = -T log J Du{x)DSh{x)e-"^^'>'^'>'^'^^^/'^, 

(4.3) 

is the elastic free energy, as a function of the macro- 
scopic deformation Aq and the average nematic director 
no- In above functional integrals a short distance (large 
momentum) cutoff is implicit, set by the lattice constant 
in crystalline materials and polymer network mesh size 
in elastomer materials and we are using units in which 
the Boltzmann constant ks is 1. 

To calculate the functional integral in Eq. (|4.3|) , we first 
expand the elastic Hamiltonian around its minimum for 
the given Xq and no, which, baring the possibility of an 
elastic instability, is given by u = and Sh — 0, i.e., a 
uniformly deformed reference state. The zero-th order 
(tree level) approximation to the functional integral is 
then simply given by the saddle point approximation, 
i.e., by setting it = and (5n = in the integrand of 
Eq. jts]), i.e., 

F[Ao,no] ~ H[Xo,no,u = 0,6h = 0] = Hmft[Xo,ho]. 

(4.4) 

This is just the Landau theory, that does not take into 
account any long wavelength fluctuations. 

To systematically address the effects of fluctuations on 
the macroscopic elasticity beyond the Landau (mean- 
field) approximation, we expand the elastic energy in 
terms of the phonon field u{x) and the director fluctu- 
ation Sh, and calculate the functional integral Eq. (|4.3p 
order by order in corresponding nonlinearities. These 
fluctuations are generically important for all soft solids. 
As we discussed in Sec. U for the particular case of ideal 
nematic elastomers, the fluctuations of nemato-elastic 
Goldstone modes lead to qualitative breakdown of per- 
turbation theory in the thermodynamic limit, as well as 
necessity of a renormalization group analysis. 

We can also study the constant traction ensemble, by 
introducing a boundary traction term, expressed in terms 
of the Lagrangian coordinates and the external traction 
tensor S^^ J(,3^A"'^),q using Eq. (j3.1bp : 

6H^ - [ Sr,{x) S^c^ix) dA^, (4.5) 



For simplicity we shall focus on the case of a vanishing 
external torque traction, "joij — 0. 

In the case where the external traction tensor S is in- 
dependent of X, we can use Gauss's theorem to rewrite 
6H lea] as 

5H^~I dVS,a.da,n = - f dVS^c^X,^, (4.6) 

a form that conflrms that the deformation gradient Ao 
and the traction tensor S are conjugate variables. There- 
fore, in equilibrium, S is equal to the nominal stress ten- 
sor, as introduced in Eq. (|3.10p . In fact this explains the 
notation we used in Eq. (|4.5p . 

The thermodynamic potential for this ensemble can 
then be calculated by summing over all elastic conflgu- 
rations with the modified Gibbs-Boltzmann weight: 

G[S] = -Tlog J Dr{x)Dh{x)e-^"+^"~^'^ , (4.7) 

As usual the thermodynamic potential G[S], Eq. (|4.7p 
and the elastic free energy Eq. (|4.3p are related by a 
Legendre transformation 61]: 

G[S] = (F[A] -FTrSA)|Ao, (4.8) 

where V is the system volume and Ao is the solution to 
the following equation 

— (F(A) - V Tr SA) = —F(X) -V^ = Q. (4.9) 

This is precisely the relation between the macroscopic 
nominal stress tensor and the macroscopic deformation 
gradient [s^ , both of which can be directly measured in 
experiments. 

B. Referential coordinates at finite temperature 

As defined in Sec. [ill the nematic reference state min- 
imizes the elastic Hamiltonian and thus precludes the 
appearance of terms linear in the invariant strain ten- 
sor V. However, in the presence of thermal fluctuations 
and elastic nonlinearities, the average position {f{x)) of 
a mass point is not the same as its nematic referential 
coordinate x. Since the nematic referential coordinate x 
and the average position {r{x)) differ only at finite tem- 
perature, their difference (to lowest order) must be pro- 
portional to T, and inversely proportional to an elastic 
modulus. 

For an isotropic solid, symmetry dictates that x and 
{r{x)) can only differ by a scalar factor and therefore the 
associated modulus must be the bulk modulus, B: 

{r{x)) — X « const. — - — x, (4-10) 
B 

where ^ is the microscopic cutoff length-scale. For ordi- 
nary isotropic rubber we have T^~'^/B <C 1 and there- 
fore the difference between the referential coordinate and 
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the average position of a mass point is negligible. In con- 
trast, for an anisotropic elastomers such as a nematic 
elastomer, the relation between (r(a?)) and x involves 
other, significantly smaller elastic constants, which are 
comparable with T^^''. Therefore, generically thermally 
induced deformation of the average position f(x) is non- 
negligible, particularly near the I-N transition. For a 
given elastic Hamiltonian, this relation can be systemat- 
ically computed and in principle is experimentally acces- 
sible. 

Because of this non-negligible difference between x and 
{r{x)) there are two (in principle equivalent) complemen- 
tary descriptions of the elastic theory of a nematic elas- 
tomer. One, most conceptually clear formulation is in 
terms of phonon modes expanded about a coordinate sys- 
tem defined by {f{x)). The advantage of such approach 
is that the computed renormalized elastic free energy is 
by construction minimized by a vanishing phonon dis- 
placement, i.e., precludes any linear terms in the non- 
linear strain tensor V. A technical disadvantage is that 
the elastic Hamiltonian (that does not include long wave- 
length fiuctuations) will then necessarily admit terms lin- 
ear in V, whose coefficients are adjusted order by order 
in perturbation theory so as to ensure that the renormal- 
ized elastic free energy is free of such linear terms [6^ . 
The latter is thus quite cumbersome to work with. An 
alternative, technically more convenient approach is to 
choose the nematic reference coordinate system x, which 
ensures that instead the elastic Hamiltonian is free of 
terms linear in the strain. The unavoidable price of this 
is that such linear terms will then be generated in the 
renormalized elastic free energy, as a signal of a shift in 
the elastic minimum configuration. These terms will be 
discussed in much more detail when we carry out the 
renormalization group analysis in later sections [63| . 

Thus, for convenience we will always choose the ne- 
matic referential coordinates x to label mass points and 
about which to expand phonon fluctuations. Neverthe- 
less, our final results in this work will of course be inde- 
pendent of the choice of the reference state, [s^ 



nonlinearities that survive at long length-scales. Further- 
more, focusing on the elastic degrees of freedom, we will 
integrate out the nematic director fluctuation in the par- 
tition function and obtain an effective theory in terms of 
only the translational degree of freedom r(x). This op- 
eration is well justified if the nematic director is tightly 
coupled to the elastic media, such that the fluctuaftions 
of the nematic director from its energetically favored di- 
rection (determined by the elastic degrees of freedom) 
are small. This condition is always satisfied [13] and is 
quite similar to the locking (at long length-scales) of the 
nematic director fluctuations to the local layer normal 
fiuctuations in a smectic liquid crystal, thereby allowing 
a well-known purely elastic description of smectic liquid 
crystals. 

A "strain-only" elastic model can be obtained by inte- 
grating out the nematic director fluctuations from the 
Hamiltonian (|2.43p . formulated in the nematic state. 
However, to clarify the role of the spontaneous symmetry 
breaking and to simplify the analysis it is more conve- 
nient to do this in a model formulated in the isotropic 
phase, i.e., expressed in terms of conformational degrees 
of freedom relative to the isotropic reference state and the 
nematic order parameter Q. Integrating out the nematic 
order parameter, Q leads to an effective elastic energy 
functional in which the effective shear modulus changes 
sign at the I-N transition [2^, thereby inducing a 
uniaxial distortion of the elastomer matrix. Expansion 
about the resulting nematic reference state then gives 
the sought-after purely elastic description of the nematic 
elastomer. 

Our goal is to deduce the most general form of the re- 
sulting elastic Hamiltonian, carefully ensuring its under- 
lying rotational and translational invariance. To ensure 
these conditions we take it to be the most general scalar 
function of the metric tensor g, deflned in Eq. (|3.4ip . 
Without loss of generality, we can express it as a func- 
tion of the following scalar quantities 

5„ = Trg" = Tr(ATA)", n= 1,2,3,... (4.11) 



C. Strain-only nonlinear elastic model 

We are interested in the long wavelength elastic fluc- 
tuations in nematic elastomers, i.e., fluctuations of the 
phonon and the nematic director fields around their 
equilibrium values. In thermal equilibrium, the Gibbs- 
Boltzmann distribution of these fluctuations is controlled 
by the coarse-grained elastic Hamiltonian functional, 
Eq. (|2.43p . augmented by the Frank energy Eq. (|2.48p 
for the distortion of nematic order. Being fully invariant 
under rotations as well as arbitrary soft deformations, 
this energy functional is rather complicated, containing 
a large number of nonlinearities, most of which are irrel- 
evant at long length-scales. To analyze effects of fluctu- 
ations, we therefore look for a simpler, minimal model 
formulation, which from the start contains only those 



which are rotationally invariant both in the reference 
space and in the embedding space. In three dimensions, 
the metric tensor has only three eigenvalues, ensuring 
that first three invariants Sn are independent |63]. 

Hence the elastic Hamiltonian is a function of Si, S2, 
and 5*3: 

H[g]=H[Si, 82,83]. (4.12) 

In the high-temperature isotropic phase, the elastic 
Hamiltonian is minimized by g = I, for which 81 = 82 = 
83 — 3, obviously corresponding to the isotropic reference 
state. In the low-temperature broken-symmetry phase, 
the elastic energy is minimized by the nematic reference 
state (NRS), characterized by a uniaxial metric tensor 

go = AjAo, (4.13) 
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where Aq, given by Eq. (|2.3p . is the spontaneous defor- 
mation accompanying the I-N transition. Inside the ne- 
matic phase, we can expand the elastic energy Eq. (|4.12p 
around the NRS: 

n.M = / [^1,52,^3] 

= Hei [S'l + 5 Si , 3° + SS2 , si + SS3] 

= f[SlSlSl]+J2^nSSn 

71 

+ i^*™„<55„,55„ + 0[<553], (4.14) 



Substituting Eqs. (|4.18p into Eq. (|4.14p and rearranging 
the series in the order of increasing powers of the La- 
grange strain tensor e, we obtain an effective elastic free 
energy for nematic elastomers with spontaneously broken 
rotational symmetry: 



1 
2 



(4.21) 



^ [Bz el^ + 2 \z± e^z en + A en ejj + 2^ efj) 



2 

•IZ JZ 



where 



<55„ = 5„-^o=Tr(g"-gJJ), 



dSjidSm 



(4.15a) 
(4.15b) 

(4.15c) 



and we have truncated the series at the second order 
in 6Sn- Since, by definition SSn's are fully rotationally 
invariant, the rotational symmetry as well as the associ- 
ated soft mode of the nematic phase are guaranteed by 
every term in Eq. (|4.14p and are not compromised by the 
truncation. 

To make contact with elasticity theory we express the 
scalars SSn in terms of the more familiar nonlinear La- 
grange strain tensor defined relative to the nematic ref- 
erence state: 



Sob 



1 f df df 

2 \dxa dxb 

i {daUb + dbUa + daU ■ dbu) , 



(4.16) 



(4.17) 



where in the last expression we used Eq. (|2.34ap to ex- 
press e in terms of the phonon field u. To this end we 
make use of the relation Eq. (|2.9p between deformation 
tensors in the isotropic and nematic reference states, find- 



ing: 



551 = 2Trloe, (4.18a) 

552 = 4Trl2e + 4Tr(loe)2, (4.18b) 
12Tr(l2eloe) + 8Tr(loe)3. 



(553 = 6Trlge 



where we defined 



In 



AoAj. 



(4.19) 



Choosing the nematic director of the NRS to be along 
the z axis, we may express the tensor Iq as a matrix: 



go 



Ci 
Ci 

Cii 



(4.20) 



where repeated subscripts z, j, k are summed over trans- 
verse components, x, y only. In above, we have only kept 
the most important terms, i.e., those that are relevant in 
the renormalization group sense, as we will verify a pos- 
teriori. The values of the eleven coefficients appearing in 
Eq. (|4.2ip are listed in Appendix El 



Not all of the eleven coefficients in Eq. (|4.2ip are inde- 
pendent to each other. It is straightforward, although 
tedious, to check that these coefficients, as given by 
Eqs. (|Aip . satisfy the following five Ward identities: 



bi 
b2 
bs 
c 



2 (a — (1 + a) aj_), 
2aBz-2{l + a) {Xz± ^ 
2aXz± -2(1 + a) A, 
2 a ^z± — 4 (1 + a) /i, 
2^2 + 2(l + a)2(A^ 
4a(l + a) iXzA_+l^z±), 



2m) 



(4.22a) 
(4.22b) 
(4.22c) 
(4.22d) 

(4.22e) 



where 



Ci 



Ci 



7T (4-23) 



are two dimensionless ratios characterizing the 
anisotropy of the nematic phase. Therefore out of 
eleven there are only six independent parameters in 
Eq. (g^. The five Ward identities, Eqs. ^j^^ reflect 
the underlying rotational symmetry in the isotropic 
reference state, spontaneously broken by the nematic 
state. 

We note that although we have derived these identities 
in three dimensions, as we explicitly show in Appendix [Bl 
they actually hold in arbitrary dimensions d > 3, gener- 
alized by simply allowing indices k, ... to range over 
all d — 1 transverse directions. The elastic Hamiltonian, 
Eq. (|4.2ip and the set of Ward identities Eqs. (|4.22p agree 
with Eq. (4.3) and Eqs. (4.4) of Ref. [l^, after appropri- 
ate redefinition of various constants. The correspondence 
is established in Table |T1 
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This work 












A 




bi 


&2 


&3 


c 


Reference [28.1 




Cia2 








2C± &3 


C1&4 











TABLE I: Correspondence between the Ward identities in this paper and those in reference [28l ]. The parameter s used in 
reference [28l | corresponds to ((^^ — Cx)/2 in our notations. 



As we have already discussed in Sec. IIVBI it is con- 
venient to choose C± and C|| in Eq. (|4.20p such that 
the nematic reference state e = minimizes the elas- 
tic Hamiltonian Eq. (|4.2ip . This ensures that the coeffi- 
cients and a± of linear terms are exactly zero, which 
further reduces the number of independent parameters in 
Eq. (|4.2ip to four. More importantly, this choice, when 
combined with the first Ward identity Eq. (|4.22ap . dic- 
tates that the shear modulus /izj_ strictly vanishes in the 
nematic phase, as we have shown in a complementary 
way in Sec. [TTl As emphasized there, a vanishing of fiz± 
without any fine-tuning is a consequence of the underly- 
ing rotational invariance of the isotropic reference state 
and is the lowest order manifestation of the correspond- 
ing soft Goldstone mode. This feature qualitatively dis- 
tinguishes a nematic elastomer from an ordinary uniaxial 
solid, whose corresponding modulus C5 in Eq. (I2.47P is 
generically nonzero. 

The other four Ward identities, Eqs. (|4.22bj|4.22ep . 
dictate that the coefficients of all anharmonic terms in 
Eq. (|4.2ip , 61 , 62 , 63 and c, are completely determined by 
the remaining 4 quadratic coefficients, the elastic mod- 
uli Bz, Xz±, As we will see in the next section and 
as is clear from general symmetry principles from which 
these arise, these relations are preserved in the presence 
of fluctuations. Using these relations, we can define an 
"effective" nonlinear strain tensor w with components: 

Wzz = ezz + 2aeli, (4.24a) 
w^j = e,j - 2{l + a) eziCzj, (4.24b) 

where a is defined in Eq. (|4.23p . and indices i and j are 
limited to the perpendicular subspace l,2,...,(i— 1. The 
elastic Hamiltonian Eq. (I4.2ip then reduces to 

Tiel = ]^Bzwlz + ]^\ wfi + \zl_ Wzz Wii + Wij ,(4.25) 

which, after supplemented with the curvature energy cor- 
responding to the Frank free energy of the nematic di- 
rector, will be the form we will use for further study of 
thermal fluctuations and network heterogeneities. 

1. Two-dimensional nematic elastomer 

The same analysis can be carried out straightforwardly 
for a two-dimensional nematic elastomer. An impor- 
tant property special to the two-dimensional case, is that 



the subspace perpendicular to the nematic order is one- 
dimensional, and thus has no shear mode (characterized 
in higher dimensions by the modulus fi) associated with 
the transverse subspace. 

Consequently, the effective elastic Hamiltonian for a 
2d nematic elastomer is given by: 

= ^ >'XX U'L + \ \y + Xxy Wyy, (4.26) 

where we chose the nematic director no to be along the 
X-axis. The nonlinear strain components Wxx and Wyy 
are 

Wxx = ex^ + 2a ely, (4.27a) 
Wyy = Cyy -2{1 + a) ely, (4.27b) 

where a — CyliCx ~ Cy)- Wc will leave the analysis of 
this very interesting two-dimensional model for a future 
investigation. 

D. Relevant nonlinearities and the minimal 
strain-only elastic model 

As we shall see later in this work, the nonlinear elastic 
model, Eq. (|4.25p is an interesting coupled combination 
of a smectic-like and columnar-like nonlinear elasticities, 
involving Uz and w^, respectively. Experience with the 
nonlinearities in the presence of fluctuations in these two 
systems [13, IH, llBl suggests that it is the former that are 
more relevant. We will verify this rigorously a posteriori 
in Sec. lV Al Thus, anticipating a stronger relevance oiuz 
elastic nonlinearities to long length-scale fluctuations, we 
substitute Eq. (|IT7)) into Eq. ((TM)) and Eq. (gSSl), and 
only keep those anharmonic terms (smectic-like nonlin- 
earities) that are proportional to (u^, m^, ii_LU^), ignoring 
all others. It is straightforward to show that this amounts 
to making the following replacements in Eq. (|4.25p : 

Wzz ^ dzUz + ^{Vi_Uzf, (4.28a) 

1 a 

Wij ^{diU-j + djUi) - - {diUz djUz), (4.28b) 

where a is defined in Eq. (|4.23p . 

In order to streamline the notations, we further rescale 
the phonon fields by 

Uz a~^Uz, (4.29a) 
uj_ a-^u±, (4.29b) 
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and absorb an overall factor a ^ of the elastic Hamilto- 
nian by redefining all the elastic constants: 

{B, , , A, m) ^ a' (S„ , A, Ai) . (4.30) 

In order to stabilize soft-mode phonon fluctuations, 
the Hamiltonian in Eq. (|4.25p must be augmented by 
the the Frank free energy, Eq. (|2.48p . using the relation 
Eq. (|2.54p to eliminate director fluctuations 5n in favor of 
the phonon fields. Furthermore, since we are most inter- 



ested in the fluctuations of the Uz phonon field, which in 
the momentum space is controlled by the pole Qz ^ q^, 
it can be easily shown that both the bend (K^) and twist 
(K2) term are (dangerously) irrelevant in the RG sense 
We will therefore ignore these two terms for the 
purpose of the RG analysis, and to simplify the notation 
will use K to denote the Ki splay modulus. 

The resulting minimal elastic Hamiltonian density has 
the following form: 



1 



K 



'Hel = -BzWzz+>^z^WzzWn + -\Wii+ ^.Wi.+ — {V ^Uz) 



ML .2 



= - (Tr ^y + C (Tr w) Wzz + ^w: 



y ■ 2 

H Wij Wij 



(4.31a) 
(4.31b) 



where the components of the (rescaled) effective strain 
tensor w now become 

Wzz = dzUz + ^iW^Uzf, (4.32a) 

Wij = ^ {diUj + djUi - diUzdjUz) , (4.32b) 

and 

Bulk Trw = Wzz+Wii, (4.33a) 

Longitudinal Shear Wzz = Wzz Trw, (4.33b) 

a 

Transverse Shear Wij — Wij — — — — 6ij (4.33c) 
I 



The elastic constants B (bulk modulus) , /iL (longitudinal 
shear modulus), /i (transverse shear modulus) and C have 
already been defined in Eqs. (|2.45p . 

The elastic nonlinearities contained in different modes 
of deformation are particularly interesting. Qualitatively 
speaking, the relevant elastic nonlinearities arise from the 
coupling between various modes of deformation with the 
soft mode. It is the fluctuations of the soft mode that 
renormalize various elastic constants and lead to anoma- 
lous elasticity. We observe that the bulk mode 

Trw = Wzz + wu ^ dzUz + diUi = V • u (4.34) 



Moduli 


B 


ML 




IJ-z± 


Descriptions 


Bulk Mode 


Longitudinal Shear 


Transverse shear 


Soft Mode 


V 


TrV 


Kn = Kn-d"'(TrV) 


Vx 


Kx 


w 


Trw 


Wzz — Wzz — (Trw) 




N/A 


e 


Tre 


ezz=ezz-d-^ (Jve) 


Sij 


Ea. (|2.39|l 



TABLE IL Various modes of deformation and their operators in different descriptions 



does not contain any relevant nonlinearities. This is due we observe that longitudinal shear mode 
to the fact that the soft deformation Eq. (|2.13p preserves 
the volume and therefore its local fluctuations do not cou- 
ple to the bulk mode. Consequently we expect that renor- 
malization of the bulk modulus is qualitatively unimpor- 
tant. On the other hand, from Eqs. and Eqs. 



Wzz = (1 - d-^) dzUz ~ d-^ d,u, + {y^uzf (4.35) 
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and the transverse shear mode 



Wi 



diUj + djU, 



diU:^djUz 



(4.36) 



do involve relevant elastic nonlinearities. We therefore 
expect that the renormalization the longitudinal and 
transverse shear moduli by fluctuations to be qualita- 
tively important. These expectations will be verified by 
explicit RG calculations in the next section. 

We emphasize the structural similarity between 
Eq. (|43Tb)) . Eq. (EISD, and Eq. dHH), with the only 
exception that the term with the coefficient in 
Eq. and Eq. ((^^ is replaced by the term with 

the coefficient K in Eq. (|4.31b[) . Indeed the elastic free 
energy Eq. (j4.3fb[) can be obtained from Eq. (|2.43p , aug- 
mented by the Frank free energy Eq. (|2.48p . by integrat- 
ing out the fluctuations of nematic director, and drop- 
ping irrelevant nonlinearities. This procedure replaces 
the term by the curvature K term. 

Finally we note that we have represented the bulk, lon- 
gitudinal shear, and transverse shear modes in terms of 
the invariant strain V, the linearized strain tensors e and 
a (Sec. HI]), as well as the effective strain tensor w (the 
current section) . At the linear order in the phonon field, 
they all agree with each other. In contrast, because the 
soft mode involves both strain deformation and rotations 
of nematic director, it cannot be expressed in terms of w 
alone since the latter does not contain 5n. Thus, the soft 
mode only appears indirectly in the current w descrip- 
tion, through an exact vanishing of fizi. and the form of 
the nonlinear strain tensor w. The relations are summa- 
rized by Table. |TT1 

The effective strain-only model, Eq. (|4.31[) will be the 
starting point for all of our further analysis of thermal 
fluctuations in the presence of elastic nonlinearities. 



V. LONG-SCALE THERMAL ELASTICITY OF 
A HOMOGENEOUS NEMATIC ELASTOMER 

In this section we use the effective elastic model derived 
in the previous section, Eq. (|4.3ip to study the long- 
scale properties of a nematic elastomer in the presence of 
thermal fluctuations. 



A. Harmonic phonon fluctuations 

We begin by studying thermal fluctuations within a 
harmonic approximation, ignoring all elastic nonlinear- 
ities in the model, (j4.31bp . As mentioned above, this 
is equivalent to going back to the harmonic theory, 
Eq. (|2.49p . supplemented by Frank free energy, and in- 
tegrating out the director fluctuation 5h. To the lowest 



order at long wavelength this amounts to the nemato- 
elastic coupling in Eq. (|2.49p simply enforcing the re- 
placement 



r + 1 

; 

r~l 



-■711 ' 



(5.1) 



in the Frank free energy Eq. (|2.48p . The resulting 
quadratic effective elastic energy is then given by 

H^,i = Jd^x i [B, {d^u^f + 2 (d^u,) (9,u,) 

+ K,{Vlu,)^ + K3id^ui_)^], (5.2) 

where we have absorbed some constant factors into the 
Frank elastic constants Ki (splay) and (bending). 
The twist term {K2) turns out to be less relevant than 
Ki and K3, and therefore is ignored in Eq. (|5.2p . Except 
for the last two terms (from Frank free energy) , Eq. (|5.2p 
is identical to the quadratic (in phonon field u) parts of 
the nonlinear elastic Hamiltonian Eq. (|4.25p . 

In terms of the Fourier transform of the phonon field 
'u{x) 



u{q) = J <rxu{x)e-''>-\ (5.3) 

the harmonic elastic Hamiltonian iJo takes the form 

J d\ TaM uM uti-q), (5.4) 

where indices a and b are summed over all d dimensions, 
and the kernel matrix Tab{<t) is given by 



i] = + fJ')qiq] + {^J-q±+ K3q^)dij, (5.5a) 



Az_L qz qi, 



r.. - B,ql+Kiqi. 



(5.5b) 
(5.5c) 



The harmonic phonon correlation functions 

{uMM-<i))o = TC7.,(g)w,(-Q')e-«»["l/^. 

(5.6) 

can be calculated using the equipartition theorem or 
equivalently by performing above simple Gaussian inte- 
gral. They are given by 

{uMub{-(!))^ = {2nY5''(q-(^)Gab{q). (5.7) 
where as usual Gabiq) is related to Tab{q) via: 

Gabiq)TbM^TSac. (5.8) 

The components of the propagator matrix Gabiq) are 
fairly complicated. For example, Gzz{q) is given by 



T 



(5.9) 
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where Cz{q) is a wavevector-dependent constant 



where 



CM - 



)l2 



(5.10) 



Since we are most interested in small q (long length-scale) 
properties, the q dependent term K3, q^/q^ in the denom- 
inator of Eq. (|5.10p is generically much smaller than the 
q independent term A -I- 2//, and therefore can be ignored 
[6g | . Therefore we may approximate Gzz (q) as 



GzM 



T 



fiql + Ki q\ 



4 ' 



(5.11) 



where 



X2 



A + 2 



= CMU.=o (5.12) 



The correlator Eq. (|5.1ip is identical to the harmonic 
phonon correlation function of a conventional smectic liq- 
uid crystal [58| with fi and Ki the modulus for layer 
compression and layer bending, respectively. In three di- 
mensions, the real space mean-squared fluctuations of the 
Uz phonon field are given by: 



{uzm^ 



(Pq^dqz 



T 



(27r)3 



oc 



—j= log — 



^J'<lz 
L 



K^ql 



(5.13) 



where L is the system size and a the small cutoff length- 
scale, set by the molecular size. The fact that {uz{r)'^)o 
diverges with the system size suggests a breakdown of 
the harmonic elasticity theory and the qualitative im- 
portance of elastic nonlinearities that have so far been 
ignored. As in the case of smectic liquid crystal, we ex- 
pect that the elasticity of three- (and two-) dimensional 
nematic elastomers is dominated by long wavelength fluc- 
tuations of the Uz phonon field. 

We can use Eqs. (|2.45p to express jl as function of B, 
C, /iL and fi. In the limit of an infinite bulk modulus, i.e., 
_B — > 00 with C, /iL and fi fixed, we find fi approaching 
a finite limit 



2{d-2) 



as B 



(5.14) 



Therefore in this limit, the correlation function of Uz field 
simplify considerably, becoming independent of the bulk 
modulus B. We will see that the same result holds for 
other correlation functions as well. Physically, for a large 
B, the bulk mode is essentially frozen out, and as a result 
does not play any role in long wavelength fluctuations. 

Matrix inversion of Tab also gives the harmonic corre- 
lations of uj_{q): 



G,j{q) ^ GMPl^Aq, 



GT{q)P7M)., (5.15) 



GM 

GT{q) 
C±iq) 



T 



C^{q)ql+K^qr 
T 



Hql + K^qi 
X + 2fi 



^z± 



Bz + Kiqi/qr 



and 



9 ' 



QiQj 
1I 



(5.16a) 
(5.16b) 
(5.16c) 

(5.17) 



are the longitudinal and transverse projection operators 
(with respect to qi) in the (d — 1) dimensional subspace 
perpendicular to z. Eq. (|5.15p is similar in structure 
to the harmonic phonon correlations of a columnar liq- 
uid crystal [ssl . [3a | . A straightforward calculation shows 
that real space fluctuation of u± remains finite in three 
dimensions. Therefore the long wavelength fluctuations 
of the phonon fields, u± are qualitatively unimportant 
as compared with those of Uz- This qualitative differ- 
ence between the phonon fields Uz and u± is a result 
of spontaneous broken rotational symmetry in an origi- 
nally isotropic system. By contrast, in ordinary uniaxial 
solids, all phonon fluctuations remain flnite in three di- 
mensions, thus leading only to quantitative corrections 
to properties of conventional solids. 

The cross correlation functions between Uz and u± are 
given by 



where 



Gz:{q) - V-\uziq)u,i-q))o, 

^ ^zijq) 

Tzziq)TLiq) - {TzAq)r 



rL(g) = (A + 2Ai) ql+K,qi, 



(5.18) 
(5.19) 

(5.20) 



and V is the volume of the system. For a sufficiently 
small q, Gzi{q) can be approximated by 



V-^{uz{q}u,{-q))o = 



(5.21) 



Xz^qzqi 



(Bz{\ + 2/.) - \l^)(ql + f <zl)(gi + f gf) 



where 



/i = (A + 2/i) - 



B, ■ 



(5.22) 



B. Naive scaling and critical dimension 

To study fluctuations in the nonlinear elastic theory 
beyond the harmonic approximation of the previous sub- 
section, one might naively hope to perform a perturba- 
tive expansion in the nonlinear elastic terms. However, 
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a standard analysis, which is relegated to Appendix [Cl 
shows that such direct perturbation theory is hopelessly 
divergent at long length scales, as already suggested by 
divergent phonon fluctuations in e.g., Eq. (j5.13p . 

Thus, similarly to systems near a critical point (but 
here applied throughout the nematic phase), to treat 
elastic nonlinearities we need to employ the machin- 
ery of the renormalization group transformation (RG) 
[69l . [toI . [nj which establishes relations between physical 
quantities (e.g., correlation functions) at different length- 
scales. These relations then allow us to extract univer- 
sal long-scale nonperturbative (in nonlinearities) proper- 
ties of the system from their perturbatively computable 
short-scale versions. 

To this end, we need to study the property of the sys- 
tem under rescaling of length-scales and coarse-graining 
(thinning) of degrees of freedom. This procedure, ex- 
ecuted explicitly in Appendix [Cl is quite nontrivial in 
the presence of nonlinearities. However, it simplifies con- 
siderably to zeroth order in nonlinear terms, becoming 
equivalent to a rescaling transformation on the "bare" 
Hamiltonian. 

Applying this to the nematic elastomer model, 
Eq. (j4.31bp it is not difficult to see that it is invariant un- 
der each of the following two rescaling operations, with b 
an arbitrary rescaling factor: 

1. Rescaling of z axis: i?||(6) 

(i3,,A,i,A,Ai) = b^B',,X',^,X',^i'), (5.23b) 
K = b^K', (5.23c) 
(Q±,0||) = iQ'^,b-'Q\). (5.23d) 

2. Rescaling of x± plane: R±{b) 

{x±, z,u±,Uz) = {bx'j^, z' ,b^ u'j^,b'^ u'J, (5.24a) 

(i?„A,i,A,A*) = &-('^+3)(i?;,A;^,A',M'),(5.24b) 

K = b^-'^K', (5.24c) 

{Q±,Q\\) = (&-iQl,Qj,). (5.24d) 

We note that in above transformations, in addition to 
system's coordinates, phonon fields and elastic moduli, 
we have introduced two large wavevector cutoffs Q — 
{Q\\,Q±} beyond which our coarse-grained model is in- 
applicable. Here Qy and Q± are cutoffs in the direc- 
tions parallel and perpendicular to no = z, respectively, 
roughly set by the inverse of the mesh-size of the polymer 
network. They also provide an ultraviolet regularization 
to the path integral representation of the partition func- 
tion 

Z = J Due-"^'/^, (5.25) 

where Hei — J d'^'x Hei, with Hei given by Eq. (|4.31b[) 
and d'^x = dx'^^dz. Thus it is understood through- 
out that functional integrals (and wavevector mode sums 



that result from these) are over phonon fields u{x) whose 
Fourier amplitudes have support \q±\ < Q± and \qz\ < 
Q||. To simplify the subsequent analysis, it is convenient 
to choose a cylindrical cutoff, in which Qn = oo (i.e., qz 
is not cutoff at all) and to denote Q± simply by Q. Uni- 
versality principle guarantees that the long length-scale 
physics is independent of the choice of the short-scale 
cutoff. 

It is convenient to consider a special combination of 
these two rescaling operations R^^{b^) R±{b): 

{x±,z,u±,Uz) = (&f^,6^z',6"^M^, m'^), 

(5.26a) 

iBz,Xz±,X,^i,K) = b^-''{B',,X',^,X',fi',K'), 

(5.26b) 

(Q±,Q||) = ib-'Q'^,b-^Q\). (5.26c) 

that changes all elastic moduli by a common factor b'^^'^ 
and therefore preserves their relative ratios and the form 
of the Hamiltonian. Let c and c' be the shorthands for 
the sets of all elastic constants before and after rescaling, 
as shown in each side of Eq. (j5.26bp . Since Hgi is linear 
in the set of couplings c, it transforms according to 

Hei[u,c] = H,i[u',c'] = H,i[u'y-^c] = b''-^H,i[u',c]. 

(5.27) 

Substituting this result into Eq. (|5.25p . we find that such 
a rescaling transformation Eq. (j5.27p is equivalent to a 
rescaling of the temperature according to 

T' = yi-d^T. (5.28) 

Thus we find that for d < dc = 3 the effective tem- 
perature grows under rescaling, demonstrating increased 
importance of thermal fluctuations. In contrast, above 
the critical dimension dc = 3, fluctuation corrections to 
harmonic theory are small at low T . 

To simplify the notation, in subsequent calculations we 
rescale all elastic moduli by temperature T so that it does 
not explicitly appear in the functional integral, but can 
be easily restored by undoing this rescaling. 

C. Renormalization-group analysis 

To treat effects of nonlinearities beyond the zeroth 
order analysis of the previous subsection we employ 
the momentum-shell renormalization group, detailed in 
Appendix [C] To summarize, we coarse-grain the sys- 
tem by separating the phonon field into high- and low- 
wavevector components, as shown in Fig. 1101 

u^{r) = v?l{r) + u<{r), (5.29a) 
Uz{r) = u>{r) + u<{r), (5.29b) 

where and have support in the momentum shell 
Q e~*' < 9± < Q, while and uf have support in the 
inner cylinder Q < q^_ < Qe~^''. We then integrate out 
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qi plane 



FIG. 11: Feynman diagrams renormalizing B, C, X, fi, and 
K. 



FIG. 10: Basic scheme of momentum shell RG, defining the 
short (ti^ ^) and long (u^ ^) scale phonon fields. 



the high-wavevector parts , it , perturbatively in the 
anharmonic terms in H^i, Eq. (|4.3ip . thereby obtaining 
a coarse-grained elastic Hamiltonian in terms of fields 



and 



with all elastic constants renormalized by 



fluctuations of the ^ fields. 

As we have shown in the preceding section, the co- 
efficients of the nonlinear terms in the elastic Hamilto- 
nian Eq. (|4.3ip are completely determined by those of the 
quadratic terms, as summarized by the Ward identities, 
Eqs. and the form of Hei, Eq. ((OT|) . Since these 

identities are enforced by the underlying rotational sym- 
metry, they are preserved by the above coarse-graining 
procedure, i.e., the perturbative corrections to various 
elastic constant must satisfy the same Ward identities. 
This observation considerably simplifies our calculations, 
as it allows us to focus on the renormalization of the har- 
monic terms. These are summarized by the Feynman di- 
agram in Fig. [Til with correlators given by the harmonic 
theory, Eqs. (|5. 915. 1515. 19p . that for K3 = reduce to: 



Kql) 



X 



z± 



(A + 2m) (A?." 

1 (Bzql^ 



Kqi)qtqj 



(A 



1 



5 it 



qiqj 



Kq'i)qi 



(5.30a) 
(5.30b) 



(5.30c) 



with /t defined in Eq. (j5.12p . The underlying symmetry 
then enforces that the corrections to anharmonic terms 
are then completely determined by the Ward identities, 
as we have explicitly verified via detailed calculations pre- 
sented in Appendix [Cl 

In order to relate the resulting coarse-grained Hamil- 
tonian with the "bare" one, it is convenient to ap- 
ply a combination of two rescaling transformation 
R\\{e'^^^)R^{e^^), Eqs. (15. 2315.241) . to the field theory of 
low-wavevector fields u'^{x), so as to restore the ultravi- 
olet cutoff for q± back to Q. Here w is an arbitrary con- 
stant that can be chosen by convenience. The transfor- 
mation of the spatial coordinates, the phonon fields, and 



various elastic constants can be read off from Eq. ()5.23p 
and Eq. (IQij) : 



[xj_,z) = [e-x'j_,e'^^'-z'), (5.31a) 
u<{x) = e(3-2-)5'^^(f' )^ (5.31b) 

uf{x) = e(2— )*''u;(x'), (5.31c) 

(i?„A,^,A,M) - M-'^-^+^'^)''(B^,A;^,A',/i'), (5.31d) 
K = bi-d+i+^)si ^ (5 3^g) 

where the choice of the phonon field rescalings are dic- 
tated by the convenience of keeping the form of the non- 
linear strain tensor, w unchanged. 

The end result of these two operations (partial tracing 
and rescaling) is a nematic elastomer Hamiltonian, iden- 
tical in form to that in Eq. (j4.3ip . expressed in terms 
of phonon fields u'{x'), the original momentum cutoff Q, 
and effective elastic moduli {B^ + SBz,Xz± + SX2±, X + 
SX, /i -I- K + SK}. Assembling corrections to the elas- 
tic constants both from rescaling and from tracing out of 
short-scale fluctuations u'^ , we obtain their one- loop flow 
equations with the scale parameter I: 



dB, 
dl 
dX^± 
dl 



dX 
U 



d fj, 
U 
dK 
~dl 



id 
id 



3- 3w)B 



(d + 3-3cj)A 

'ipd 



3uj) X^± 

iBz - A,i) (X,i_ -X-fi) 



(Bz - A,j_)2(5.32a) 



(5.32b) 



id + 3- 



(5.32c) 
1 



3 w) /i 



— - {d-l-uj)K 



ipd 



(5.32d) 



(5.32e) 



8y/]nO iX + 2 fl) 
[B,{X + 2^) + 12^I{X + ^l)K -ACfi-C^] , 



where 



ipd 



2{2TT)d-i 



(5.33) 
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and Qd-i is the surface area of a d— 1-dimensional sphere. 
We note that (keeping consistent with a one-loop calcu- 
lation) in above we have set o? = 3 in the calculation of 
all diagrammatic corrections. 

The above flow equations simplify considerably when 
formulated in terms of parameters {B , C , fih , ■ Using 
Eqs. dSlSl), we find: 

dB 

u 

dC 

U 

dl 
d fi 

u 

dK 

~dT 

where various "q exponents encode diagrammatic correc- 



(d + 3- 


3 LO 


-Vb)B, 


(5.34a) 


(d + 3- 


3 LO 


-riL)C, 


(5.34b) 


(d + 3- 


3 LO 


- VL)Hh, 


(5.34c) 


(d + 3- 


3 LO 




(5.34d) 


(d - 1 - 


^ + 




(5.34e) 



tions and are given by 



VB = VL = '^9L, V±^'^9±, (5.35a) 



Vk = 



<?L {g± (8x - 24yV^ + 27) - 3<?l {y^ - l)) 



8 {9g±x + 9l{x- 6yy/x + 9)) 



(5.35b) 



In above, the gi,, g±, and two ratios x and y are all 
dimensionless and are given by 

5L = —^=^, g± = fT7^ ' (5.36a) 



B '■ 



c 



(5.36b) 



The flow equations for these four dimensionless cou- 
pling constants can be calculated from Eqs. (I5.32|) and 
Eqs. (|5.36all5.36b|) : 



<3L^^ gi ^ 

dl 16 (9g^a; + ffL (.T - 6j/V^ + 9)) (95L (y2 - 1) - g± (4a; + 12yV^ + 9)) 

[9g2(y - l)(y + 1) {-Ax + 24yVS + 9 (y^ - 5)) 

+2gLff_L (972\^?;3 _ 9(g4^ ^ 135)2;^ - 24V^(a; + 36)?; + 82;^ + 522x + 1377) 

+.gi (-27632;?;^ + 6V^(38a; + 153)y + 4a;(89x + 495) -I- 2106) ] , (5.37a) 

dgi_ , .g± 



dl 16(9g±a; + gL(.T-6j/V^ + 9))(g±(4a; + 12j/V^+9)-9gL(2/'-l)) 

[ - %\gl (2/2 _ 1) V \%glgAy - 1)(2/ + 1) (31x - lUy^/i + 144) 

^2 /'oncn^„.2 



-K5L.gi (3069a;?;2 + 114Vx(2a; - 9)y - 4(x(55a; 477) + 567)) - 18^1 a; {Ax + \2y^ + 9) ] > (5.37b) 



da; 
dl 



-4 5La;(l -2/2), 

f 1 /I 2\ 

-77 = -o 5Ly (1 - 2/ ), 



(5.38a) 
(5.38b) 



D. Solution of renormalization-group equations in 
three dimensions 



where e = 3 — d is the small parameter controlling the 
flow equations expansion, and as before the I dependence 
of the coupling constants is implicit. The initial condi- 
tions for the flow equations, (|5. 37115. 35)) are determined 
by the bare values of corresponding elastic moduli: 



5l(? = 0) 
x{l = 0) 



gx_{l = 0) = 



(5.39a) 
(5.39b) 



Naturally, we are most interested in three-dimensional 
nematic elastomers, corresponding to e = 0. Deferring 
the analysis of the d < 3 regime to Appendix IC 5[ we 
set e = in the right hand sides of Eqs. (|5.37all5.37bp . 
which leads to a vanishing of all linear (in gi^ and g±^ 
terms. Also, as we discussed earlier, rubber is nearly 
incompressible and is therefore characterized by a bare 
value of the bulk modulus B that is much larger than 
other moduli, C, ^l, and /i. Thus, using Eqs. (|5.36p we 
see that the initial values x{l = 0) and y{l — 0) are much 
less than one. Also, as we will prove later in this section, 
small x{l) and y{l) both flow to zero as / ^ 00. Conse- 
quently it becomes asymptotically exact to directly set x 
and y to zero in the flow equations for the coupling con- 
stants (7l(0 and g±{l), Eqs. (|5.37aj|5.37bp . which greatly 
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simplifies these equations: 



dl 
dl 



gi^ {5gl + 34g^gL + 26ffi) 

16(.9L + g±) 
gi_ {gl + 325i.9L + 28.gi) 



16(gL + 



(5.40a) 
(5.40b) 



To solve for these two equations, we define a new variable 
a{l): 



.9l(/) 
5±(0 



(5.41) 



From Eqs. (|5.40p we can then derive the flow equations 
for and tr: 



du 
H 
dgh 
dl 



gh 

= cr; 



4 

gl (20^2 



156(7+ 177) 



16 (4cr2 + 8(7 + 3) 



(5.42) 
(5.43) 



Anticipating that, for large I, ghil) asymptotically flows 
to zero as 



T 

9l{1) « y, 



(5.44) 



with T a positive constant (which we shall determine be- 
low), we easily see from Eq. (I5.42p that (t{1) also mono- 
tonically flows to zero: 



a{l)Kl- 



(5.45) 



Going back to Eq. (|5.4ip , we find that ct ^ implies that 
the ratio between the transverse and longitudinal shear 
moduli flows to a universal value 2: 



g±ii) i^{i) 



gL{l) fihil) 



2, as / 



oo. 



(5.46) 



Because of the slowness of a{l) decay (power-law in I 
and logarithmic in a length-scale) we expect that it will 
be difficult to observe this universal ratio in a real exper- 
iment. 

In the asymptotic regime where a^l) is small, we may 
simplify the flow equation for ^l, Eq. (|5.43p . by setting 
(7 to zero: 



dfL 
dl 



59 2 



which is solved by 



" 597- 



(5.47) 



(5.48) 



We have therefore verified the assumption (|5.44p . with 
T = 16/59. (5.49) 
0, we find the asymptotic 



Combining this with a{l) - 
flow of g±{l) to be 
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(5.50) 



We can now explicitly show that two dimensionless 
ratios x{l) and y{l) indeed both flow to zero for large 
I. Substituting the asymptotic solution Eq. (|5.48p into 
Eqs. (|5.38aj|5.38bp . we flnd: 



dx T 2n _ 

Ji = -47"^^-^)"-47"' 



dy T 



V 



T 

87^' 



with approximate asymptotic solutions 

y{i) « 2/or^/« 



(5.51a) 
(5.51b) 

(5.52a) 
(5.52b) 



indeed flowing to as anticipated. 

Armed with the results in Eqs. (|5. 4815. 5015. 52p . we 
then obtain the asymptotic behaviors of the rj{l) expo- 
nents as deflned in Eqs. (|5.35p : 



m(i) 



4^0 ,-(i+T/4) 

59 

''^(')=5^' 
38 
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(5.53a) 
(5.53b) 
(5.53c) 



Even though all 'q[l) exponents vanish as / oo, ?7_b(^) 
does so qualitatively faster than all others. As we shall 
discuss next, this leads to a qualitatively different scaling 
behavior of the renormalized bulk modulus as compared 
to other renormalized elastic moduli, and is responsi- 
ble for an exact asymptotic incompressibility of three- 
dimensional nematic elastomers. 



E. Anomalous elasticity of a thermal 
(homogeneous) nematic elastomer 

The renormalization group analysis of the previous 
subsection allows us now to relate the (difficult to com- 
pute) long-scale correlation functions, computed with 
bare elastic moduli to the corresponding functions at 
short-scales, but with renormalized elastic moduli. The 
latter are determined by the solution of the flow equa- 
tions (|5.32p . that are related to the microscopic moduli 
through the initial conditions. In contrast to the pure 
thermodynamics, which follows from the partition func- 
tion and thus only affected by transformation of the elas- 
tic moduli, correlation functions are also affected by the 
RG rescaling of the coordinates and the phonon fields. 
In reciprocal space, these transform according to: 



q^{l)e-\ 

e('^+2-)'^i(g(0), 



Kim), 



(5.54a) 
(5.54b) 
(5.54c) 
(5.54d) 



where at this point the anisotropy exponent lu is arbi- 
trary. 
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As an example, let us explicitly consider the RG trans- 
formation of the two-point correlation function of the 
phonon. Using the transformations of fields and wavevec- 
tors above, we find: 



where for later convenience we chose the arbitrary rescal- 
ing parameter w to be a; = 2. 

These flow equations are easily solved to give: 



(5.55) 

where again we have used c as a shorthand for all pa- 
rameters of the elastic Hamiltonian, that also enter the 
correlator Gzz (q, c) • 

Utilizing a multiplicative transformation of S'^{q) = 
HQz)^'''^~^H'1±) under rescaling (|5.54p we find: 



(C(/),ml(/),MO,^(0)-( 



C* III fi* 



19 T 



(5.56) 



Similarly analysis gives the transformations of other 
phonon correlators: 

GzMc) = e('^+4"2-)'G,,(<f(0,c(0), (5.57a) 
G,,(g,c) = e('^+=^^3'^)'G„(g(0,c(0), (5.57b) 

that then give the renormalized two-point vertex func- 
tions Tab, the inverse of Gab- 

TzMc) = e-('^+3-")'r,,(g(0,c(0), (5.58a) 
r,,(g,c) = e-('^+4-2.)'r,,(9(/),c(/)), (5.58b) 
r,,(g,c) = e-('^+5-3.)'r,,(g'(0,c(0). (5.58c) 



Now we may use the asymptotic solutions Eqs. (I5.53P 
for r\ exponents in the flow equations (|5.34p that control 
the large I behavior of ciX): 



dB 
~dl 



59 ;(i+Tr/4)' 



(5.59a) 



^(G,ml,m) - (G,/iL,/i), (5.59b) 



dK 

m 
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K. 



(5.59c) 



;T/4 ' fC/A ' IT /A ■ 

(5.60) 

where /i* = 2/i£ due to Eq. (|5.46[) . The slow decay of 
the effective G, /iL, /i and growth of K with length-scale 
contrasts with the long-scale behavior of the bulk mod- 
ulus that flows to a nonzero and nonuniversal value B* 
determined by yo as well as the initial (bare) value of B: 



B{1) 



B* 



as I 



(5.61) 



We can now use these relations to determine the 
wavevector dependence of these renormalized vertex 
functions by choosing the flow parameter / such that 



9±(0 
9.(0 



q_i_ e 



qz e'^ 



Qz — - 

<11 



(5.62a) 
(5.62b) 



This choice ensures that the rescaled vertex functions on 
the right hand side of Eqs. (|5.58p are easy to compute as 
they are evaluated at the large wavevector Q and with 
small nonlinear couplings gL,j_(0; as latter flow to for 
large l{q±) ~ logQ/gj^ (small q). Namely, in this regime 
the effect of anharmonic fluctuations is negligible and 
we can replace the right hand sides of Eqs. (|5.58p by 
harmonic vertex functions. Focusing on Vzzisi-i'^ as an 
explicit example, we have: 



r,,(g,c) - e-4'r,,(q(Z),c(Z)) = e-4' 

„-4/ 



Bz{l)qz{lf + K{l)q^{l) 

4 



(5.63) 



where we have used Eq. (|2.45p (with d = 3) to express Bz in terms of B, C, and /il. Using Eqs. (|5. 6015. 6115. 62p . we 
finally find the renormalized vertex function Tzz {q, c) to be given by 



B* 



4 4 
B* + {-C* + -iil) 







log — 




q± 





ql ql^Q^ + K* 
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log — 









■G* 



log- 



ql + Kl 



log- 



9-L 



(5.64) 
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The other two vertex functions can be similarly calculated, 



B* 



B* 



1 ^* 2 , 



2^* 1 * 

--C + - Ur 

3 9^^ 



1 Q 

log — 



1 Q 

log 


4 ■ 
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mj + 




1 ^ 
log — 
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9J- 






9-L 





r 



(5.65) 



where we have restored in Fy . Comparing these 
renormalized vertex functions with their bare forms, 
Eas. l|5.5p , we can interpret the effects of long wavelength 
elastic fluctuations as wavevector- dependent renormal- 
ized elastic moduli: 



B{q) - B*, 
log- 



Ki{q) cx 



log — 



(5.66a) 
(5.66b) 

(5.66c) 



We observe that the renormalized shear moduli /i and 
/iL are singular functions of the wavevector and vanish 
in the long wavelength limit. In contrast, the renormal- 
ized bulk modulus remains finite, while the renormalized 
splay constant diverges in the long wavelength limit. As 
advertised, these results show that in the thermodynamic 
limit an ideal nematic elastomer is effectively strictly in- 
compressible. They also imply an absence of a linear 
stress-strain response (i.e., a breakdown of Hooke's law) 
even for a nonsoft shear deformation, as (due to a vanish- 
ing of iJ.{q) at long wavelengths) the shear stress vanishes 
faster than linearly with a vanishing strain. Finally due 
to Eq. (|5.46[) , the ratio of two renormalized shear moduli 
goes to a universal number: 



lim — — - = lim 



9l{1) 



= 2. 



(5.67) 



The above anomalous behavior begins to manifest it- 
self roughly when the asymptotic values of coupling con- 
stants (as given by Eq. (|5.48p and Eq. (|5.50p ) become 
comparable with their bare values, given by Eqs. (|5.39p . 
This happens at ^ = T, which satisfies 
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(5.68) 



where we restored the temperature dependence. This 
value of I* in turn determines the two crossover length- 
scales and via 



which are thus given by 

C-L = Q^^ cxp 
^11 = exp 



59^<i fii^ksT j ' 



(5.70a) 
(5.70b) 



.2r 



(5.69) 



with fi defined by Eq. (|5.14p . 

It is interesting to estimate the order of magnitude of 
these two important crossover length-scales and ^y. 
In three dimensions, ipd = l/47r. The typical value of 
Q^^ (mesh-size of the polymer network) is a few nanome- 
ters. At room temperature, ksT is roughly 4 x 10~^^ J. 
The splay constant Ki is usually 2 — 4 x 10~^'^N for 
low molecular weight nematic liquid crystals, but may 
be several times larger than this for liquid crystalline 
polymers. The shear moduli fi and /iL may vary signif- 
icantly near the vulcanization transition, with a typical 
range of 10'*Pa— lO^Pa. Substituting these numbers into 
Eq. (|5.70p . we find a wild range of these two crossover 
length-scales: 

IQ-^m < C_L < lO^m, lO^^m < ^|| < lO^^m, (5.71) 

where the lower limit corresponds to large shear moduli. 
Therefore for soft elastomers, the anomalous effects of 
thermal fluctuations may be impossible to observe, while 
for hard elastomers, thermal fiuctuations may become 
important even at the scale of the polymer network mesh 
size. The ratio between two length-scales and also 
exhibit a wide range and depends sensitively on the shear 
moduli. This extreme sensitivity of and to shear 
moduli is of course due to the exponential functional form 
in Eqs. (|5.70p . which is a consequence of the marginal 
irrelevance of temperature in three dimensions. 



VI. LONG-SCALE ELASTICITY OF A 
HETEROGENEOUS NEMATIC ELASTOMER 

A. Elastomer heterogeneity 

In all preceding analyses we have treated nematic elas- 
tomers as homogeneous continuous elastic media. How- 
ever, as all rubber and gels, nematic elastomers are ran- 
dom polymer networks that are only statistically homo- 
geneous and isotropic at the macroscopic level. At the 
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microscopical level, their elasticity is due to a heteroge- 
neous polymer network, made by random crosslinking of 
a polymer melt, a process called "vulcanization" or "gela- 
tion" . Thus, for a complete understanding of the elastic- 
ity of nematic elastomers, it is essential to examine the 
nature and consequences of the network heterogeneity. 

Experimental evidence unambiguously shows the 
importance of network heterogeneity. Elastomers 
crosslinked in the isotropic phase, upon lowering temper- 
ature below the I-N transition point of the uncrosslinked 
system, are observed to exhibit a poly-domain "nematic 
phase" , where a nematic order freezes locally into a 
micron-size domain structure, with no long-range orien- 
tational order. When the system is uniaxially stretched, 
these randomly orientated nematic domains align as the 
strain deformation exceeds a threshold value [T^l . When 
the strain is removed, the original poly-domain pattern 
is restored [t^. On the other hand, nematic elastomers 
crosslinked under anisotropic condition, i.e., inside the 
nematic phase or under a uniaxial stretch, do exhibit a 
long-range nematic order, but are characterized by semi- 
soft elasticity, with a small but finite shear modulus fiz± 
[21I [2^ . When heated back to high temperature, the 
system exhibits a para-nematic phase with a weak rem- 
nant nematic order. Thus, such network permanently 
breaks the rotational symmetry with the isotropic phase 
no longer accessible. That is, the system "remembers" 
the anisotropic crosslinking conditions, i.e., its history 
of formation. The mechanism underlying this memory 
effect is poorly understood. 

To understand how the history of formation affects 
the properties of nematic elastomers (or more generally 
speaking, vulcanized solids), a detailed characterization 
of elastic heterogeneities is needed. A semi-microscopic 
level of understanding of network heterogeneity may be 
obtained from the vulcanization theory. As pointed out 
by de Gennes (t^I and Edwards (t^, elastomers (gels) 
are "frozen" random systems. Their physical properties 
depend on both the state within which the measurement 
is taken, and the preparation state of the system. A com- 
plete study of these systems necessarily involves two en- 
sembles: a preparation ensemble and a measurement en- 
semble. The vulcanization theory [zl, [t^ captures both 
of these and therefore has the potential to yield a bet- 
ter statistical characterization of elastic heterogeneities 
in elastomeric materials, as well as their connections to 
the crosslinking process and elastic properties of the re- 
sulting amorphous solid. One recent theoretical work [t^I 
along this direction has indeed showed some success. At 
this stage, however, it is beyond the Landau approach 
that we take in this work and will therefore not be dis- 
cussed further. 

Within the framework of continuous elasticity theory, 
two types of quenched disorders can be identified: ran- 
dom elastic constants and random internal, or residual 
stress. Internal stress is a stress configuration of solids 
that does not vanish even if all external tractions are re- 
moved. It is a generic property of random solids and sys- 



tems with topological defects. At the equilibrium state, 
the internal stress has to satisfy 

= 0, (6.1) 

which is essentially the condition for mechanical equi- 
librium. This implies that random internal stress het- 
erogeneity only couples to the linearized strain at the 
quadratic order: 

J d'^'x aljc^j = J df^x a\j-{diUj + djUi + diu ■ dju) 

^ j d''x]^a%dS-djU, (6.2) 

where is the nonlinear Lagrange strain tensor defined 
relative to the true equilibrium state. 

DiDonna and Lubensky [78^ have recently studied the 
effects of random internal stress and random elastic con- 
stants in macroscopically deformed conventional solids. 
They found that random elastic constants, but not ran- 
dom internal stress, induce nonaffine components of the 
deformation field. These authors have also introduced 
short-range correlated random stress field into an oth- 
erwise homogeneous solid, by, e.g., starting from a tri- 
angular lattice of particles connected by central force 
springs, and randomly change the equilibrium length of 
each spring. They found that, after re-expanding the 
elastic energy around the true equilibrium state, the sys- 
tem exhibited both random elastic constants and random 
internal stress. 

In this paper we shall take a phenomcnological ap- 
proach similar to that of DiDonna and Lubensky [78|. 
We start from an ideal homogeneous nematic elastomer 
and introduce short-range (uncorrelated) random tensor 
fields that couple to strain and nematic order parame- 
ter. By re-expanding the elastic Hamiltonian around the 
ideal nematic reference state (i.e., around the true ground 
state for the disorder-/ree system) and integrating out 
the fluctuations of the nematic director field, we obtain 
an effective strain-only model for heterogeneous nematic 
elastomers, where the nonlinear Lagrange strain tensor 
is coupled to a random symmetric tensor field, which we 
shall refer to as the random initial stress. Such a stress 
tensor field, <t generically contains both the longitudi- 
nal part (satisfying V x cr^ = 0) and a transverse part 
(satisfying V • cr^ = 0). However, since we shall only 
consider the linear coupling between the random initial 
stress and the linearized strain tensor (see Eq. (|6.2[) ). our 
model ultimately only captures the fluctuations in the 
longitudinal part of the initial stress, but not the trans- 
verse one. Although we might argue that the latter is 
less important as it only couples to the nonlinear part of 
the strain, a more detailed analysis and justification of 
this approximation is necessary, but will not be studied 
here. 

Our analysis below shows that these initial random 
stress fluctuations qualitatively modify macroscopic elas- 
tic properties of heterogeneous nematic elastomers, lead- 
ing to a zero-temperature analog of anomalous elasticity 
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discussed in the previous sections for thermal homoge- 
neous elastomers. In fact, we will find that these random 
initial stress effects strongly dominate over purely ther- 
mal fluctuations and already below five (and therefore 
in the physically most interesting three) dimensions lead 
to shear moduli that vanish as power-laws of a prob- 
ing wavevector, and to a breakdown of Hooke's law, re- 
placed by a strictly nonlinear stress-strain relation for 
arbitrarily small strain deformation. These anomalous 
elastic properties are consequences of the subtle inter- 
play between the elastic heterogeneities and the long 
wavelength fluctuations of soft deformations. Neverthe- 
less, we find the the long-range nematic order is stable 
against these random initial stress fluctuations in three 
dimensions, at least within the approximation of ignoring 
transverse internal stress fluctuations, and within the ap- 
proximation of one-loop renormalization group analysis. 
We emphasize that our analysis of the ordered state pro- 
vides a necessary stability condition for the existence of a 
mono-domain nematic elastomer in the presence of weak 
quenched disorder. Finally, our way of modeling nematic 
elastomer heterogeneity, though simple, is far from well 
justified. A justification of this approach must come from 
a more detailed heterogeneous model of rubber elasticity 
and ultimately from comparison with experiments. 

In the following sections we present a detailed analysis 
of the fluctuations of the initial stress in nematic elas- 
tomers and their effects on the long-scale elastic prop- 
erties of the systems at low temperature. A summary 
of these results have already appeared in Ref. [1^ . We 
leave other important and challenging questions such as 
for example the nature of the isotropic- nematic transition 
in heterogeneous elastomers (i.e., whether it survives or 
is replaced by a crossover) to another publication |79i] . 



B. Heterogeneous strain-only elastic model 

We are interested in nematic elastomers crosslinked un- 
der isotropic conditions, i.e., inside the isotropic phase 
and in the absence of shear strain. Under these con- 
ditions, the lowest order elastic energy terms that in- 
volve network heterogeneities, elastic deformations, and 
nematic order, and at the same time satisfy all relevant 
symmetries are: 



H^l = - I (TX 



Tr^Fi(X)A'^A^ 



(6.3) 



-^Tr( F2(A)A^Q(X)Aj 



where, as defined in SecHH X is the isotropic referential 
coordinate and A the deformation gradient defined rela- 
tive to the isotropic reference state. Fi(X) and F2(A^), 
encoding network heterogeneity, are two random tensor 
fields in the isotropic reference space. They transform as 
scalars with respect to rotations in the embedding space. 



in agreement with the fact that they describe the in- 
trinsic properties of random network. Since the nematic 
order parameter Q is a tensor in the embedding space, 
it can couple to the tensor field F2 only through the de- 
formation gradient A, hence the structure of the second 
term in Eq. (|6.3p . Since the system is isotropic at the 
macroscopic level, the spatial correlation of the tensor 
fields Fi(X) and F2(X) must be isotropic. To simplify 
the analysis, we shall also assume that they are Gaussian 
and short-range correlated. We may also say something 
about the microscopic origin of these two terms. Fi(af), 
directly coupled to the Lagrange strain tensor, may be 
simply due to the fluctuations of crosslink density. F2(af), 
coupled to the nematic order through the deformation 
gradient, may be due to the orientational effects of rod- 
like crosslinkers on local nematic order, as is illustrated 
in Fig. [TH 
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FIG. 12: A cartoon for a random tilting field coupled to a ne- 
matic order of a nematic elastomer. In (A), a crosslink (the 
black rod) connects two polymer chains. Neighboring liquid 
crystalline mesogenic units tend to align along this rod-like 
crosslink, an effect that can be modeled as a random tilt- 
ing field coupling to nematic order (shown as blue ellipses). 
In (B), this random tilting field rotates with the polymer 
network. Therefore elastic deformation of the network also 
changes this random tilting field seen by local nematic order. 
This property is captured by Eq. (|6.3p . 



For a given realization of quenched disorder fields 
Fi(A) and F2(A), in contrast to the homogeneous case 
of Sec. ini the isotropic reference state (IRS) defined by 

r{X) = X, Q = 

is generically not the true ground state, i.e., it does not 
minimize the total heterogeneous elastic energy, which 
is the sum of Eq. (|6.3p and the elastic energy for ideal 
nematic elastomers Eq. (12.431) . This raises the question 
of physical significance of the IRS in a realistic hetero- 
geneous elastomer. Indeed IRS has no physical signif- 
icance; it just provides a conceptually convenient but 
arbitrary reference point around which we expand the 
elastomers free energy. In such an expansion, the low- 
est order two terms involving random heterogeneity are 
shown in Eq. jO]) 
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The true ground state of a disordered nematic elas- 
tomer can always be found by minimizing the total elas- 
tic energy. Hence for a given realization of the disorder, 
there is a well-defined relation between the IRS and the 
true ground state. One could in principle re-expand the 
elastic energy around the true ground state. However, 
due to the nonlinear nature of the elastic system, this 
analysis is rather messy and the final result is not illumi- 
nating. 

We re-expand the disorder Hamiltonian in Eq. (j6.3p 
around the ideal nematic reference state (NRS), which 
is the ground state for a disorder-free system as defined 
by Eq. ^M^, and is related to the IRS by Eq. (j^ . 
As usual we choose the (uniform) nematic order in the 
NRS to be along z axis so that ng = z in Eq. (|2.6p . As 
in Sec. |TT1 we shall use the notation x for the position 
of mass point in the ideal nematic reference state, and 
define the "phonon field" u{x) to be the displacement 
from this NRS, i.e., 



- u{x). 



(6.4) 



To further simplify the analysis, we integrate out the 
fluctuations of nematic director 5h around the ideal ne- 
matic reference state. After a tedious but conceptually 
straightforward calculation, which we relegate to Ap- 
pendix [Dl we find that, to a linear order in the phonon 
field u, the most relevant part of the disorder Hamilto- 
nian Eq. (|6.3p reduces to 



d'''xaab{x)sab{x), 



where 



£ab{x) = {daUb{x) + dbUa{x)) . 



(6.5) 



(6.6) 



are the components of the linearized symmetric strain 
tensor, defined relative to the ideal nematic reference 
state. The quenched random tensor field (t{x) as a func- 
tional of Fi and F2 is given by Eqs. (|D10|) . Since it is 
coupled to the linearized strain, we shall call it the ran- 
dom initial stress, to make a distinction with an internal 
stress field in the equilibrium state which satisfies the 
mechanical equilibrium condition daCfab ~ ^- Obviously, 
as mentioned above, only the longitudinal part of the 
random initial stress (Jab{x) contributes to Eq. (|6.5p . 

It is further shown in Appendix [D] that, inside 
Eq. (j6.5p . the most relevant terms are aizSiz^ where i 
and j only take values x and y. Other terms, such as 
OijSij and a^z^zz, are less relevant (in the RG sense) and 
therefore can and will be neglected. The statistics of 
the disorder fields Uzi (x) is assumed to be Gaussian with 
short-range correlations: 



^{x)aziix') ^ AAS.j S'^ix - x'). 



(6.7) 



In the framework of the strain-only description, the 
total elastic Hamiltonian is then given by Eq. (|4.3f ap 



for the pure system, augmented by the leading disorder 
terms in Eq. (j6.5p : 



H 



d^x 



^Bzwl^ + \zl_WzzWu + ^Aw| -I- fiWijWtj 



+ ^(Vi«.)' + ^(dlu^f - <yz.{S)ez,{x) 



, (6.8) 



where the effective strain tensor w is defined in 
Eqs. 

Since the Hamiltonian (j6.8p contains random stress 
terms linear in Siz , it is clear that the NRS with r = a' is 
not the true ground state. Rather the latter is given by 



fo(f, cr) = x + uo{x,a), 



(6.9) 



where uq (x) describes the deformation of the ground rel- 
ative to the ideal nematic reference state due to the pres- 
ence of quenched disorder (random stresses). The physi- 
cal phonon field, which we denote by df{x), is defined as 
the displacement from the true ground state ^811]. It is 
related to the fictitious "phonon field" u{x) by 

r{x) — fo{x,a) + Sr{x) ^ x + u{x), (6.10) 
Sr{x) = u{x) — uq{x , a) , (6-ff) 

The harmonic thermal fluctuations around the 
disorder-dependent ground state fo (x, tr) can be char- 
acterized by the disorder averaged thermal correlation 
function of the physical phonon field 5r(x), 



G'{x-y) = {5r{xW{y)) = {u - {u)){u - {u)),{Q. 12) 

where we use angular brackets and over-bar to denote 
thermal and disorder averages, respectively. Note that 
G""" is a tensor with components G^f^. On the other hand, 
the quenched fluctuations, i.e., sample-to-sample varia- 
tions of the ground state relative to the ideal NRS can 
be characterized by the quenched correlator 



G {x~y) = UQ{x,a)uo{y,a) = {u){u), (6.13) 

where is also a tensor with components G^^. As 
shown in Appendix [El G"^ and G^ can be calculated 
using the replica method [s^ . 

Even though the choice of the ideal isotropic and ne- 
matic reference states are rather arbitrary, the distortion 
field uo{x) has important physical significance, as it is 
closely related to the local nonafhne deformation field for 
a macroscopically strained elastomer. To see this, let us 
consider applying a macroscopic deformation A, with A a 
constant matrix. Because of the network heterogeneities, 
the elastomer does not deform affinely as in the idealized 
homogeneous (disorder-free) case. The new ground state 
configuration subject to the macroscopic strain can be 
parameterized by 



ro(f , a,A)—A-x + uo{x, a, A), 



(6.14) 
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where we have shown exphcitly the dependence of the 
ground state on A. 1x0(2?, cr, A) can be found by mini- 
mizing the total elastic energy Eq. (|6.8p subject to the 
constraint of macroscopic deformation A. The local non- 
affine part of the deformation field can then be charac- 
terized by: 



t(x, (7, A) 



fo(f, cr, A) 
uo{x, (T, A) 



A • ro{x, a) 

- A • uo{x, fj), (6.15) 



where we have used Eq. (|6.9p and Eq. (|6.14p in the last 
equality. t(x, cr, A) is generically nonzero for a system 
with quenched disorder. The correlation function of the 
nonaffine deformation field t can also be calculated us- 
ing the replica method and is found fW\ to be linearly 
proportional to the quenched correlator G'^: 



t{x,a)t{y,a) (X uoix,a)uo{y,a) = G (x-y). (6.16) 

As we will demonstrate shortly, the renormalized 
quenched correlator (and hence the nonaffinity corre- 
lation function) will exhibit scaling that is distinct from 
that of the thermal correlator G^J. Physically, this is 
a reflection of the strong sensitivity of the nonaffine re- 
sponse of the ground state to an external traction. The 
resulting macroscopic strain deformation is a nonlinear 
function of the stress even for an infinitesimal value of 
the stress. 



C. Replica trick and harmonic theory 

The fields azi{x) are quenched random variables, i.e., 
their values are fixed for a particular sample by the net- 
work heterogeneities, and in contrast to the phonon and 
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nematic director fields do not fluctuate on experimental 
time scales. We make a standard assumption that the 
system is self-averaging, namely, that physical quantities 
for a typical system coincide with corresponding disorder 
averaged quantities. 

The disorder average can be performed using the stan- 
dard replica trick, which we review in Appendix [E] As 
explained there, for an arbitrary nonzero integer n, we 



define an n-replicated Hamiltonian iJ„ 



by 



exp (~Hn[u^ , 



]/T) = nexp(-i?P"]/T) 



exp I - ^ H[u"] j (6.17) 



\ a = l / 

To compute physical quantities, at the end of their cal- 
culation we will analytically continue all replicated n- 
dependent quantities from integer n to positive real val- 
ues < n < 1, and will take the replica limit n ^ 0. 

To calculate the disorder average in Eq. (|6.17p . we use 
Gaussian, zero-mean field identity 



together with Eqs. (|6.5l6.7p . we find 



(6.18) 



al3 



Consequently, we find that the replicated Hamiltonian 
is given by 



I (Trw")2 + C (Trw") w^^ + ^ + ^^g. 



iifi(Vi<)2 + ii^3(afui)^ 



2A 



f d'xY^el 

a,/3 



(6.20) 



where a and /3 are replica indices summed over 1, . . . , n. 

We first study the harmonic part of the replicated elastic Hamiltonian Eq. (|6.20p . which is given by: 



a, 13 



where Tio is defined in Eq. (|5.2p . In Fourier space, it can be written as 



(6.21) 



(6.22) 



a/3 ah 



where indices a and h are summed over all d possible values x, y, . . . , z, Tahio) is given by Eq. (j5.5p . and J is a 
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d X d matrix with components Jab- 



Jij 




(6.23) 


J iz 


= Jzi = ^ Ag^gi, 


(6.24) 


Jzz 




(6.25) 



In replica space one can think of J as either an n x n 
matrix with all elements equal to 1 or simply as a scalar 
equal to 1. 

The harmonic correlation functions (propagators) of 
the replicated model can be easily calculated: 



1 

v' 



'<{q) uU-q))o = GlM + Gi(q), 

(6.26) 

where the thermal correlators G^j, are defined by 
Eq. and are given by Eqs. (|5.9l5.15l5.19p . The har- 
monic thermal correlators are linear in temperature and 
independent of random stress variance A. The harmonic 
quenched correlators are given by 



Gf^, = T-i(GJG)„ 



6 ' 



(6.27) 



and are linear in the disorder variance A and independent 
of T. It is easy to check that in the n ^ limit G"j^{q) 
is indeed the inverse matrix of the kernel in Eq. I|6.22p : 

E E '^ab [^bc{q) - JbM] - T Sac Sc.,. (6.28) 

/3 b 

Let us first look at the quenched fluctuations of the 
phonon around the ideal nematic reference state in real 
space: 



{u.{f))l = G?/^(f=0)= / f^iGfM 



J ^!^A[q,GUq)+q,G.AQ)? (6.29a) 

d'^q Aql 
{2nr ifiq^ + Kqi)^ ^^'^^^^ 

^ i^^-"*), ford<5, (6.29c) 



where in going from Eq. (|6.29ap to Eq. (|6.29bp . we have 
only kept the most (infra-)divergent term, proportional 
to Gl^. The result in Eq. (|6.29cp shows that quenched 
fluctuations of field diverge with the system size for 
d < 5, suggesting a break down of the harmonic theory 
at long length-scale. Therefore our model of a heteroge- 
neous nematic elastomer, Eq. (|6.20p has an upper-critical 
dimension d^ = 5. This should be contrasted with a ho- 
mogeneous thermal elastomer that we studied in Sec. IVl 
where the critical dimension is d^ = 3. 

On the other hand, thermal fluctuations of the phonon 
Uz about the ground state for a specific realization of 



quenched disorder are given by: 

d'^q 



(2 7r)' 
d<^q 



T 



Gf.m 



{2T:Y[iql 
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(6.30) 



and as found earlier are far weaker, finite for d > 3. 
However, as we will see shortly, thermal fluctuations of Uz 
around the disorder-renormalized state are in fact finite 
at and even below three dimension. 

Similarly, we can show that the most divergent part of 
the quenched fluctuations of u± about the ideal NRS are 
given by 



GtHr-- 



0) 



d-^q 
(2 7r)' 



■GtM 



= A 



d'^q Aq\ 
d'^q 



[Gl + {d~l)G^ 

„2 



{d-l)ql 



{jlql + K^qtY (/.gi + if3gf)2_ 

(6.31) 



where, again, we have only kept the most infra-divergent 
terms. One can readily see that this integral diverges 
with the system size only when d < 7/2, similarly to a 
randomly pinned columnar liquid crystal and vortex lat- 
tice in a disordered magnetic superconductor [H, [s^ . As 
in the case of thermal fluctuations of a homogeneous elas- 
tomer, we again find that the fluctuations of Uz phonon 
field are qualitatively more important than those of u±. 
The difference in critical dimensions for these two phonon 
fields justifies our model in Eq. ()6.8p . where we only kept 
the more relevant smectic (as opposed to columnar) non- 
linearities. 



D. Renormalization group analysis 

We will therefore focus on the larger quenched phonon 
Uz fluctuations. To simplify the analysis, again we drop 
the dangerously irrelevant term K3{dzU±Y in ths elas- 
tic Hamiltonian, Eq. (|6.8p . Furthermore, we will neglect 
the coupling between the quenched random stress fleld 
azi and the u± phonon fields, as we have just seen that 
fiuctuations of the latter are subdominant to those of the 
Uz field. These simplifications amount to making the fol- 
lowing replacements in the elastic Hamiltonian (|6.20p : 



0, 



a 
£■ E- 
^iz^iz 



i(V^<).(V^j.f). 



(6.32a) 
(6.32b) 



Similar to Sec. El we rescale elastic constants B, G, 
/iL, 1^ and K = Ki by temperature T, and at the same 
time also scale the disorder variance A by T^, such that 
T does not appear in the formalism. We shall restore the 
T dependence after the technical RG analysis. 
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To study the effects of tfie network lieterogeneity and 
the elastic nonlinearities beyond harmonic approxima- 
tion, we perform a momentum-shell RG analysis. We 
use the same cylindrical ultraviolet cutoff scheme {Qi_ — 
Q,Q\\ = oo) and rescaling transformations Eqs. (j5.3ip as 
in Sec. |Vl After ignoring the irrelevant terms according 
to Eqs. (|6.32p . the harmonic correlators of the replicated 
elastic Hamiltonian, Eq. (|6.26p reduce to: 



G,, (q) {S^p + Aql C iq)) , (6.33a) 
G,^{q) {5^p + A gi G,M) , (6.33b) 
G,M + Agi C{q) G,,iq}, (6.33c) 



with Gabi<f) the harmonic thermal correlators given by 
Eqs. ([5301) . 

Detailed calculations of the diagrammatic corrections 
to various model parameters {B^, Azj_, A, A} are 
quite involved and are relegated to Appendix |F] Their 
flow equations are listed in Eqs. (jF13p . Similarly to the 
case of thermal fluctuations, we find that the RG fiows 
of these parameters are controlled by four dimensionless 
coupling constants, defined as follows: 

gL = —F^=, g-L = -^=^, (6.34a) 
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c 



y = 



(6.34b) 



where ipd is defined in Eq. (|5.33p . In contrast with 
the case of thermal fluctuations discussed in Sec. |V] 
(Eqs. (|5.36p ). however, here the two coupling constants 
gL and gj_ are proportional to the disorder variance. A, 
instead of temperature, T. The flow equations of all elas- 
tic constants are given by: 

dB 

HJ 

dC 
dl 

dl 
d/i 

dl 
dK 

IT 
dA 

'dl 

where anomalous 77 exponents are functions of dimen- 
sionless coupling constants gL, g_L, x, and y, given in 
Eqs. (IF22IF23I) of Appendix [H 

As discussed in the context of homogeneous elastomers, 
the bare values of two dimensionless ratios x and y are 
much smaller than one. Furthermore, as we shall show 
in Appendix [f1 they flow to zero exponentially as long 
as d < 5. Therefore to a good approximation, which be- 
comes asymptotically exact, we can set them to zero in 
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all flow equations, thereby considerably simplifying cal- 
culations. In this limit, the 77 exponents reduce to: 



(6.36a) 
(6.36b) 
(6.36c) 
(6.36d) 
(6.36e) 



The RG flows of the dimensionless coupling constants, 
gL, g_L, X and y are also calculated in Appendix iFl and 
in the limit that both x and y approach zero reduce to 
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dgL 
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5 gL (4 gL^ + 44 gLg_L + 51 g_L^) 

32(2gL + 3gj.) 



(6.37a) 



g± (-4 gL^ + 188 gLg± + 261 gi^ 



32(2gL + 3gi) 



(6.37b) 



where e = 5 — d. 

Below five dimensions, e > 0, Eqs. (|6.37[) admit four 
fixed points, Gaussian (G), Smectic (S), X, and Elas- 
tomer (E), which we list in Table IIIII Also shown in 
Table Hill are the 77 exponents for all elastic constants, as 
defined in Eqs. (j6.36[) . The fiow pattern of gL and gj_ 
under a renormalization group transformation for d < 5 
is shown in Fig. [T51 

It is particularly interesting to note that at fixed point 
S, where gj_ vanishes, various 77 exponents are identi- 
cal to those of a smectic liquid crystal confined in ran- 
dom environment of e.g., an aerogel matrix, discovered 
by Radzihovsky and Toner J4] . This is not merely a 
coincidence. It is clear from Eq. (|6.8p that if the trans- 
verse shear modulus = vanishes (that is, g_L), then 
u± phonon fields can be integrated out completely. The 
resulting effective model is identical to that of a smectic 
liquid crystal with one-dimensional phonon field Uz and a 
shifted compressional modulus, as well as a random tilt- 
ing field with variance A. This is precisely the model ex- 
tensively studied by Radzihovsky and Toner jsi], [13, HI] , 
and a recovery of their randomly-pinned smectic flows 
and exponents here is a nontrivial check on our calcula- 
tions. 

We are not aware of any physical system that is de- 
scribed by the fixed point X, characterized by gL = 0. 
We note that rjB vanishes at all four fixed points. Conse- 
quently the bulk modulus B does not acquire any anoma- 
lous dimension, with the underlying physical reason for 
this already discussed in Sec. |Vl 

As shown in Fig. [T3] only the fixed point E is stable in 
both gL and gj^ directions. It therefore controls the long 
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length-scale elasticity of a heterogeneous nematic elas- 
tomer, which, based on this fixed point shall be analyzed 
in detail in the next section. 
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TABLE III: Fixed point couplings and rj exponents for het- 
erogeneous nematic elastomer, x* = j/* = at each of these 
fixed points. 




FIG. 13: Renormalization group flow diagram for dimension- 
less couplings gL and of a heterogeneous nematic elas- 
tomer. Gaussian fixed point G of a harmonic elastomer is 
unstable to elastic nonlinearities, flowing to a globally stable 
fixed point E that controls long-scale properties of a heteroge- 
neous nematic elastomer. The fixed point S at gx is identical 
to that of a randomly pinned smectic studied by Radzihovsky 
and Toner jsH,!!!]. 
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where, to facilitate later discussion on renormalized 
stress-strain relation, we have also introduced an exter- 
nal nominal stress Sia coupled linearly to the deforma- 
tion gradient, which we first discussed in Sec. IIVI (c.f. 



VII. ANOMALOUS ELASTICITY OF 
HETEROGENEOUS NEMATIC ELASTOMERS 

As we discussed in previous section (and is true more 
generally), fluctuations associated with local network 
heterogeneity are dominant over thermal fluctuations. 
This is also supported by the difference between the crit- 
ical dimensions of the homogeneous and heterogeneous 
nematic elastomers. We further note that the coupling 
constants gL and g_L defined in Eqs. ()6.34|) depend on 
quenched disorder A but not on temperature T. On the 
other hand the thermal coupling constants and g±, 
whose flow we studied in Sec. |V] for homogeneous case, 
are less relevant. In fact we shall see below that the 
fixed points we have identified in the preceding section 
are finite-disorder, zero-temperature fixed points. Qual- 
itatively speaking, the universal, long-scale elastic prop- 
erties of a heterogeneous nematic elastomer are deter- 
mined by its ground state, which sensitively depends on 
the particular realization of disorders and external trac- 
tion, i.e., macroscopic strain deformation. Thermal fluc- 
tuations around the ground state remain finite and are 
qualitatively unimportant. Similar scenario also appears 
in the equilibrium physics of random field Ising model 
and smectic liquid crystal confined in random geometry. 



A. Zero-temperature fixed point 

To illustrate the physics of the new fixed point E that 
we found in the preceding section, let us first summarize 
the RG analysis for the disordered model in a slightly 
different notations. Let us first restore the temperature 
dependence in the replicated Hamiltonian Eq. ()6.20|) by 
some appropriate rescaling transformation: 



(7.1) 



(7.2) 



Eq. (|4.6p ). The relations between the two set of coeffi- 
cients appearing in Eq. (|7.ip and Eq. (|7.2p are shown in 
Table HVl 
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Eq. dLl 
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1 . The crossover exponent j/y for temperature is given 
by 



TABLE IV: Correspondence between the coefficients in 
Eq. dZIII and Eq. fT!^. 



The form of in Eq. (|7.2p leads to form of model 
parameters that are slightly more convenient for a study 
of consequences of the flow equations. Firstly, we ob- 
serve that since both /x^ and K are mapped into 
in Eq. (|7.2p . in order to preserve the Hamiltonian's 
functional form under the RG flow (a convenience for 
later analysis), the anisotropy scaling exponents lo in 
Eqs. (|6.35p must be chosen such that and K flow 
exactly the same. This condition gives 



Vk) 



(7.3) 



We shall see that this uj is precisely the physical 
anisotropy exponent which appears in all renormalized 
correlation functions. Since rj^ and tik are positive we 
observe that at long-scales the 2 : 1 anisotropy between 
xj_ and z of the harmonic theory is reduced below 2 for a 
fully nonlinear heterogeneous nematic elastomer. Using 
Eqs. (|6.35p and Table IIVI we can readily derive the flow 
equations for the new parameters of the rescaled Hamil- 
tonian in Eq. (fT^ : 
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^VK + VA 


A = 0, 


(7.4d) 



In above we have used the property of 77 exponents de- 
duced in the previous section, such that tjb 0, and the 
fact that at the critical point E the following relation 
holds: 

1 5 

{5-d--TjL- -m + Va) = 0, (7.4e) 

guaranteed by the underlying rotational invariance and 
the form of the nonlinear strain tensor w, Eq. (I4.32p . 
The flow of the external traction Sia , which we have not 
derived, can be inferred from other flows using equilib- 
rium fluctuation-dissipation relation, as we shall shortly 
demonstrate. 

From the flow equations (17. 4p we can draw several im- 
portant conclusions about the long wavelength physics of 
the elastic Hamiltonian Eq. (|7.2p . 



(7.5) 



Since both rjL and rjK are positive for c? < 5, in 
three dimensions, i/t is always negative. Therefore 
T flows to zero as / — *■ 00, and thermal fluctuations 
are irrelevant at long length-scales in physical 3d 
elastomers. 

2. xb flows to infinity, indicating that in the thermo- 
dynamic limit the bulk mode freezes out at long 
length-scales. For the same reason the cross cou- 
pling term xc is irrelevant at long-scales. 

3. v = fJ-/ fJ^L flows to a constant value as / 00. Thus 
a nematic elastomer exhibits universal Poisson ra- 
tios at long length-scales. 

With T flowing to zero and disorder variance A flowing 
to a constant, the long wavelength physics is controlled 
the zero-temperature, finite-disorder fixed point E. That 
is, long-scale properties of heterogeneous nematic elas- 
tomers are dominated by the ground state energy for a 
given random initial stress (j{x), a nontrivial problem due 
to the elastic nonlinearities of the model. In the scaling 
regime, v and A flow to constants, while xb and xc play 
no role since the bulk mode is frozen out. Consequently, 
all physical quantities become functions of only two pa- 
rameters temperature T and the external stress 5^- 



B. Correlation functions 

We now proceed to combine RG flow equations with 
a matching scale analysis to deduce long-scale properties 
of a nematic elastomer through the behavior of renor- 
malized correlation functions. To this end we look at the 
RG transformation of the renormalized (i.e., with all ef- 
fects of fluctuations and disorders included) correlation 
functions. For simplicity, we shall consider the zero ex- 
ternal stress case and take G"f (g) as an example. Using 
Eq. (|5.56p we have 

G^f (gl,g„T) = e('^+3-")'Gf (qe',g, e"',Te^-'), 

(7.6) 

where the exponent yr is given by Eq. (|7.5p and the 
anisotropy exponent lu fixed by Eq. (|7.3p . Using the 
results in Appendix [El we can separate the replicated 
correlation into a quenched part G^ and a thermal part 
G^: 



= {Ua){ub)+ (UaUb) - (Ua) (Ub) 5afj- (7.7) 

The thermal correlator is expected to be proportional to 
temperature T, with all other T dependence subdomi- 
nant due to the irrelevance of T at the fixed point E. 
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The quenched correlator is expected to be indepen- 
dent of temperature, and as we have discussed in Sec. lVIi 
is directly related to the nonaffinity correlation function. 

Combining Eq. (|7.6p with Eq. (|7.7p and proportionality 
of with T we obtain 



Gt{q^,q.) = e('*+3-)'Gf,(gle',g,e-'), 



(7.9) 



where we used Eq. (17. 3|) and Eq. (|7.4ap . From above, 
it is already clear that it is the dangerously irrelevant 
temperature, T that leads to different scaling behaviors 
of G^ and G'^. 

Now by taking e' — Q/q± in Eq. (|7.9p we find 



GzziQ^^qz) = qj_ 
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where 



(7.10) 



(7.11) 



and for simplicity of notation we have set Q — I (used 
units of the uv-cutoff). In the limit q^ — s- 0, G^^ must be 
a function of q± only. Therefore we have 



GL(gl,g. = 0)^gI^+''-r$-(0)«gI 



(7.12) 



Likewise, in the limit q± —> 0, Gj^ must be a function of 
only. This requires that as x — *■ oo. 



which leads to 



(7.13) 



(7.14) 



We can use the renormalized correlation function to 
calculate the real space fluctuations of the Uz phonon: 



{(Uz-Mf) 



dqzd"^ ^q±„T ^ 
-^^^GzAq±,<lz) 

dqzd'^-^qj_ ^-A+rjK^T ( Qz 



Similarly, we can deduce the renormalized quenched 
correlator Gf^- Setting = q^ in Eq. (|7.8p and using 
the Ward exponent identity, Eq. (|7.4ep . we find 

Gf.(<U,g.) = g-(«-2"-+"-)$f(^). 

qr-''~'^''^^-\ qz»ql. ^'-^'^ 



As argued above, the nonaffinity correlation function 
tziojtzi—Q) (defined in Eq. (|6.15p ) is expected to exhibit 
the same scaling behavior. 



Breakdown of linear response theory: nonlinear 
stress-strain relation 



The effective renormalized shear (and other elastic) 
moduli can be extracted from the phonon renormalized 
correlation functions. To this end, we introduce a uniax- 
ial traction SdzU^ in the elastic Hamiltonian Eq. (|6.20p . 
It can then be shown that, for an incompressible and uni- 
axial solid, the renormalized longitudinal shear modulus 
/iL is related to the renormalized correlation function G^^ 
via the following Kubo-like formula: 



To determine whether thermal fluctuation are divergent 
in a macroscopic limit, power-counting is sufficient. To 
this end we observe that the integral is dominated by the 
region where ~ , which gives 



{iuz~{uz)f)^qi, 



(7.16) 



S={d-3) + ^{VK-VL) = ^>0, ford = 3. 

Therefore thermal fluctuations are indeed finite in three 
dimensions for the renormalized heterogeneous elas- 
tomer, even though they are divergent with system size 
for a homogeneous one. 



= T 1 lini ql Gjz(g, 5). 



(7.18) 



Using Eq. ((7HP and Eq. (|7?^ . we find that the hnear 
shear longitudinal modulus (i.e., — > 0) vanishes in the 
long-scale limit according to 

fiL{S ^ 0) cx lim g^^/" = 0. (7.19) 

Analogous result can be derived for the linear transverse 
shear modulus. This implies that the elasticity of our 
system is strictly nonlinear for arbitrary small stress. 



C 



Us) 



00 



( g ) • (7.15) 



Critical Regime 

FIG. 14: External stress introduces new stress-dependent 
length-scales £,±,z(S), that diverge with a vanishing stress ac- 
cording to Eq. (|7.22|l . On longer scales the critical scaling of 
elastic moduli and the associated strictly nonlinear anoma- 
lous elasticity is cut off by ^±,z. The short-scale cutoffs 
for the critical regime are set by a combination of network 
mesh size (uv-cutoff) and the bare elastic constants. 



The origin of nonlinear elasticity is straightforward to 
understand within the RG formulation. The external 
traction Sia is a relevant coupling whose value grows un- 
der RG rescaling, thereby flowing away from the nematic 
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elastomer fixed point at Sia = 0. The scale at which a 
weak external stress rescales to a large, maximum (as de- 
fined by elastic model's uv cutoff and bare elastic mod- 
uli) effective stress defines a new length-scales £,±^z{S), 
beyond which the RG flow is cut off and the associated 
critical scaling breaks down. For a nonzero S, the scal- 
ing form of the thermal correlation function, Eq. (|7.9p is 
extended to be 

Gliq^, q.,S) = e^-^-^-^'GLCgl e\q. e-',Sey- 

(7.20) 

where ys is a positive exponent (calculable within our 
e RG expansion; see below) governing the flow of the 
dimensionless stress S (measured in units of a maximum 
stress for which the elastic model remains valid) under 
the RG transformation. Choosing S'e^"'' = 1, we obtain 



ys $ 



This automatically gives us length-scales 



(7.21) 



(7.22) 



relative to which all other scales (e.g., l/q±^z) are mea- 
sured. Setting q± to zero in Eq. (|7.2ip . we find that the 
correlation function behaves very differently in the large 
and small qz^z regimes. For qzCz{S) 3> 1, we must re- 
cover the critical scaling critical scaling form (|7.14p valid 
at S* = 



GliO,qz,S) 



for qzUS) » 1 (7.23) 



On the other hand, for qz^ziS) ^ 1, we must have 
GUO,qz,S) = q;25^(2.-4+,.)^ for 9,^,(5) «1 
= q^^S~'^'-/y% (7.24) 

whose form is dictated by the requirement of a qz inde- 
pendent shear modulus in Eq. (|7.18p . Using this inside 
Eq. (jTT^ . we find 



(7.25) 



that vanishes with a vanishing stress, 5—^0. 

To calculate the unknown exponent ys, we may study 
RG fiow equation for S near the fixed point E. However, 
we can also obtain the value of ys from the scaling form 
of the elastic free energy density, which is a function of 
T and S near the fixed point: 

/(T, S) = e"('^"i+")'~^^ ' /(Te^^', S-e^s'). (7.26) 

This scaling form is different from the corresponding form 
of thermal critical phenomena by a factor of e^^^K The 
difference is again due to the fact that long-scale proper- 
ties of the elastomer are controlled by a zero-temperature 
fixed point, E, with T a dangerously irrelevant coupling. 

As we have discussed above, the fact that the thermal 
exponent j/t is negative implies that thermal fiuctuations 



are not important at low temperature and thus the sin- 
gular part of the free energy, to the leading order, is inde- 
pendent of temperature, reducing simply to the ground 
state energy for given disorder realization and external 
traction. Thus neglecting this subdominant T depen- 



dence in Eq. (|7.26p and choosing / such that 



we find 



fiS) cx 



(d—l+uj+yT) 



(7.27) 



Using the fundamental definition of the inverse (nonlin- 
ear) shear modulus as the second derivative of the free 
energy with respect to the external stress S we obtain: 



l+Lj+yT-2ys) 



Ml(5') 

Comparing this with Eq. (|7.25p we find 



ys 



which then gives 



^l(S') cx S^-"k 



(7.28) 



(7.29) 



(7.30) 



This agrees with the result we obtained from a comple- 
mentary argument in Ref.[2^, based on the fact that the 
external stress S experiences only a "trivial" (not inde- 
pendent of other couplings) renormalization, analogous 
to an external magnetic field in an Ising model (0* the- 
ory). 

To the leading order in the e-expansion, we can read 
off the values of various exponents from Table Hill 
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which in three dimensions lead to 



(7.31a) 
(7.31b) 
(7.31c) 

(7.32) 



While this violation of linear strain-stress relation is quite 
weak in three dimensions (the exponent is small), it is 
qualitatively quite significant. 

As we have just seen from Eq. (|7.23p and Eq. (|7.24p . 
critical fluctuations of a nematic elastomer are sup- 
pressed by the external traction beyond the stress de- 
pendent crossover length-scales £,±{S) and £,z{S). Similar 
effect also arises from the application of a magnetic field 
along the nematic director, which adds a new term into 
the nematic elastomer Hamiltonian: 



6H„, =-^ja I d''x{h-hf ^^h" 



d'^xSnixf. (7.33) 



After integrating out the nematic director fluctuations 
using Eq. (|5.ip as we discussed earlier, SHm reduces to 



j d'^x iy^Uzf 



less relevant. (7.34) 
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This can be further transformed (using Wzz = dzUz + 
i (V_lM2)^) to reduce SH,n to 

6Hjn — —-jlah^ j d'^x (dzUz) + less relevant, (7.35) 

after shifting the reference state appropriately. Therefore 
7a ft.^ is equivalent to the stress component S along the 
director, which we considered in Eq. (|7.18p . With this 
mapping in hand, the following results straightforwardly 
follow: 

1. A magnetic field along the director introduces two 
crossover length-scales 

^^{h)^h-i, iz{h)^h-Tt, (7.36) 

beyond which the critical fluctuations are sup- 
pressed. 

2. A magnetic field along the director induces sponta- 
neous strain along the director, which scale as 

ezz h^^^^-^K (7.37) 

3. If we apply both a magnetic field h and an external 
traction S along the director, then the stress-strain 
relation is linear for S 'jah^ and is nonlinear for 

These predictions should in principle be experimentally 
testable. 



D. Stability of the anomalous nematic phase 

The validity of above analysis and the associated sta- 
bility of the critical phase requires that the nematic or- 
der is long-ranged. This translates into a constraint that 
elastic uniaxial anisotropy survives thermal fluctuations 
and network heterogeneity, a condition given by 

w = 2-i(r? + 7?K) > 1. (7.38) 

Using our one-loop exponents approximation in Table [TTTl 
we find that the condition (|7.38p is satisfied as long as 

d > d/c = 17/56. (7.39) 

This is clearly satisfied by the physical case of d = 3, in 
which 



fluctuation Sn is massively tied to the asymmetric lin- 
earized strain through Eq. (|2.39p . the fluctuations of Sn 
can be estimated by 

WW) « (r- - 1)-' {qliluzim (7.41) 
+ 2rqzq^{uz{q)u,{~q))+r^ Qz {\^±m^)) ■ 

Given the scaling RG analysis above, the dominant 
contribution on the right hand side is given by quenched 
fluctuations of the Uz phonon: 

WW) « ql{uz{q)){uzi-q)) = qlCtXq) 

where we have used the scaling form of Gf^, Eq. (|7.8p . 
In real space, we have 

WW) « j d'-'q^dqzql^'+'--^ (^^) . (7.43) 

A straightforward power-counting shows that the above 
integral scales as and therefore converges as long 

as w > 1, in agreement with the result in Eq. (|7.38|) . 

We have so far been ignoring the nonlinearities as- 
sociated with phonon, since they are less relevant 
at the Gaussian fixed point than the Uz nonlinearities. 
However, we have also seen in Eq. ()6.3ip that the har- 
monic quenched fluctuations of u± are actually divergent, 
where the unrenormalized correlator is used. To check 
whether these quenched fluctuations are still divergent 
in the renormalized theory, i.e., at fixed point E, we need 
to use the renormalized correlator Gf^^j, in Eq. (|6.3ip . 
Without detailed calculation, we observe that the renor- 
malized disorder variance A diverges and the shear mod- 
uli /i, fi vanish as g — > 0, according to our RG analysis. 
On the other hand, the bending constant K3 does not ac- 
quire anomalous dimension from Uz fluctuations. There- 
fore we deduce that the renormalized fluctuations of uj_ 
becomes stronger, i.e., more divergent than the naive har- 
monic analysis at the Gaussian fixed point. This suggests 
that we may not be able to ignore the columnar elastic 
nonlinearities associated with u±. However, a simultane- 
ous treatment of both smectic and columnar nonlinear- 
ities is a challenging open problem that we leave to the 
future research. 



E. Universal Poisson ratios 



ujd=3 « 1.574, (7.40) 

is significantly above unity, as required for nematic state 
stability. 

The same condition, Eq. (|7.38p can be obtained in a 
complementary way by looking at the real space fluctu- 
ations of the nematic director. Recalling that a director 



We conclude our analysis of long-scale anomalous elas- 
ticity, by observing from Table IIIII that the fixed point 
values of two coupling constants gL and g_L satisfy a re- 
lation: 

^ = 6. (7.44) 
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Combing this observation with the definitions of gL and 
g_L, Eqs. (|6.34p we find that the ratio between the long- 
scale effective transverse and longitudinal shear moduli 
approaches a universal constant: 



lim 

/ — >oo 



m(0 



6. 



(7.45) 



On the other hand, since 775 = at the fixed point E, we 
also have 



hm — — = 0. 

l^oc B{1) 



(7.46) 



That is, nematic elastomers are strictly (not just quan- 
titatively) incompressible at long length-scales. The ex- 
istence of these universal ratios characterizing a critical 
nematic elastomer phase is an analog of a well-known 
universal amplitude ratios at a continuous phase transi- 
tion. 

Let us now revisit the experiment shown in Fig. [9l 
where we impose a fixed strain deformation e^x in the 
X direction and let the sample relax without any macro- 
scopic reorientation of the nematic director. The other 
strain components are still given by Eqs. (j2.55p . if we re- 
place all bare elastic constants by the renormalized ones. 
Using Eq. (|7.45p . we find two universal Poisson ratios: 



-yy 



12 

r 



(7.47) 
(7.48) 



This prediction should be experimentally testable. 



VIII. SUMMARY AND CONCLUSIONS 

In this paper, we have developed a theoretical frame- 
work for the nonlinear elasticity of uniaxial nematic liq- 
uid crystalline elastomers, a fascinating class of materi- 
als which has an internal orientational order as a conse- 
quence of a spontaneous symmetry breaking. We have 
analyzed the soft Goldstone modes associated with the 
spontaneously broken rotational symmetry. In a strain- 
only description, nematic elastomers resemble conven- 
tional uniaxial elastic solids, but with a strictly vanishing 
shear modulus We have developed a model elas- 

tic free energy, in terms of an invariant strain tensor, 
which completely encodes the soft modes. We have also 
discussed its connection and difference with the popular 
neo-classical theory of nematic elastomers. 

We have discussed the angular momentum transfer be- 
tween the translational and orientational degrees of free- 
dom in nematic elastomers, characterized the antisym- 
metric part of the Cauchy stress tensor. We have also 
established the connection between stress tensor, couple- 
stress tensor, and the elastic free energy, and have calcu- 
lated these quantities explicitly using two model elastic 
free energies. 



We have also developed a complementary strain-only 
elastic model of nematic elastomers, and have derived 
a set of (rotational symmetry-enforced) Ward identities 
relating the coefficients of all anharmonic terms to those 
of quadratic ones. Using this minimal model, we an- 
alyzed the long-scale properties of ideal nematic elas- 
tomers and found that due to soft modes, the long-scale 
elastic properties are qualitatively modified by thermal 
fluctuations and polymer network heterogeneities, with 
the latter dominating over the former. 

Our key finding is that thermal and quenched fluc- 
tuations lead to the "anomalous elasticity" phenomena, 
familiar from a number of other soft-matter contexts, 
such as smectic and columnar liquid crystals, polymer- 
ized membranes and putative spontaneous vortex lat- 
tices in ferromagnetic superconductors. Specifically, as 
a result at long scales nematic elastomers are charac- 
terized by singular, length-scale dependent shear elas- 
tic moduli, a divergent splay elastic constant, long-scale 
incompressibility, universal Poisson ratios and, a non- 
Hookean (nonlinear) stress-strain relation down to ar- 
bitrary small strains. Furthermore, we show that these 
properties are universal, and are controlled by a non- 
trivial zero-temperature flxed point, constituting a qual- 
itative breakdown of classical elasticity theory. Thus, 
nematic elastomers constitute a stable "critical phase", 
characterized by universal power-law correlations akin to 
a critical point of a continuous phase transition, but ex- 
tending over an entire phase, as illustrated in Fig.[T31 We 
have also found that for weak disorder and low thermal 
fluctuations the anomalous elasticity stabilizes long-range 
nematic elastomer order in three dimensions. Such ori- 
entational order in a 3d disordered system with purely 
orientational degrees of freedom, e.g., spin system or a 
nematic liquid crystal, is known to be impossible. Thus a 
stable orientational order in nematic elastomers is a con- 
sequence of stabilizing interplay between orientational 
and elastic degrees of freedom. 

However, experimentally, at low temperature, ne- 
matic elastomers crosslinked under isotropic conditions 
always exhibit a polydomain nematic director pattern 
with short-range nematic order, in apparent contradic- 
tion with our e-expansion results. There are a number 
of possible explanation for this disagreement. Firstly, 
and most likely in our view, real nematic elastomers are 
quite likely characterized by strong heterogeneity, which 
may not be captured by our perturbative RG analysis. 
It is also possible that e = 2 is too large to trust the e- 
expansion. Although qualitatively subdominant, colum- 
nar nonlinearities may play an important quantitative 
role. Also our analysis neglected the transverse part of 
the random stress. This naturally couples to the ne- 
matic topological defects, and therefore for strong dis- 
order may lead to a proliferation of the latter, thereby 
destroying long-range orientational order. This problem 
is being actively studied by the authors. Finally, it is also 
quite possible that nematic order is indeed stable in real 
weakly-heterogeneous elastomers, but (as in the case of 
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the random- field spin systems [85|) requires long equili- 
bration times or "field-cooling" to realize. If the nematic 
orientational order is indeed unstable to disorder, but 
with a long orientational correlation length, our analysis 
and predictions are expected to be valid in the resulting 
polydomain nematic for a range of length-scales between 
the nonlinear elasticity scale ^nl and the orientational 
correlation length-scale ^nem , with latter diverging in the 
weak heterogeneity limit. Further theoretical and exper- 
imental studies are clearly needed in order to clarify the 
true nature of the nematic elastomer ground state. 

Our present work leaves open a number of interest- 
ing questions. One very important one is the fate of 
the isotropic-nematic transition in a heterogeneous elas- 
tomers. Also, we exclusively focussed on the ideal case, 
ignoring the fact that all mono-domain elastomers are 
crosslinked under pre-stretched conditions, that imprints 
a weak uniaxial anisotropy, that was recently shown to 
lead semi-soft elasticity. Thus for a direct compari- 
son with experiments our theory must be extended to 
a detailed treatment of semi-soft elasticity. Interesting 
work on this subject was recently carried out by Ye and 
collaborators, [l^ It is quite clear from their and our anal- 
yses that imprinted anisotropy cuts off the critical phase 
(anomalous elasticity) phenomenology at long scales and 
weak external traction. We leave further detailed work 
on these subjects to future investigations. 
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APPENDIX A: ELASTIC MODULI IN THE 
EFFECTIVE STRAIN-ONLY MODEL 

Here, for completeness we list all eleven coefficients 
appearing in the effective elastic Hamiltonian Eq. (|4.2ip : 



ai_ - 2Ci ($i+2Ci$2 + 3Cl3'3), 
a, = 2Cf ($i + 2Cf $2 + 3Cf$3) , 

Mz± = 4CiC|| (2 $2 + 3 (Ci+Cfi) $3), 
3 

= 4 ^(z.j)<i>«,C,f^'V8C|<i>2 + 24Cf<I>3, 

= 4^(z.j)<I'.,CifC?', 

A = 4^(z.j)$,,Cf+^'\ 
M = 4Ci ($2 + 3Ci3'3), 

61 = 8ciC|| (3*3 + 2*12 + 3 (ci + cf) *i3 

+ 4C|f$22+6C|| (Ci+2C') $23 

+ 9cf (ci + cf) $33) , 

b2 = 8ClCf (2*12+3 (Ci+Cf) $13+4 Ci $22 

+ 6 ci (2 ci + C|) $23 + 9 ci (ci + c') $33) , 

63 = 24ClC'$3, 

C - 8ClC| (4*22 + 12 (Ci+Cf) $23 

+ 9 (ci + C|)' $33) . 

APPENDIX B: WARD IDENTITIES REVISITED 

In this appendix we re-derive the Ward identities 
Eqs. (|4.22p via an alternative, more straightforward ap- 
proach, whose advantage is that for d > 3 it is indepen- 
dent of space dimensions. 

To this end, let us consider a c?-dimensional ideal uni- 
axial nematic elastomer. Using Eq. (|2.13p . it is easy to 
see that the most general expression for soft deformations 
with constraint n = no is given by 

no ^ no, A = Aa„ O Ar^\ (Bl) 

where O is an arbitrary rotation and Ang is given by 
Eq. (|2.12p with n replaced by no: 

Afi„ =Cil+(C'-Ci)^iorio, (B2) 

Let us choose the coordinate system such that no is 
parallel to the z axis and consider a rotation O in the xz 
plane. All other d— 2 axes are unchanged by the rotation 
and therefore do not need to be considered. Therefore 
we only need to keep track of (xx), (xz), {zx), and {zz) 
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components of the tensors O and An^ , and can express 
them as 2 X 2 matrices: 



Given that is a soft deformation, Eq. (|B5|) . all terms 






Cn 



(B3) 



O 



acting in this 2d subspace. Using Eq. (jB2p and Eq. (|B3|) . 
we find the Lagrange strain tensor corresponding to this 
soft deformation is given by 




sin^6i 



-2^sin20^ 



(B4) 



where we have expanded up to order of 6^. 

For this particular soft deformation (|B4p . the first 
three terms of the elastic free energy, Eq. (|4.2ip are given 
by 



1(1 
2 V 



4 r 



(B5a) 
(B5b) 
(B5c) 



which are all at order of 9^ . It is also clear that all other 
terms explicitly shown in Eq. (I4.2ip are at least of order 
6"^, and with others (not shown) are even of higher, at 
least 6^ order. Therefore, in order for the elastic energy 
Eq. (|4.2ip to vanish at order of 6^ , we must have 







(r-1) 



a± e 

2/)2 



4r 



2 



(B6) 



2r 



1 



1 



which can be seen to be identical to the first Ward iden- 
tity, Eq. (|4.22ap . after using the relation Eq. (|4?23)) be- 
tween r and a. 

We can then make an appropriate choice of the ne- 
matic reference state such that the first three terms in 
Eq. (|4.2ip vanish simultaneously. The remaining eight 
terms, all of order of 9"^, can be reformulated into 

1 2 1 

n = -B, {e,^ + 2a, e^-) - A (e„; + 2q_l e^J 

+ A^_L {ezz + 2a2 e^j) (eu + 2a± e^j) 
+ 11 {e^j + 2a'^eziezjf + 



-2 2 

^ ^zi ^zj 



(B7) 

The new coefficients a^, aj^, a'^ and c can be deter- 
mined by comparing coefficients of every terms between 
Eq. HZl) and Eq. g^B): 





= 2azBz + 


2q:_l A3j_, 


(B8a) 




= 2az Xzj_ ' 


F 2a^ A, 


(B8b) 


h 


= Aa'j^ji, 




(B8c) 


c 


= 2alBz + 


4 az aj_ \zi__ + 2 a]_ A 






+ 4(al)V 


+ c 


(B8d) 



in Eq. (jB7p order by order in 9 and in particular to or- 
der 0* must vanish identically. Furthermore, the elastic 
energy must be nonnegative for arbitrary strain defor- 
mation. The only way to satisfy both conditions is by 
requiring that every single term in Eq. (jB7p vanishes for 
the soft deformation Eq. (IB4p . This imposes stringent 
constraints on coefficients a^, a_L, a', and c as follow: 



1 



Ctz 

c 



r-1 



= a, = -- 



1 



0. 



(B9) 

(BIO) 
(Bll) 



Under these conditions, Eqs. (|B8p reduce to the last four 
Ward identities in Eqs. (|4.22p with iiz± set to zero, while 
the elastic energy Eq. (|B7p reduces to Eq. (|4.25p . 



APPENDIX C: DERIVATION OF RG FLOW 
EQUATIONS FOR AN IDEAL HOMOGENEOUS 
ELASTOMER 

In this appendix we present a detailed derivation of 
renormalization group flow equations (|5.32p near three 
dimensions for elastic constants characterizing a nematic 
elastomer. We will also show explicitly that the form of 
Hamiltonian Eq. (j4.3ip is preserved by the renormaliza- 
tion group transformation, as we have already argued in 
Sec|V] based on symmetry grounds. 



1. Momentum shell RG 

The principle of renormalization group transformation 
(RG) is to transform a Hamiltonian H for a fluctuating 
field u{x), and characterized by a set of coupling con- 
stants c and a uv wavevector cutoff Q to a Hamiltonian 
of the same form for a coarse-grained field u'{x'), the 
same momentum cutoff Q, and characterized by new set 
of couplings constants c'. The goal is to establish a rela- 
tion (RG flow) between coupling constants c' and c under 
such successive RG transformations. Armed with the RG 
flow of the coupling set c, as well as the relation between 
the original field u{x) and the coarse-grained one u'{x') 
allows one to extract the scaling behavior of correlation 
functions of u{x) . 

A momentum-shell RG (MSRG) consists of two steps: 
(i) an integration of short wavelength fluctuations of the 
field u{x) (i.e., coarse- graining), whose Fourier transform 
has support in the momentum shell Q e~^^ < q < Q, giv- 
ing a coarse-grained field theory with uv wavevector cut- 
off Q e~^\ and a set of modified coupling constants; (ii) a 
rescaling of coordinates x and field u[x) so as to restore 
the uv cutoff to the original value of Q. The rescaling 
transformation (ii) leads to further modification of the 
coupling constants c. Although this step is in principle 
unnecessary, it facilitates the establishing of the relation 
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between coupling constants c and c', characterizing the 
two Hamiltonians. 

For a uniaxial system such as a nematic elastomer, 
there are generally two different momentum cutoff 
and (5_L, where subscripts || and _L denote directions par- 
allel and perpendicular to the nematic director, respec- 
tively. However, the principle of universality guarantees 
that the long length-scale physics of the system is inde- 
pendent of different choices of momentum cutoff. It is 
convenient to choose a cylindrical cutoff scheme, where 

(Qy = 00,Q_L = Q). 

For most applications it is sufficient to focus on the 
partition function Z[c, Q], given by 

rQ 

Z[c,Q] = / Pwe^^I"'^!, (CI) 



where 



where the functional integral is over field u{x), whose 
Fourier transform has support in the region \q±\ < Q, 
and we have explicitly shown the dependence of the par- 
tition function on the uv cutoff Q and the set of coupling 
constants c. 

To carry out the coarse-graining procedure, as illus- 
trated in Fig. IIOI we decompose the field u{q) into a sum 
of high- and low-wavevector parts u^{q) and u^{q), 

u{q)^u>{q)+U<{q), (C2) 

where u^{cj) has support in the outer cylinder shell of 
the q space: 



U>{q) 



^0, if Qe-'<\q^\<Q, 
= 0, if 0<|ql|<Qe-', 



(C3) 



while u^{q) has support in the inner cylinder of the q 
space: 

^<(-^N0, if Qe-'< 1,-11 <Q, 



Because of the translational invariance, the harmonic 
part of the Hamiltonian can be separated into higher and 
lower momentum pieces, giving 

H[u] = Ho[u< + u>] + Hi[u] 

= Ho[u<] + Ho[u>] + Hj[u\, (C5) 

where Hi[u\ contains all the anharmonic terms in u{r). 

This allows one to separate the functional integral into 
a product of two parts, 



, (C6) 



where we have also introduced a short hand for functional 
integral over and . Carrying out the functional in- 
tegral over the high-wavevector field u'' , we obtain an- 
other field theory with a wavevector cutoff Q e^"^' and a 
modified Hamiltonian defined by: 









Du< / Du> 







Z[c, Q] 



.-Ho 



-Ho[u>]-Hi[u<+u>] 



/O 1 



-Ho[S>,{X}] 



(^>0> 



Z^ 



(C8) 
(C9) 



This allows us to define a coarse-grained Hamiltonian 
as a functional of and new coupling constants c + SgC 



H[u<,c + 6gc] = Ho[u>] - log(e-«^[«"+""'^l)>, (CIO) 



in terms of which the partition function (|C7p can now be 
written as 

Z[c,Q] = Zo>y"^e-^["^'^+*«=l = Z>Z[c + SgC,Qe'^'] 

(Cll) 

This thereby transforms a computation of Z into a func- 
tional integral over a coarse-grained field (a?) , momen- 
tum cutoff Qe^^'' and coupling constants c + SgC. 

For convenience, we furthermore rescale the spatial co- 
ordinates [13] such that the wavevector cutoff is restored 
to Q: 



x = x'e^\ q^q'e~^K 



(CI2) 



This leads to further transformation of the coupling con- 
stants, which we denote as 5rC [sst : 

Z[c,Q] = Z>Z[c + 6gC + 5rC.,Q]. (C13) 

This allows us to summarize the effect of the RG trans- 
formation in terms of a transformation between coupling 
constants, given by 



c + 5c = c + SnC + Sj-c, 



(C14) 



with SC — 6gC + SrC proportional to an infinitesimal SI. 
The RG transformation can therefore be summarized by 
ordinary differential equations for c{l), describing the RG 
flow of coupling constants under successive infinitesimal 
renormalization group transformations. 



2. Derivation of RG flow equations, Eqs. (|5.32|) 

The key non-trivial step in the RG transformation is 
the calculation of the second term in the right hand 
side of Eq. (jC10|) . Near the upper-critical dimension, 
due (above which nonlinearities and fluctuations are no 
longer important), this step can be carried out perturba- 
tively in nonlinearities (which amounts to be perturba- 
tively in e = due — d) using cumulant expansion: 



(C7) 
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log(e-«^P'=+«^^{Ml)> = {Hj)> - \ {{Hj)> - {{Hj)> f) + ■■■. 



(C15) 



The total nematic elastomer Hamiltonian is given by Eq. (|4.31ap . with the effective strain tensor w defined in 
Eqs. (|4.32p . The harmonic part Ho[u] is given by: 



The nonlinear part Hi[u\ contains all cubic and quartic terms given by 

Hl[u\ — -ffcubic + -ffquartic, 

HqunMc = Bijki{diUz){djUz){dkUz){diUz), 



where 



1 



Bijki = + + Bz-'i A^j_) {5ij5ki + SikSji + SiiSjk) ■ 

These can be conveniently represented diagrammatically, as illustrated in Fig. 1151 



(C16) 



(C17a) 
(C17b) 
(C17c) 



(C18a) 
(C18b) 



Hi[u] = 





FIG. 15: Feynman diagrams for cubic and quartic nonlinearity, respectively. The solid lines represent field diUz, while the 
wiggly line represents field 



We calculate the cumulant expansion (|C15|) up to 1-loop order, with every term in this expansion representable by 
a Feynman diagram. We refer the reader to reference [s^ for a introduction of graphical representation of a cumulant 
expansion. Since we are after a correction that is a functional of coarse-grained fields, with only high-wavevector fields 
integrated out, all Feynman diagrams have harmonic propagators of as internal lines and u"^ as external lines. 

A basic property of a cumulant expansion, is that no disconnected diagrams appear in the expansion. Furthermore, 
momentum conservation requires that all so-called one-particle- redwciWe diagram, which can be separated into two 
pieces by cutting one internal line, vanish identically. An example is shown in Fig. 1161 Therefore we only have to 




FIG. 16: A one-particle-rerftJcifcZe diagram. Because the internal line carries a high momentum p > Q e while the external 
Une carries lower momentum q < Q e~^' , the diagram vanishes identically under MSRG by momentum conservation. 



keep track all one-particle- irrerfMci6/e diagrams in our calculation. 

The first-order cumulant in Eq. (jC15[) can be represented as the sum of four Feynman diagrams shown in Fig. 1171 
It is easy to see that the first two diagrams are just Hj[u'^]. The last two diagrams generate terms such as 



j d'^x A,j[u<] and j d'^x (V_lz?<)^ 
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FIG. 17: Non-vanishing Feynman diagrams for the first cumulant in Eg. lClSl Other diagrams either vanish or are only constants 
independent of fields u^. External lines represent field uf (solid lines) or Aij[it^] (wiggly lines). Internal lines represent or 



which do not appear in the original Hamiltonian Eq. (j4.3ip . However, they can be expressed as a linear combination 
of Wzz and wa and therefore can be eliminated by an appropriate shift of the nematic reference state, as we will show 
explicitly later in this appendix. 

By calculating the second order cumulants we find corrections to all quadratic terms in the Hamiltonian, Eq. (|4.31ap . 
which can be represented by the diagram in Fig. [TlJ Because the functional form of the elastic Hamiltonian Eq. (|4.3ip 
is preserved by the renormalization group transformation, these are only terms that we need to calculate. Rotational 
symmetry (Ward identities) ensure that from these corrections to quadratic terms we can infer corrections to all 
anharmonic (cubic and quartic) terms appearing in Eq. 14.311 

The part of the second order cumulant that renormalizcs elastic constants Bz, \z^_, A and /i is given by 



Aij[u^]Au\u 



' ^ PiPjPkPl Gzzipf, 



(2 7r)« 



where 



df^P Q 



d-l 



SI / dfld- 



dpz 



(C19) 



(C20) 



is the integral over the momentum shell (— oo < pz < oo, Qe^*' < p± < Q). dVld-i is the differential surface element 
of a unit sphere in a d — 1 dimensional space. The integral of the product pi pj pk pi over the d ~ 2 sphere can be 
easily calculated to give: 



d^d-i Pi Pj Pk Pi = ^i-^ {SijSki + SikSji + SiiSjk) ■ 



(C21) 



Therefore Eq. (|C19|) reduces to 



dpz 



-oo ifip^z+KQ^r 

^-—(^{Bz-Xz±)Hdzuf)^+2{Bz^Xz±){X + pi-Xz±){dzuf){W^-u<) 



d-3 



SI 



2{2TT)d-^ (d2_l)/f^ 



+ {x + ^l-Xz±f{s/^■u<f + -^i'{^,u<f 



(C22) 



According to the momentum shell RG scheme outlined above, the preceding set of diagrammatic correction is to be 
equated to: 



-{SgBzidzU<f + 2SgXz^{dzU<){d,u<) + {SgX + <5<, ) ( 9. « f ) ^ + Sg^i{^,u<f}, 



(C23) 



which gives graphic corrections to elastic constants Bz, X^, B± and fi as shown in the second terms, on the right 
hand sides of Eqs. (|5.32all5.32dp . 



49 



The part of the second order cumulant that renormahzes the splay constant K is represented by the wavevector 
dependent part of the Feynman diagram shown in Fig llll The corresponding correction to Hamihonian Eq. (j5.2p is 
given by 



1 1 \ d 

-5gKq\_\u,Xq)\^ = -2'^ ■ 2 g»gjgfcg;|nz(g)p- ^^^^^^ J 

[{A,.m[u{p + q)]Ajn[u{-p- q)])o {PrnPn\Uz{p)\'^)o 
{Aijn[u{p + q}]{pn + qn)Uz{p + q))o {PrnUz{p}Ajn{p)f)o] 

na^iQ'^-^Sl {Bz (A + 2 ^) + 12 /i(A + ^) - 4 ^ C - C 



(2 7r)'* 



-ql\uM\\ (C24) 



2(2 tt)''-! 16(B_l + ^i)VK? 

and can be readily computed using Mathematica. Since we are after a result that is lowest order in e, the graphical 
correction can be approximated by its value in c? = 3. From above we obtain the graphic correction to K used in the 
second term on the right hand side of Eq. (|5.32ep of the main text. 

Finally, we apply the rescaling transformation R^^{e'^ ^'■)R±{e^'') to restore the momentum cutoff to Q. From 
Eqs. (|5.23p and Eqs. (|5.24p it is easy to see that the corrections to various elastic constants from this rescaling are 

{6rBz,6r\z±,SrX,SrfJ-) = (d + 3 - 3 w) , Az_L , A, 

(C25a) 

SrK = {d-l~Lu)KSl. (C25b) 

Assembling these graphical and rescaling corrections to all elastic constants, we obtain the RG flow equations (|5.32p 
of the main text. 



3. Generation of linear strain terms 



The vanishing of linear (in w) terms in the "bare" elastic Hamiltonian Eq. (|4.31ap is a consequence of our choice 
of the nematic reference state. However, for this choice of the bare Hamiltonian, due to thermal fluctuations and 
elastic nonlinearities, thermal averages (dzUz) and {diUi) do not vanish in equilibrium. This implies that the nematic 
reference state, around which we expand the elastic Hamiltonian, is not the true equilibrium state in the presence 
of thermal fluctuations. This manifests itself through a generation of a nonzero linear terms azWzz + ci±Wii by the 
coarse-graining. 

However, Ward identities, Eqs. (|4.22ap . dictate that we must also generate a term /iz_L(Vj^Wz)^, with /i^j^ satisfying 
Eq. I4.22al This guarantees that such terms assemble into a nonlinear strain w term, and can therefore be shifted 
away by re-expanding the coarse-grained Hamiltonian around the true (thermal fluctuations corrected) ground state. 
Thus around the new state, both linear coefficients az,±, as well as the shear modulus fiz± are guaranteed to vanish. 
In this section, we explicitly verify that this is indeed the case. 

In the process of renormalization group transformation, cubic nonlinearities Aij[u\diUzdjUz generate terms linear 
in the phonon field u, the so-called tadpole diagrams. At one-loop order, this is represented by the third diagram in 
Fig. [17] and is given by 

-±-P—^-—= (-(A + ^,- A,^)V^ • + {Bz - Xz±)dzuz) . (C26) 



At the same time, the fourth diagram in Fig. [T7| as It is straightforward to see that the sum of these two 
well as a part from the diagram in Fig. [TT] also generate corrections can be written in the form of 



a quadratic contribution 



SazWzz + 5a±_wu (C28) 



with 



_ ^d-iQ^-^Sl {Bz - Xz±) 

,Q^'Hl 1 " 2(2 TT)'^-! 8/?^ ' ^^ ^^ 

(C27) 2(2^)'^-! ■ 
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These linear terms in Eq. (IC28P can then be eliminated 
by an appropriate shift of the nematic reference state. 

4. Renormalization of cubic and quartic elastic 
nonlinearities 

As we discussed in the main text, rotational invari- 
ance requires that perturbative corrections to cubic and 



quartic nonlinearities be related to those of the quadratic 
terms found above, so as to preserve the form of the non- 
linear elastic Hamiltonian. 



Diagrammatically, corrections to cubic terms 
Aij[u\diUzdjUz can be represented by Feynman di- 
agrams in Fig. [THl Detailed, somewhat technically 
involved calculations, that we omit here give: 



iQ SI 1 / 
2(2 tt)" 64 y'K-^ fl ^ 
+ (2 (A - \zi_) [Bz - 2 + A) + 2 (S, - 3 + 2 A) /i -f /i')(V± • 
+ fi^idiUj + djUi)'j {diUz){djUz). (C31) 
Interpreting this as a correction to elastic constants of a coarse-grained cubic terms from Eqs. (jC17|) and Eq. (|C18[) 
^[{SgBz - SgXzj_){dzUz)Sij + {6gXzj_ - (5gA)(V_L ■ uj^)Sij - 6gn{diUj + djUi)]{diUzdjUz), (C32) 

we find that these corrections, 6gB, 5gXz±,- ■ ■ are identical to those obtained from the quadratic term above. 
Similarly, fluctuation correction to the quartic term, represented by diagrams in Fig. 1191 is given by 

n,^,Q^-^6l (2 (B. - 2 A.^ + A)^ + 4 (B. - 2 A.^ + A) ^ + 3 ^,^) 
6gH,^^rUe - 2(2^)._i 256 ' ' ^^^^^ 

which when identified with 

i {dgX + 2 6gtl + SgBz - 2 5gXzl_) {V z)\ (C34) 

again gives the same corrections to elastic moduli as quadratic and cubic terms. 

Thus, as required by symmetry we have demonstrated explicitly that the functional form of the elastic Hamiltonian 
Eq. ()4.3H) (encoded through the nonlinear form of the strain tensor w) is preserved by thermal fiuctuations and 
therefore by the coarse-graining RG transformation. 

5. Solution of RG flow equations for d 7^ 3 For d > 3, e < and it is clear that only the Gaussian 



In this appendix, we solve the renormalization group 
equations (|5.37p for d ^ 3. 
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FIG. 19: Feynman diagrams renormalizing quartic term 



fixed point g£ = = is stable. This means that elastic 
nonlinearities are irrelevant, and the harmonic approxi- 
mation becomes asymptotically exact, with subdominant 
corrections (from elastic nonlinearities), that can be com- 
puted in a controlled perturbative expansion in elastic 
nonlinearities. 

For d < 3, it is clear that elastic nonlinearities destabi- 
lize the Gaussian fixed point, with two nonlinear, dimen- 
sionless couplings ghil) and g±{l) expected to flow to a 
finite fixed point. We first look at the flow equations for 
the dimensionless ratios x{l) and y{l), Eqs. (|5.38p . As 
discussed in Sec. El for physical elastomers, the bare val- 
ues of X and y are much less than unity. Furthermore, for 
a finite and positive gh{l) it is easy to see from Eqs. (|5.38p 
that both x{l) and y{l) flow to zero as I oo. Therefore 
to a very good approximation, which becomes asymp- 
totically exact, we can set x{l),y{l) to zero in the flow 
equations for gi, and g_\_, Eq. (|5.37ap and Eq. (j5.37bp . 
which then drastically simplify to: 



dl 

dg± 



gL (5gg + 34g^.gL + 26gl) 
l6(gL + g±) 

gi_ {gl + 32g^.gL + 28ffi) 
16(.gL + 5-L) 



(C35a) 



(C35b) 



Solving these two equations we find 4 fixed points, that 
we list along with the corresponding set of rj exponents 
in Table El 

The flow pattern of and g± under renormalization 
group transformation for d < 3 is shown in Fig. 1201 It is 
interesting to note that the critical exponents r/L and rjx 
at fixed point S are identical to those of smectic liquid 
crystal in 3 — e dimension, first studied by Grinstein and 
Pelcovits [s^l . Since this fixed point S is attractive for the 
bare value of the coupling g± = (proportional to the 
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gL 


gl 
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4e 
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^9 
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TABLE V: Fixed point and corresponding dimensionless cou- 
plings and Tj exponents, with x* = = at all fixed point. 



transverse shear modulus /i) a nematic elastomer with 
/i = shares the same long wavelength "anomalous elas- 
ticity" with a smectic liquid crystal. This is not merely a 
coincidence, as it is clear from Eq. (|4.3ip that for a van- 
ishing fx, the phonon field u± can be integrated out. The 
resulting effective model is exactly a smectic liquid crys- 
tal with a one dimensional phonon field Uz and a shifted 
compressional modulus. Another way to understand this 
result is to note that with a vanishing in-plane transverse 
shear modulus fi, our elastomer model Eq. (|4.3ip is essen- 
tially a stack of liquid membranes with Uz the inter-plane 
displacement. Physically, this is clearly equivalent to a 
smectic liquid crystal. 

Fig. 1201 shows that this smectic fixed point, S is un- 
stable to a turning on a finite in-plane shear modulus, 
with the fiowing toward a globally stable fixed point E 
at finite couplings g± and grmL ■ It thus controls the long 
length-scale physics of ideal uniaxial nematic elastomers 
below three dimensions. 

Although as noted above, physical elastomers are usu- 
ally characterized by small x and y parameters, it is inter- 
esting to examine a special case of 2/(0) = ±1. As is clear 
from Eqs. (|5.38p for y = ±1 neither x{l) nor y{l) fiow, 
giving two unstable fixed lines (0<a;<oo,y = ±l) that 
correspond to very special systems. Setting y — ±1 in 
Eqs. (j5.37p . we obtain considerably simplified flow equa- 
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9, 

X 



tions for gL and g± for this special case: 



S 



FIG. 20: Flow diagram for coupling constants ql and g± of an 
ideal homogeneous nematic elastomer. 5* is the smectic fixed 
point discovered by Grinstein and Pelcovits '30] . The globally 
attractive fixed point E controls the long-scale properties of 
an ideal homogeneous nematic elastomer, characterizing its 
universal anomalous elasticity. 



dl 



dg± 

dl 



gl Agi^ (3 T V^) + g± (89a; T 210^^ + 234) 



16 (ffL (3 T V^f + Q9±x 



gl {18g±x + 5L (55x T 222Vi + 252)) 
16 (^L (3 T V^f + 95±^ 



(C36a) 



(C36b) 



I 

Solving these two equations, we find following fixed lines parameterized by x: 

1. Line of unstable fixed points: 

{gl = 5l = 0,?/* = ±1) 
These two lines are unstable with respect to all directions. 

2. Line of mixed fixed points: 

(.9L = 4e,.gl = 0,y*=±l) 
These two lines are unstable in the g± direction and stable in the gi^ direction. 

3. Line of mixed fixed points: 

(5L=0,5l = 8e,2/*=±l) 
These two lines are unstable in the direction and stable in the g± direction. 

4. Line of stable fixed points: 

, 8(19xt6v/^+9) , , 

9l2;=F 222Vi+ 252^-^ ^ 

These two lines are stable with respect to both and g±. rj exponents for various elastic constants explicitly 
depend on a continuous parameter x: 

2(19xtQVx + 9) 



91xt 222Vx + 252' 

(C37a) 

= 16(8xT24Vi^+27)^ 
' 91xT 222Vi+ 252 ^ ' 

Therefore systems characterized by — C^/B = 1 exhibit a continuous family of universal long length-scale 
elastic properties. It is not clear to us at the moment what physical system is characterized by y = 1 and therefore 
shares properties given by these fixed points. 



APPENDIX D: MODEL OF HETEROGENEITY based disorder model in Eq. 
IN NEMATIC ELASTOMER 

H^ = - I d'^X [TrFi(X)A'^A + TrF2(X)A'^Q(X)A] , 

In this appendix we derive the disorder Hamiltonian 
Eq. (|6.5p starting with a phenomenological, symmetry- 

(Dl) 



53 



where A is the deformation gradient relative to the 
isotropic reference state and is defined in Eq. (|2.ip . Fi 
and F2 are random tensor fields in the reference space. 

Using Eq. (|2.9p to express A in terms of the deforma- 
tion gradient relative to the nematic reference state, A, 
and defining two new tensor fields Fi{X) and F2{X) by 



Fi — Jq AoF^Aq . 



(D2) 



Choosing the direction of the nematic order in NRS 
to be the z axis, uq = z, to first-order in Sh (where 
6n-z = 0), the nematic order parameter Q for the current 
state is given by 



where 



Q = Qo + S (z Sn + 6h z), 



(D5) 



where Jq is a Jacobian factor 



Jo = det 



dx 



(D3) 



we can express the disorder Hamiltonian Eq. (|D1[) as 



Ha 



TrFiA'A + TrFaA'QA 



(D4) 



Qo = {zz - -I), 



(D6) 



Using Eq. (j2.34bp and Eq. (|D5|) . we find that up to 
linear order in u and 5h^ 



TrFiA^A = const. + 2(Fi)ah£a6 



(D7) 



and 



TrFsA^QA = const. + 2S{F2)^, [ ^ {d,u,) - ^d,u, + 5h^ + ^^"^^ ^^ ^(A)..^.^. - \s{F2\jd,u,, (D8) 

I 



where, as usual, indices i and j are limited to the d — 1 
dimensional subspace perpendicular to z axis. 

We proceed to integrate out the fiuctuations of the 
nematic director 5h. To the lowest order this amounts 
to making the replacement Eq. (|5.ip in Eq. (jD8|) . It is 
interesting to see that this replacement exactly cancels 
the antisymmetric strain components azi in Eq. (jPSp . so 
that as the final result, the disorder Hamiltonian only 
depends on the linearized symmetric strain 

- J d'^xY^ aabix)eab, (D9) 

a.b 

where the random stress field (t{x) is given by 

a.. = 2{F,)zz + ^^^-^S{F2)zz, (DlOa) 
d 

a,, = 2(Fi)., -^^(Fa).,, (DlOb) 

/ d—2 r + 1\ 
= a,, = 2(Fi),,+ f^- + — -^j5(F2),,. 

(DlOc) 

We note again that because random stress couples lin- 
early to strain, the ideal nematic reference state defined 
by r(x) = X is not the real ground state in the presence 
of this quenched random stress. 

At three dimension, a symmetric tensor field cr, gener- 
ically contains a longitudinal part satisfying V x er^ = 0, 
and a transverse part satisfying V- cr^ = 0. In Eq. (|D9p . 
however, because <t is coupled to the linearized strain, its 
transverse part cr^ has no bulk effect to the system, since 



it disappears after a simple integration by parts. There- 
fore only the longitudinal part of the random stress are 
included in our model. On the other hand, if we were 
to keep higher order terms in linearized strains in our 
derivation, we would have found the reduced disorder 
Hamiltonian to be 

Hd = - J d'^xaab{x)eabix), (Dll) 

rather than Eq. (|D9p . That is, the random initial stress 
field cTy is coupled to the nonlinear Lagrange strain, in- 
stead of the linearized strain . This result is of course 
guaranteed by rotational symmetry: the linearized strain 
is not rotational invariant in the embedding space, while 
the nonlinear Lagrange strain is. In the resulting model, 
then the transverse part of the random initial stress, afj 
will couple to the nonlinear part of the Lagrange strain as 
shown in Eq. (j6.2p . In this work, however, we shall only 
analyze the simplified disorder Hamiltonian Eq. (jD9p . i.e. 
we shall completely ignore the transverse part of the ran- 
dom stress. 

Because of the vanishing of quadratic term in the 
elastic energy Eq. (|5.2p . it is not difficult to see (using 
power-counting in previous sections) that the compo- 
nents azi are the most relevant part of the random stress 
tensor field aab- We therefore focus on these components, 
neglecting all others: 

Hd^-Jd''xJ2 > (D12) 

i 

where again i is restricted to the d — 1 dimensional sub- 
space perpendicular to z. The random field azi{x) is 
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assumed to be Gaussian with a short range correlation: 

APPENDIX E: REPLICA TRICK 

For completeness, in this appendix we briefly review 
a convenient technology, the so-called "replica trick" for 
treating field theory in the presence of a quenched ran- 
dom field a{x). Thus we consider a general heterogeneous 
system described by a Hamiltonian H[u^ a{x)], where 
u — u{x) is the physical degree of freedom, while (j{x) 
denotes all quenched random parameters. We are inter- 
ested in the quenched average of the free energy 



(El) 



(}F ^ ~\ogZ = -log J Due-f^"^^^^^ 
and of other physical quantities (correlation functions) 

(E2) 



{0[u]) 



DuO[u]e-l^"i^-^'i, 



as average representation of a heterogeneous sample char- 
acterized by a{x). Above, we use (O) for the average over 
thermal fluctuations (functional integral over u) , and use 
O to denote the average over disorder realizations a. The 
main difficulty in the calculation of Eq. (|E1|) and Eq. (jE2p 
is due to the log and inside the disorder average. 

To overcome these obstacles, we deflne a n-replicated 
Hamiltonian 

as well as the associated partition function Z„ and the 
free energy F„ by: 





n 






a=l 


(E3) 


Zn = 


a = l 


(E4) 




-T logZ„, 


(E5) 



where n is an integer. Even though the replicated theory 
is only well-defined for integer n, we analytically continue 
it to a real number < n < 1 so that F„ becomes a func- 
tion of a continuous variable n. The key advantage of this 
replica method is that the resulting replicated Hamilto- 
nian Hn[u°'] is disorder- free and can therefore be studied 
using standard tools for treatment of a homogeneous field 
theory. Using the identity 



T log Z 



T 

lim — log , 

rwO n 



(E6) 



we easily see that the disorder averaged free energy F is 
related to the replicated free energy Fn by 

F^\vai-Fn. (E7) 
n— ►o n 

A potential problem of the replica technique is the break- 
down of commutability of the thermodynamic limit with 



the replica n ^ limit, on which identity Eq. (jE6[) is 
based. It can be shown that for problems (like the one at 
hand), where no real glass physics emerges (as in a spin 
glass problem) , replica trick is innocuous and is simply a 
convenient way to throw out unphysical diagrams 90]. 

The replica technique can also be used to calculate the 
disorder- averaged correlation functions, such as Eq. (jE2[) . 
To see this, we let (O") = (0[w"]) be the "thermal" 
average of a replicated analogue of an observable O in a 
replicated theory. We can show that it is related to (O) , 
the same physical quantity averaged over both thermal 
fluctuations and quenched disorders, by 



J 



-Hn/T 



where H'^ 



= (Z")-iZ" {0[u]) {0[u]), as n 0(E8) 
H[u", a], and in the last step we have used 

Z^=l + nlogZ+ >1. (E9) 

Using similar technique we find 

jAj, (ElOa) 



{A){B), (ElOb) 
(A"'B^ ■■■C) = {A){B)---{C), (ElOc) 

where a, /3 • • • , 7 are assumed to be all different. 

The correlator matrix G"^ of the replicated theory, 
defined by 



(Ell) 



therefore gives 



G' 



■a(3 
ah 



{Ua){ub) 



(E12) 



for a ^ (3. It can be expressed as a sum of two parts: 



G 



a(3 
ah 



Gah^ap + Gq^, 



where 



Gah {UaUb) 



{Ua){ub) 



[Ua - {Ua))iub - (Ub)) 



(E13) 



(E14) 



is a thermal correlator, characterizing thermal fluctu- 
ations around the random ground state for a given 
quenched disorder realization, while 



{Ua){ub) 



(E15) 



gives the quenched correlator, characterizing the sample- 
to-sample fluctuations of the ground states. To the lead- 
ing order, it can be shown that G^j^ is linear in the dis- 
order variance A but is independent of temperature T, 
while G^^jj is linear in temperature. 
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APPENDIX F: DERIVATION OF RG FLOW EQUATIONS FOR HETEROGENEOUS ELASTOMER 

In this appendix we present details of the rephcated RG analysis to treat nonlincarities and random stresses of a 
heterogeneous elastomer. The procedure we follow is quite similar to that in Appendix O for a homogeneous system, 
but now applied to a replicated heterogeneous elastomer Hamiltonian, (|6.20p . To this end, we need to calculate the 
cumulant expansion Eq. (IC15|) for the replicated elastic Hamiltonian Eq. (|6.20p . The harmonic part of the Hamiltonian 
is given by Eq. (I6.22p . with the temperature rescaled away, and set to zero. The nonlincarities are given by 

n 

Hi[u"] = (^y ^^u'^^,u'^ + B^Jkl a,<a,<afe<aiM?) , (fi) 

a = l 

where [u] and Bijki are given in Eqs. (jClSp . We may represent these nonlincarities by the Feynman diagrams 
shown in Fig. 1211 with the understanding that the replica index a is always summed over. 



Hi[u'...u"] = 




FIG. 21: Feynman diagrams for cubic and quartic nonlinearity in the replicated elastic Hamiltonian Eq. (|6.20p . The replica 
index a is summed over. Solid lines represent while the wiggly line represents ylijlii"]. 



As usual, we shall use an internal line (straight or wiggly) for a high-wavevector harmonic propagator of the Uz 
phonon field G"f . From Eq. (|6.33ap . we see that this propagator is a sum of two terms. The first term is diagonal in 
replica indices describing thermal fluctuations. The second term is independent of replica indices and is proportional 
to the disorder variance A, describing quenched fluctuations. We can graphically represent this harmonic propagator 
(a solid line) by a sum of two lines, corresponding to thermal (T) and disordered (A) contributions, illustrated in 
Fig. El 

^ T ^ A 



FIG. 22: The replicated propagator is a sum of a thermal term (marked with T) and a quenched term (marked with A). 

As in the thermal case, we only need to keep track of Feynman diagrams that renormalize quadratic terms in the 
elastic Hamiltonian Eq. (|6.2ip : underlying rotational symmetry (Ward identities) guarantee that nonlincarities have 
identical renormalization, so as to preserve the form of the nonlinear strain tensor w. We start with a Feynman 
diagram on the left hand side of Fig. I23i which renormalize elastic moduli -B^, Az_l, A and /i. The corresponding 
correction to the elastic Hamiltonian is therefore given by 

= -J2 ^'^'(2*^)^ {d - l){d + 1) ^^''^''^ ^ ^'''^'^ ^ ^'^^'''^ [^'']Aki [u'^] X 

/oo 
dp, [5cp{Gzz{p) + 2 ApiC(p)3) ^ A2p4 g^^(p)4] _ (P2) 
-oo 

Using the decomposition shown in Fig. 1221 tl^G Feynman diagram in the left side of Fig. [23] can be separated into 
three pieces, each proportional to A°, A^ and A^ respectively. This is illustrated on the right hand side of Fig. [531 K 
is clear that the first two terms (first two Feynman diagrams on the right hand side of Fig|23l) renormalize the elastic 
constants B^, A^j^, A and ^, while the last piece generates irrelevant contributions that will therefore be neglected. 
The first piece is independent of the disorder variance A and therefore describes thermal fluctuations. According 
to our previous analysis, it is less relevant than the second disorder-dependent piece and can therefore be neglected. 
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FIG. 23: Feynman diagram renormalize quadratic couplings of the elastomer model. Each internal line represent a propa- 
gator, which is itself a sum of two terms ( T + A). The first two diagrams on the r.h.s. renormalize Bz, Azx, A and fj,. The 
last diagram generates terms only irrelevant terms that are therefore neglected. The first term in r.h.s. comes from thermal 
fiuctuations and therefore is also irrelevant. 




FIG. 24: Feynman diagram that renormalize disorder variance A and splay constant K. 



We thus focus on the second piece in Fig. [531 which corresponds to the second term in the integrand in Eq. (jF2p . 
Consequently Eq. (|F2|) reduces to 

+ 2{3X^-6\,±X + 3piX + 3Xl^ + ^l^-3X^±^l){^^u'^f+^l\^^uf)^]. (F3) 

where fi is defined in Eq. (j5.12p and we have set d = 5 in subsequent calculations. Interpreting Eq. (jF3p as a 
coarse-graining correction of the harmonic part of the Hamiltonian, i.e., equating it to 

^ Ki?.(a,<)2 + 2<5,A,^(a,zi^)(9,<) + {dgX + Sgfimu^)^ + Sgfiid^u'^)'] . (F4) 

we obtain the diagrammatic corrections to elastic moduli B^, Xz±, X and /i to be: 

S,B. - -i^.Sl'-^^j=^, (F5) 
8 y^K^ n 

AX /A/ 3A(i?.-A,^)(2A + A^-2A,^) 

dgX,j_ = ^PdOl , (F6) 

16 \/ K'^ \i 

Feynman diagrams renormalizing A and K are shown in Fig. 1241 As before, we can expand each diagram in power 
of disorder variance A. It is not difficult to see that the part linear in A renormalizes the splay constant K while the 
part proportional to A^ renormalizes A itself. 
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The corresponding analytic expression (in momentum space) is given by 
-i ■2-2Y,{<l^<m{<lk<{-q)) r 

+ {A,[u'^ip}]piu.{~p))>{ip, + + q)Aki[u^{p + q)])>) 

= (^9^ ^i'^"/' + + irrelevant, (F9) 

where, as illustrated in Fig |24i we have used p to label the large momenta that run over the momentum shell 
Q e^'' < p < Q, while q is the small momentum below the shell. 

We use Mathematica to calculate the left hand side of this involved expression and expand it in powers of the 
external momentum q. As is indicated on the right hand side of the same equation, we only need to keep track of 
two types of terms: terms that are diagonal in replica indices and proportional to q'^, thereby renormalizing K, and 
terms that are are independent of replica indices and proportional to q^ , thereby renormalizing the disorder variance 
A. Furthermore, we keep only corrections to K and A from quenched fluctuations, as thermal fluctuations are less 
relevant. This ensures that SgK/K and SgA/A are proportional to the disorder variance A. The final results are 
given by 



We follow this coarse-graining with the rescaling transformation R^^{e^^^)R±{e^'') to restore the uv momentum 
cutoff back to Q. The infinitesimal rescaling leads to the following corrections to the coupling constants: 

{6rB,,6rXz±,SrX,drfi) - {d + 3 - 3 Lu){B,, X,±, X, fi)6l, (F12a) 
SrK = {d-l-Lu)KSl, (F12b) 
SrA = {d+l-Lu)ASl. (F12c) 

Putting together the graphical and rescaling corrections we find the one-loop (lowest order in e) RG fiow equations 
of elastic constants and disorder variance A: 

— — = {d + 3-3uj)B^-ijd ^ r—^ , (F13a) 

d I 8 V ^ 

dX.i_ , , 3A(i?,-A,^)(2A + /i-2A,^) 

= (d + 3-3a;)A.^-l-^. , (F13b) 

dX , A(6(A-A,^)^ + 6m(A~A,^)+m^) 

- = (d + 3-3a;)A-V^, , (F13c) 

^ = (rf + 3-3.),-^.^^^, (F13d) 
'-^ ^ id-l-.)K + ^MB.iX + 2,)+20,iX + ,)-X..iX..+4,)) 

^ ^ + (F13f) 

dl 32Vk^ 



where here we defined 
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We then define four dimensionless coupling constants gL, g_L, x and y by Eqs. (|6.34p . Their flow equations can be 
obtained fi-om Eqs. (|F13p . facihtated by Mathematica to be: 

dgh gL eL(gL,g_L,a;,?;) 

■dT - ^^^-16>.(g.,g.,.,,)' 

-df = ^^^-T6-*,(g„g,,x,,)' 
^ = -IgLvii-y') (F15d) 

where e = 5 — d and 

eL(gL, g±, y) = -2500gi^y* + 6000giV^y3 + 63500g£gi x/^j/^ _ lOOOOgi^y^ 

-115000g£giy2 _ 600g3xy2 _ 207250gLgia;2;2 - STSOOg^g^xy^ 
+150000gia:3/2y + 9600gLgia:^/2y - 4320g£gi.T3/2y - GOOOg^V^y 
+176250gLgi^/iy - 36500g£g± + 12500gi^ + 159375gLgi 
+60000gia;2 + 24000gLgia:2 + 576g^g_La;2 + 137500g^g_L 

+600g^2: + 93750gia; + 220750gLgi2; + 38400g^g_La;, (F16) 
^LigL,g±,x,y) = lOOOglV^y^ ~ 2500gly^ ~ lOOglxy^ - 6150gLg±xy^ 

+9000gix3/2y - 720gLg±x3/2y - lOOOg^v^y + 4500gLg±V^y 
+2500g£ + 3600gix2 ^ 96gLg^x2 + 3750gLg± 

+100g^2; + 5625gi2; + 6300gLg±2;, (F17) 
eAgL,g±,x, y) = -2500g3y4 + 64500g2g^ Vij/^ + 5000gi^y2 _ nysoog^g^y^ 

-213400gLgia:y2 - 700g£gia;2/2 + 159000gia;3/2y - 45120gLgi.T=^/'y 
+180750gLgiV^y - 64500g^g± V^y - 2500gi^ + 163125gLgi 
+63600gi.T2 + 2496gLgi.T2 + 117500g£g_L + 99375gia; 

+193300gLgi.x + 700g£gix, (F18) 
*±(gL, g±, 2:, y) - lOOOg^V^y^ - 2500g^y2 _ loog^^y^ _ 6150gLg±a;y^ 

+9000gix3/2y _ 720gLg±.T3/2y - lOOOg^V^y + 4500gLg±%/^y 
+2500g£ + 3600gia;2 + 96gLgix2 + 3750gLgj. 

+100g^x + 5625gix + 6300gLgj.2;. (F19) 

As noted earher the bare values of x and y for typical elastomers are much smaller than unity. Furthermore, below 
five dimension, we expect gL(0 to fiow to a finite positive value gL- These considerations, together with Eq. (|F15cp 
and Eq. (|F15dp indicate that x{l) and yil) start from small values and flow to zero exponentially according to: 

x{l) w xe^^sIVa^ (F20a) 
y{l) w 2/e~3giV4_ (F20b) 

Consequently we can set them to zero in the flow equations for gL and gj^ . This leads to a considerable simplification 
of Eq. (|FT5a]) and Eq. (|FT5a)) : 

dgi. 5gL (4g^ + 44g^gL + Slgj) 

d I 64gL + 96g_L 

dgi_ g± (-4g! + 188g^gL + 261gi) 

-^ = ''- 64gL + 96g, ■ (^2^^) 

Using Eqs. (|2.46p . Eqs. (|F13p . and Eqs. (|6.34p . we can derive fiow equations for elastic constants B, C, /il, and 
fi. The results are shown in Eqs. (I6.35p . The exponents rjB, rji^, r]j_ are the same as in the first three equations of 
Eqs. (|6.36p . which we duplicate below: 

3 3 1 

'7B = g 5L, VL = gL, V± = Yq 9±- (F22a) 



59 



r]K and ?7a are more complicated but can be straightforwardly found with the crutch of Mathematica to help with 
the algebra: 

^ 5 (^10 {y^ -l)gl + g± (24x - SOyyS + 175) gL + lOOgjx) 

16X2(75gia; + 2gL(.T-10tjV^+25)) ' ^ ' 

^ gL (3g± (16x + 40yVi + 25) - 50gL {y^ - l)) 
''^ 75grx + 2gL (x - lOy^ + 25) ' ^ ^ 

In the limit x and y approaching zero, these reduce to the l ast two equations in Eqs. (|6.36|). 
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